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SUMMAR\ 


A general outline is given of possibilities for favorable Ve 
interference between supersonic aircraft components rhe 
interference between a three-dimensional wing with swept lead 
ing edges and a small reflecting surface in biplane arrangement 
is discussed in terms of linearized flow theory Relationships 
are developed which facilitate the determination of the inter 
ference drag of a large class of such three-dimensional configu 
rations without the need for determination of the flow field. The 
method is applied to several biplane configurations of practical 
interest. It is shown that the pressure drag of the wing can be 
greatly reduced or that the volume near the root section of the 
wing can be considerably increased without penalty in pressur« 
drag. The three-dimensional biplane is then discussed from the 


practical point of view 
SYMBOLS 


area of integration 
1,7 ])'/2 
Cy = drag coefficient based on plan-form area of wing 
= interference drag coefficient based on plan-form area of 
wing 
= pressure coefficient 
= chord at root section 
= drag 
= interference drag 
= dimension shown in Fig. 9; see Eq. (80 
= quantity proportional to intensity of source distribution 
= distance between planes of wing and plate 
tangent of sweepback angle of associated edge 
= Mach Number 
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Mach Number of undisturbed streat 

k/B 

dimensionless distance in Fig. 10 

dimensionless distance in Figs. 14 and 18 

radial coordinate in cylindrical coordinate 

area of integration 

area of reflecting plate plan form 

area of wing plan form 

distance shown in Fig. 9; see Eq 

length of ar¢ 

ratio defined by Eq. (47 

x component of perturbation velocity 

velocity of undisturbed stream 

ry component of perturbation velocity 

> component of perturbation velocity 

streamwise Cartesian coordinate 

lateral Cartesian coordinate 

normal Cartesian coordinate 

oblique coordinate used in Appendix; see Fig. 8(a 

oblique coordinate used in Appendix; see Fig. 8 

variable used in Appendix 

normal Cartesian coordinat: 

lateral Cartesian coordinate 

cylindrical angular coordinate 

inclination of surface with respect to plane of wing In 
plane of constant y 

Mach angle 

streamwise Cartesian coordinat« 

density of undisturbed stream 

maximum thickness at root section relative 


perturbation velocity potential 


INTRODUCTION 


8 be: rOTAL DRAG produced by the motion of a body 
at supersonic speeds is considered to be comp sed 
of (1) the friction and form drag directly related to 
viscous effects near the surface of the body, (2) the vor- 
tex drag due to the lift distribution, and (3) the wave 
drag due to the lift distribution and to the distribution 
of the thickness or volume of the body. In the non- 
viscous theory of supersonic flow it is the last two that 
appear as a variation of momentum at infinity. All 
types of drag are sensibly affected by the interference 
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effects resulting from the presence of more than one 
body in the stream, and it has been known for some 
time that wave drag can be notably reduced by judicious 
use of interference effects. It was shown in 1935 that 
the wave drag due to thickness of a two-dimensional 
profile can be eliminated by the introduction of either 
another wing or a reflecting surface, these being de- 
signed to establish a wave pattern in which no waves 
are propagated away from the system.’ The first 
configuration is usually referred to as the Busemann 
biplane. 

The purpose of this report is, first, to outline different 
types of favorable wave interference in three-dimen- 
sional systems and, second, to discuss in some detail 
the use of interference of the biplane type for the re- 
duction of wave drag due to thickness of three-dimen- 


sional swept wings. 
FAVORABLE WAVE INTERFERENCE IN SUPERSONIC FLOW 


In supersonic flow, shock waves are produced either 
as the envelope of compression waves or by discon- 
tinuities in the shape of stream surfaces imposed by the 
presence of the body in the flow field. The presence of 
shock waves produces an increase in the entropy of the 
stream and, therefore, a loss in momentum correspond- 
ing to wave drag. In the first approximation the loss of 
momentum per unit area through a given stream tube 
extending from minus infinity to plus infinity is pro- 
portional to the product of the entropy increase and the 
free-stream pressure. When a body travels at super- 
sonic speeds both compression and expansion waves 
are produced by the motion. The compression waves 
are normally cancelled by the expansion waves at large 
distances from the body but, in the general case, the 
compression waves will form an envelope before being 
cancelled by the expansion waves and thereby produce 
drag. In order to reduce entropy losses the cancella- 
tion of shock or compression waves by means of either 
interfering expansion waves or interfering surfaces 
must be accomplished as close as possible to the body. 

If the interference is analyzed in terms of small per- 
turbation theory, these nonlinear features of the wave 
pattern disappear, and the momentum losses appear 
differently. Let the system be enclosed in a control 
surface consisting of a right circular cylinder whose 
axis is parallel to the direction of motion; then let the 
lateral curved surface and the downstream end face 
recede to infinity (see reference 2 for proper limiting 
procedures). By application of the momentum theo- 
rem to this control surface, the following expression is 
obtained for the wave drag: 


Dp = | uv,rddodx (1) 
« 7 S 


The integration is extended over the infinite lateral 
surface S of the cylinder. It is also found that the 
cylindrical velocity perturbation components are mu- 
tually related by the expression? 


u = —(v,/B) 
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on the cylinder. No contribution to the wave drag 
appears on the downstream end face of the control sur- 
face. In the case of lifting systems, however, a con- 
tribution to the drag from this surface does arise which 
is, in fact, the vortex drag of the system. It is re 
marked that the problem of minimizing the wave drag 
of a system is not, in general, compatible with the 
problem of minimizing its vortex drag, so that the sum 
of the two forms of drag is the quantity to be subjected 
to minimization in a general treatment of arbitrary 
lifting systems. 

In terms of the linearized theory, favorable wave 
interference in a given system is obtained when meas 
ures are adopted which decrease the magnitudes of the 
disturbances reaching the remote lateral curved sur- 
face S, thereby diminishing the value of the integral in 
Eq. (1). Two general methods for accomplishing this 
end present themselves. These are discussed below. 

The type of favorable 
achieved as follows. An 
which is induced by a modification of the primary sur- 


interference is 
field, 


first wave 


interference flow 
face or by a suitably chosen complementary surface 
in the system, is superposed upon the flow field induced 
by the primary surface of the system. The second 
surface is designed so that the primary induction reach- 
ing the remote control surface is attenuated by the 
interference induction. Then deviations of the pres- 
sure from the free-stream value throughout the flow 
field are, on the average, reduced; variations in area of 
individual stream tubes are correspondingly suppressed. 
Because no one-to-one correspondence exists between 
points on the surface of the system and points on the 
remote control surface as in two-dimensional flow, care 
must be exercised that attenuative overlapping of pri- 
mary and interference induction over one part of the 
control surface is not offset by reinforcement of the 
two flow fields over a second part of the control surface. 
That is, the induction due to the primary system at a 
given point on the remote surface can be reduced, say 
to zero, by producing a suitable interference disturb- 
ance somewhere on the surface of the system. No 
assurance exists, however, that this interfering flow 
field will combine with the primary flow field to produce 
similar attenuation of induction at other points on the 
control surface. These favorably interfering flow fields 
are established by designing, positioning, and contour 
ing the several components which constitute the air- 
craft—e.g., wing, body, tail, and inlet—-such that they 
are complementary, the flow fields of each dovetailing 
to produce minimal total induction on the remote 
cylinder. These geometric modifications are limited 
by the geometric, aerodynamic, and other practical re- 
straints imposed on the system. The total volume and 
lift of the system might well be prescribed; the struc 
tural integrity must be maintained; important non- 
linear and viscous effects attendant upon undulatory 
fuselage contouring, for example, must be considered. 
Restraints such as prescribed volume and lift are im- 
portant to a formal analysis in order to rule out, for 
example, the trivial result that the wave drag is a 
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minimum when there is superposed upon the primary 
system its negative. Itisapparent that the entire con- 
figuration must be designed as an entity if the possi 
bility of important drag reductions by favorable wave 
interference is to be fully realized. 

This favorable interference is perhaps more con- 
veniently envisaged in terms of a direct determination 
of the drag by pressure integration over the surface of 
the system. (Strictly speaking, the pressure integra- 
tion gives the sum of wave drag and vortex drag. 
Then favorable interference between two components 
is established when the flow field induced by each im- 
It is clear that this inter- 
ference is a the 
position, and of the contours of the two components. 


parts thrust to the other. 
function of the design, of relative 
The drag of the two components is given by the sum 
of the integrals over each surface of the self-induced 
pressure and of the interference pressure. It is quite 
possible for the sum of these four terms to be less than 
the self-induced drag of either component. 

A second type of favorable interference is also pos 
sible. Disturbances originating from the primary 
surfaces are intercepted by interference surfaces before 
they arrive at the remote control surface. These dis- 
turbances may be cancelled on the interference surface 
or reflected back to the primary surface where they are 
No induction at infinity due to these 
disturbances then appears. With interference of this 


type, complete elimination of the wave drag due to 


then cancelled. 


thickness is possible in principle. In terms of pressure 
distributions over the system, the effect of this type of 
interference is to balance the high pressure distributed 
over the upstream projected area of the system with 
high pressure distributed over the downstream projected 
area. 

Drag reduction by means of the first type of favor 
able interference recently has received considerable 
attention. Stimulated by the implications of the tran 
sonic area rule* and utilizing the concept of the equiva 
lent source position,* a number of investigators 
have made systematic mathematical studies of favor 
able wave interference at supersonic speeds by con 
sidering the effectiveness at the remote control surface 
distributions of disturbances 


of planar and_ linear 


representing wing-body combinations. Fuselage con- 
touring is proposed which attenuates wing-body dis- 
turbances at infinity. Wing-body combinations which 
exhibit minimum drag under prescribed restraints are 
discussed. Experiments indicate that the significant 
drag reductions predicted by the theory are often 
possible, although the theoretical findings are subject 
to considerable practical interpretation. 

Other investigators have observed that considerable 
favorable interference of this type between aircraft 
components can be realized by considering directly the 
nature of the pressure field induced by each component 
on adjacent solid boundaries. Attempts are made 
to orient and contour the surfaces of components ap- 
pearing in a given system so that each component re- 
ceives a thrust contribution from the pressure fields 
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other These engineering 


which are predicated on the direct 


induced by components. 
considerations, 
rather than the indirect formulation of drag, provide 
a relatively simple approach to wave drag reduction in 
but from 


which absolute optima are not usually forthcoming. 


which practical restraints are prominent 


As an elementary example, consider a protuberance 
consisting of the triangular pyramid ABFC which is 
followed by a semi-infinite right triangular prism (see 
Fig. 1). The protuberance rests on the infinite x-y 
plane in an axial supersonic stream as shown in Fig. |. 
The protuberance might represent, for example, an 
mounted on a circular cylindrical 


aircraft canopy 


fuselage of very large radius. It is proposed to con 
tour the x-y plane in the neighborhood of the canopy in 
order to establish a favorable interference. Suppose 
that a recess or indentation is produced in the plane 
whose shape is again described by a triangular pyramid 
semi-infinite triangular prism. 


followed right 


The canopy then rests in the recess as shown in Fig. 1. 


by a 


If the surface inclinations are regarded to be sufficiently 


winglike, the geometry is represented in terms of 


linearized wing theory (see for example, reference S 
by a uniform distribution of sources over ABC and a 
over ABD and ACE. 


above the x-y plane is 


uniform distribution of sinks 
If the height of the point F 
7c and if the depth of the point G below the x-y plane 
is 7c, then the inclination in the planes parallel to the 
v-z plane of surfaces ABF and ACF is 72 and the like 
ABD and ACE is The 


surfaces downstream of the line x c exhibit no such 


inclination of surfaces — Ts 
inclination and thus require no source distributions. 
The sweep of the leading edges of the canopy is given 
k, B> 1 and of the recess by n, k, B. The 
drag of the system, supported on the exposed surfaces 
It is evi 


by m2 


of the pyramids only, is readily calculated. 
dent that the pressure field induced by the exposed 
recess is of opposite sign to that induced by the canopy 
and that a favorable interference has been established 
On the other hand, the selfinduced 
With a suitably chosen 


between the two. 
drag of the recess is positive. 
recess, the drag of the system can be reduced to a value 
below that of the isolated canopy. 

From the computed drag of the system the optimum 
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0 2 3 4 5 6 
n, 
Fic. 2. Relative drag reduction of canopy configuration for 


optimum inclination of indentation; c, m2, and 72 prescribed 


values of the recess parameters m and 7 have been 
computed for two cases. First regard the values of 
N», Tz, and ¢ as prescribed; then the height of the point 
F on the canopy above the x-y plane is fixed. The 
minimum drag of the system for optimum value of 7; 
may be obtained in terms of m by differentiation. In 
Fig. 2 is presented the relative drag reduction as a 
function of , for various values of the canopy sweep 
parameter #:. The quantity Cp represents the drag 
of the isolated canopy for the same values of m2, 72, and 
c. The quantity AC» is defined as the drag of the iso- 
lated canopy minus the minimum drag of the system. 
The corresponding optimum ratio 7;/72 is also plotted 
versus , for various values of m. in the same Figure. 
The optimum value of », for the prescribed m. may be 
noted from the stationary point on each curve for 
ACp/Cp. This value of 7 is only slightly less than the 
corresponding value of m. It is seen that the maxi- 
mum value of ACp Cp increases with increasing 7» 
whereas the associated value of 7;/72 decreases with 
The maximum value of ACp,Cp is 
When the recess is optimum 


increasing Me. 
about !/5 when nm. = 6. 
for a given canopy, it is seen from Fig. 2 that 7; is 
several times 7:. But, because , is only slightly less 
than m2, most of the recess is covered by the canopy and 
the actual indentation depth is much less than 7¢. 
Suppose now that mm, (7; + 72), and ¢ are prescribed; 
in this case the total height of the canopy (71 + 72)¢ is 
fixed (this is the length FG in Fig. 1) but the canopy is 
permitted to sink into the recess without constraint. 
The calculations just described for the first case may 
be repeated for this case. In Fig. 3 is presented ACp 
Cp versus , for several values of m). and the optimum 
value of 7,/(7; + 72) versus m, for the same values of 
m.. For given m and for m < m2, it may be seen that 
the drag reductions tend to be large and that the pro- 


tuberance extends above the x-y plane. When », = 
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M2, the drag of the system is zero and 7 = (7; 
We thus arrive at the trivial result wherein the pro- 
tuberance is completely submerged in the recess and 
the x-y plane is completely free of disturbance. This 
simple problem illustrates on the one hand that large 
drag reductions can be obtained by use of favorable 
interference and on the other hand that care must be 
exercised to impose realistic restraints on geometric 
modifications. 

A classical example of the second type of inter 
ference is the two-dimensional Busemann_ biplanc 
The configuration, consisting of two symmetrically 
arranged wings, 1s illustrated in Fig. + in the nonlinear 
instance and in Fig. 5 in the linear instance. The 
waves produced at the leading edge of each wing are 
reflected from the opposite wing, as shown in the 
figures, and are then cancelled on the back of the 
original wing. The wave drag of the system (to the 
second approximation) is zero because no waves are 
propagated to infinity. The pressure distribution on 
each wing is symmetric. An axially symmetric coun- 
terpart, consisting of a body enveloped by a cowling 
(see Fig. 6), was proposed in reference 9 and analyzed 
in reference 10. 

It is of interest for later application to discuss further 
the nature of the linearized flow field in Fig. 5. The 
disturbances produced by the surface AB are trans- 
mitted unchanged to the flow in the region between 
parallel waves AE and BF. If the surface EF has the 
same inclination as AB, no reflections of the waves AE 
and BF occur. The pressure forces on the two sur- 
faces are equal and opposite so that their resultant is a 
pure couple. The surfaces DE and BC bear a similar 
relationship, supporting a couple equal and opposite 
to the first. It is evident that the flow field in the re- 
gion ABCNTA is unaltered if the stream surface TN 
is replaced by a flat reflecting plate of zero thickness 
and if the upper wing is subsequently removed. Then 
the waves AT and BN, produced by the surface AB, 
are reflected from the plate TN as waves of the opposite 


n=2 3 o C 
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Fic. 3. Relative drag reduction of canopy configuration for 
optimum inclination of indentation; c, m2, and (7; + 72) prescribed 
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Fic. 5. Linear Busemann biplane 


family. There is produced along the surface BC a 
pressure disturbance of the same sign and a deviation 
of opposite sign as that existing along surface AB. 
Since the inclination of BC is equal and opposite to 
that of AB, no reflections of the waves TB and NC 
occur. The axial pressure forces on the surfaces AB 
and BC balance each other and their resultant normal 
force balances the resultant normal force supported by 
the plate. 

The remainder of the present report is concerned 
with analogous favorable interference between three 
dimensional swept wings and reflecting surfaces in bi 
plane arrangement. Reflecting surfaces will be pro- 
posed which reduce the wave drag due to thickness of 
the given wing. The necessity of dealing with the re 
flection of three-dimensional waves makes an analysis 
of the problem considerably more involved than the 
analogous two-dimensional treatment. 


DETERMINATION OF THE INTERFERENCE EFFECTS 
BETWEEN THREE-DIMENSIONAL SURFACES 


The interference between a wing and a reflecting sur- 
face were analyzed in reference | for the case of two- 
dimensional supersonic flow. In reference 11 two iden- 
tical wings of rectangular plan form oriented in the 
Busemann relationship were analyzed by linear theory 
to determine the tip effects due to finite aspect ratio. 
The more general case of a three-dimensional wing of 
variable chord interacting with a reflecting surface 
can also be analyzed by means of linearized wing 
theory. 

Consider, for example, the biplane configuration 
shown in Fig. 7. Let the geometric angle of attack of 
both wing and reflecting surface be zero. Suppose for 
simplicity of discussion that all edges of each plan form 
are supersonic and that a source-sink distribution, of 
intensity proportional to the local slope Oz Ox of the 
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to each considered 


isolated, in the usual manner. 


surface, is used represent wing, 
The pertinent velocity 
and pressure field may then be calculated from the 


well-known relations 


g(x, ¥, 2 
({ (E, n)d&dn 
») 
FIV = 3 B2(y — 9 Biz —¢ 
Ux, ¥, 2 Oy Or 3 
W(X, V, 2 Og Oz { 
C; 2u/V; ») 


where the integration is to be extended over the planar 
regions 7 Which appear inside the fore Mach cone from 
the point P(x, y, 2). 

The distribution of w in the plane of the reflecting 
surface due to the wing may be found from these rela 
tions. An additional distribution of sources must be 
added in the region of the reflecting surface to cancel 
this upwash and thereby to maintain the boundary 
shape, the source intensity at each point being propor 
tional to the value of w at the same point. The pres 
sure distribution over the reflecting surface, due to the 
source distributions in this surface and in the forewing, 
may then be calculated. Similarly, the distribution 
of downwash over the back of the wing due to the com- 
plete source distribution in the reflecting surface is 
to be computed. An additional source distribution must 
be added on the upper wing surface to cancel this down 
wash. The pressure distribution over the wing may 
then be computed, and can be considered in two parts: 
the pressure distribution induced by the sources repre 


senting the isolated wing, and the pressure distribution 























Wing-reflecting surface configuration 
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Fic. 8(a). Flow field of elemental source distributions 


induced by both the complete distribution of sources 
representing the reflecting surface and the interference 
sources in the afterwing. The first leads to the drag 
of the isolated wing and the second to the interference 
drag of the wing. If the wing lies in a region influenced 
by the flow behind the trailing edge of the reflecting 
surface, then the influence of this wake on the wing 
must also be considered. 

This tedious procedure is unlikely to suggest the eri 
teria by means of which superior configurations—4.e., 
those giving substantial favorable interference, might 
be proposed. It would appear that integral relation- 
ships, providing information on the total amount of 
upwash and on the total pressure force without recourse 
to determination of the flow field, would implement 
the development and analysis of suitable wing-reflecting 
surface configurations. 

Such relations are in fact suggested by an exami 
nation of the linearized flow field of the two-dimensional 
Busemann biplane which is illustrated in Fig. 5. As 
previously discussed, the drag of the system is zero be- 
cause the disturbances produced by the surface DE are 
transmitted to the surface BC, along which there re- 
sults a balancing thrust and, collaterally, the closure of 
the lower wing without reflection (surfaces AB and 
EF have a similar relationship). For generality we 
superpose a spanwise component to the undisturbed 
stream shown in Fig. 5; according to simple sweep 
theory the spanwise velocity remains everywhere the 
same and the perturbation flow field is unchanged. 
Let k be the tangent of the resulting sweepback angle. 
We then consider further the relationship between the 
complementary surfaces DE and BC with the aid of 
Fig. S(a). The surface a is the half of the chord plane 
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DF in Fig. 5 which supports the inclination DE, and 
the surface b is the half of the chord plane AC in Fig. 5 
The streamwise 


chord of surfaces a and b is denoted by c in the Figure. 
for the pur 


which supports the inclination BC. 


Surfaces a and b are supplied in Fig. S(a 
poses of reference and orientation. The inclination DE 
in Fig. 5 is represented in Fig. S(a) by distributing over 
the lower face of the surface a a source distribution 
of constant intensity proportional to /. 

Identify on b a streamwise element of area EF whose 
chord is c and whose elemental width is dy. For con- 
venience the leading and trailing edges of this element 
(and of similar elements to follow) are taken parallel 
to the y axis; the discrepancies are of higher order. 
The fore Mach cones from points E and F result in the 
hyperbolic traces TAU and SBW, respectively, in the 
plane z = h. These hyperbolas are tangent to the lead 
ing and trailing edges of plane a at points A and B, 
respectively. The perturbation velocity field estab- 
lished over the element EF (and, indeed, everywhere in 
the region between the plane of the leading edges of a 
and b and the plane of the trailing edges of a and b) by 
the constant inclination existing over the lower surface 
of a follows immediately from two-dimensional flow 


considerations as 


u —af/BV 1 — n’ 6) 

w= —nof (7 

in terms of the wind axes shown. Here 1 is defined 

as the ratio of kto B. These perturbations are induced 

by disturbances appearing in the region GBH of the 

plane a. We can then form the two area integrals 
over the elemental area /F 


Ucpu dl = —afcdn BWV | ne = 
J EF 


wopu del 
EF 


—nfcdyn 

Let the region of plane a between points G and H 
be subdivided into a lateral array of streamwise ele 
mental areas such as AB and A’B’, each having the 
uniform 


Each supports a 


The induction over EF, given by 


elemental width dn. area 
source distribution. 
Eqs. (6) and (7), can be regarded as the summation of 
the contributions of each of these elemental distribu 
tions, since at least part of each influences EF. The 
0 influenced by the points 


regions of the plane s = 
A and B of the elemental source AB are defined by the 
hyberbolic traces OER and PFOQ of their respective 
Mach cones; these are tangent to the leading and trail- 
ing edges of b at points E and F, respectively. The 
element AB influences the plane b in the region CED. 
The element A’B’, whose lateral distance from AB is 
n, produces a similar pair of hyperbolas which are dis- 
placed laterally a distance y and longitudinally a dis- 
tance ky from the first pair. Let E’F’ be an element 
of area displaced from element EF by 
tance (—yn) and a longitudinal distance (—kyn). Be- 


a lateral dis- 


cause the physical phenomena depend on only rela- 
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| the «-y plane from points J and J’, respectively. 


tive distances and not on the absolute coordinates, the 
induction over E’F’ due to AB is the same as the in- 
duction over EF due to A’B’. Suppose, as before, that 
the region of plane b between points C and D is sub- 
divided into a lateral array of streamwise elemental 
areas like EF and E’F’. 
effects, the integral of the induction over all the ele- 


Then, by summation of like 


ments of area which compose the region CED due to 
the elemental source distribution AB is equal to the 
integral of induction over the elemental area EF due 
to the array of elemental source distributions com- 


posing the region GBH. Thus, 


UAR dA = Ucpu dA 10) 
7 CED / EF 
WARB dA = | WGBH dA 11) 
/ CED / EI 


But the right-hand members of Eqs. (10) and (11) are 
given by Eqs. (8) and (9), respectively. We therefore 
obtain the relations 


Uap dl 


7 CED 


~xfcdn/BV1 —n (12) 


Wa pd = —afcdn (13) 
e/ CED 


which give the integral of induction over region CED 
due to AB, 
intensity is given by /, whose chord is given by c¢, 


the elemental source distribution whose 


and 
whose width 1s dy. 

Consider a point L located on the streamwise axis of 
the elemental area EF at an arbitrary distance down 
stream of the point F. Construct a line segment 
through point L parallel to the edges of surface b such 
that it extends from one branch of the hyperbola OER 
from two-dimen- 


L due to 


to the other branch. It is evident 


sional considerations that the induction of 
the above-mentioned lateral array of elemental source 
distributions extending from point G to point H on 
surface a is zero. From the same argument applied 
above, it follows that the line integral of the induction 
AB along the 


due to elemental source distribution 
line segment through L is zero. That ts, 
°N\I 
Uay ds = 0 (14 
JN 
*\I 
Wands = 0 (15) 
JN 


It follows from Eqs. (14) and (15) that the area integral 
of induction due to element AB is zero over any region 
circumscribed by the two branches of hyperbola OER, 
by the line segment CD, and by any line parallel to 
CD which is downstream of CD. 

Consider a point I on the element AB whose distance 
downstream from point A is &. 
point 1’ whose distance downstream from point A is 
{+ dt). Let KJ and K’J’ be line segments parallel 
to CD which are tangent to the hyperbolic traces in 
Apply 


Consider an associated 


NG FL 


OW FIELDS d 


the theorems expressed in Eqs. (12) and (13) to obtain 
the integral of induction over region J‘EK’ due to ele 
ment AI’. Subtract from this the integral of induc 
tion over region JEK due to the element AI. This 
difference gives the integral of induction over the el 
mental area JJ’K’K due to the elemental source dis 
tribution AI’. Noting that the induction due to ek 


ment AI’ is the same as that due to element AB in this 


region, and dividing by the elemental length dé, one 
finds 
*J 
Uap dy tldn BY | n* 16) 
JK 
“J 
WAR dy ahdn 17) 
e K 


It is seen that the line integrals are independent of the 
distance £; thus the value of the integral of induction 
along all oblique line segments like KJ which are within 
region CED is the same. Because Eqs. (16) and (17) 
supply information on the distribution of induction 
due to the element AB in the region CED, they suppl 
ment Eqs. (12) and (13). 

The relations derived provide information on the in 
duction due to a source distribution, whose constant 
intensity is proportional to /, over the elemental area 
cdyn. It is clear that this disturbance can be regarded 
as a line source of constant linear intensity propor 
tional to fdy, provided that the distance h is finite 
Therefore, the relations derived still hold if the element 
AB is rotated about its streamwise axis to an arbitrary 
Moreover, the be made 
The now be 
recast in a readily applicable form with the aid of Fig 
S(b). 


position. distance ¢ may 


arbitrarily small. above relations can 


Consider a rectangular element of area dédw 











on some surface supporting a source distribution. Let 
i 
a Fa f . 
MEDIAN LINE ae ) 
SEGMENT ——_ : 
y / 
dw y 
F ‘ 
y 
x 
ff 








F1G. &8(b Integral relations 








8 [JOURNAL OF THE AERONAUT 
Let 


6u and éw be the x and z components, respectively, of 
the perturbation velocity induced in the x-y plane by 


the lateral edges of the element be streamwise. 


the elementary source distribution / supported by the 


element of area. Then Eqs. (12), (13), (14), and (15 
may be written 
6u dl rf dE dw BY | 1 IS 
JA 
éw dl F arf dé dw 19 
JA 
6u dS 0 20 
Ow ds 0 (?] ) 


where » = k B, where the areas A and § are related 
to the hyperbolas and oblique line segments as shown 
in Fig. S(b), and where the minus sign applies in Eq. 
(19) when ¢ > 0 and the plus sign applies when ¢ < 0. 
An analytical derivation of these integral relations 1s 
given in the Appendix. 

For certain classes of configurations, Eqs. (1S) 
through (21) can be used to determine the total lift or 
drag and the integral of the downwash induced by a sys 
tem of disturbances without computing the flow field 
of these disturbances.* In the present paper they will 
be used for the analysis of the interference between 
the parallel surfaces of three-dimensional biplanes. 
Because only unswept trailing edges will be considered, 
the value of the parameter » will henceforth be zero. 

Consider the evaluation of the drag of the previously 
discussed wing-reflecting surface configuration illus- 
trated in Fig. 7 by means of the integral relations just 
developed. In addition to the aforementioned char- 
acteristics of the configuration, assume that the re- 
flecting surface is a flat plate. Suppose that the pre 
liminary step requiring the representation of the iso 
lated wing by a distribution of sources and the calcu- 
lation of its drag has been accomplished and that the 
Let the distribu- 
tion of sources on the wing which influence the re- 


flecting plate be ABC. The integral of the upwash 


interference drag is to be computed. 


* Those relations above which give the integral of the pressure, 
as well as some additional information on the center of pressure, 
Although 


no applications are presented, the possibility of the application 


are briefly indicated in reference 12 for the case » = O 


of those integral relations to cruciform wing arrangements is 
mentioned the 
author indicates that the relations he gives are generalizations of 
13. The latter 
expressions give the value of the integral of source-induced pres 


Proofs are not presented in reference 12, but 


similar relations which are derived in reference 
sure over the plane which contains the elemental source, as well as 
the corresponding center of pressure. They are used in reference 
13 to compute the lift and center of pressure of planar wings; 
several interesting planar systems with both swept and unswept 
trailing edges are considered. The different point of view pre 
sented herein provides information on the integral of downwash 
as well as pressure and facilitates repeated application of the rela 
tions in multiple reflection problems for the purpose of the calcula 


tion of drag. 


ri 
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NCES JANUAR Y, 
from region ABC over the plan form of the reflecting 
taken with the 
provided that the two 


plate may be computed from Eq. (19), 
21 


hyperbolic traces (which sensibly merge 


plus sign, and from Eq. 
in the plane 
of the reflecting surface from each elemental source in 
region ABC intersect the plan-form boundary of the 
plate only along a straight trailing edge. We assume 
that the leading and trailing edges of the plate in Fig. 7 
conform with this requirement. Thus the shape of 
the plate is given by the curved envelope of hyperbolas 
emanating from points on the leading edge of the re 
gion ABC (coinciding with the leading edge of the 
wing) and by the straight line which comprises the 
envelope of hyperbolas emanating from points on the 
straight line BC. It is further assumed here that the 


trailing edge of the plate (and, consequently, line 


BC) is unswept. The leading edge of the plate is seen 
to extend forward to the undisturbed flow; it may, of 
course, be extended farther upstream with no effect. 
In order to satisfy the boundary conditions, a dis 
tribution of sources is to be added over the plan form 
of the plate which cancels this upwash with downwash, 
By superposition of the effect of each 
(19) 


point by point. 
elemental source in region ABC, given by Eqs. 
and (21), it follows that 


| (ds = —| 


« 


Wlate dS 
plate 


WapcdS = | fdS (22) 
~ plate 7 ABC 


Let the region DEF in the plane of the wing be 
constructed from the envelopes of hyperbolas originat- 
ing from points on the periphery of the plate in the 


manner described above. Let DE’F’ be the entire 
region of the wing influenced by the plate. The in- 


tegral of the downwash induced by the plate on the 


afterwing may be determined as above, providing 


again that the pair of hyperbolas associated with each 
elemental source in the plate intersect the plan-form 
boundary of the wing only along a straight trailing 
Then the 
total source strength in the plane of the wing necessary 


edge, which is here also assumed unswept. 


to cancel this downwash is 


fdS = 
DEF 


/ plate 
If the line EF falls downstream of the trailing edge of 


i ¢@s (23) 


the wing, only that part of the plate influencing the 
afterwing should be considered in the present calcula- 
tion. If the wake of the plate is not to influence the 
wing and if, concurrently, the full effect of the plate on 
the wing is to be obtained, then the plane of the trailing 
edges of plate and wing must be inclined at the angle 
u to the x axis so that EF coincides with E’F’. Such 
an orientation is henceforth assumed. 

The integral over DEF of u due to the plate source 
distribution may be obtained by superposition of the 
effects of each elemental source in the plate given by 


Eqs. (18) and (20). Just as above, 
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USE OF INTERFERING FLOW FIELDS Y 
eflecting ‘ - eee — 
with the | Upiate dS = —(mr B) | fdS (24) | 
~/ DEF 7 plate ¥ 

the two | | 
1e plane In order to find the corresponding integral effect due A y 
ource in | to the sources in the afterwing, Eqs. (18) and (20) | 

of the are used for the degenerate case where the distance ¢ c 

assuine in Fig. S(b) is zero. Thus, | 
in Fig. 7 ~ ° C 
hape of | | Uper dS = —(m B) fds (25) 

J DEF J DEF 
perbolas 
the re- | The value of the integral of u on the wing due to the x*B/ yorh 
> . - . - . . ~ | ' 

of the | complete interference is found by combining Eqs. ' 

ises the 22), (23), (24), and (25) to obtain 
| — 
> ON the ar ‘ -| 
i dj tan 

that the | udS = —2(mr B) fds (26) ; / | \ + . 
ly, line 7 DEF 7 ABC PLANE ah PLATE 
e IS 2een The interference drag is defined by the formula a a 
may, of / \ / | \ 
o effect. Dj/(piV1?/2) = -—(2 rn Au dS (27 

a dis- 7 DEF 


an form 
wnwash, 


If the inclination A is considered constant over region 


qj D 
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: ; |} DEF, Eqs. (26) and (27) then give, for the interference 7 ; 
ot each . ‘ ‘ ‘ ie - 
lrag of the wing, C ff n'k 
qs. (19) ? / fs \ \ ta ’ 5 
. 2,2 
D,/(p:Vi?/2) = (40d Bi) | fdS (28) PLANE 2-0 _| \x-B/y*h latin 
J ABC AF \ E™ 
“— ; / \ / . 
where \ is generally negative over DEF. The drag of 
the reflecting plate, having no inclination, is neces- Fic. 9. Interaction regions when region ABC is triangular 
S (22) | sarily zero. With the specified geometric restrictions, 


wing be 


riginat- 


——————E 


the interference is evidently quasi-two-dimensional; 
itis seen that the interference drag may be readily cal- 


culated. The integral relations also indicate, by in- 


cussed in a preceding paragraph in reference to the two 
dimensional problem. 
Some applications of the method discussed to the 


» in the | spection, that no lift is introduced into these systems a4 a . co ntaid 
; } : determination of the interference drag of wings and 
e entire | when the plate is added. tara : 
a. 3 ; ; i reflecting plates in combination are given in the follow- 
he in- It is possible to generalize somewhat the restriction ? ; TI 6 ; d . 
_ r _ ' ing section. le configurations considered are vari 
-on the | that \ be constant over DEF. With use of Eqs. (18) 8 ais a 8 
ations of Fig. 2. The plates are assumed to have zero 


roviding 
ith each 
an-form 


21), only the lateral distribution of the induction 
from each elemental source is unknown. By applying 


the above procedure to three successive lateral ele- 


thickness; the geometric angles of attack of the wings 
and plates are likewise zero. The size and orientation 


of the plate used with each wing are governed by the 


trailing | mental strips in region ABC, in the plate, and in re- 5 fh, ii : é J 5 a 
7 ; pies : , size of suitable interaction regions ABC and DEF ex- 
hen the | gion DEF—each separated by the distance Bh in the Sod ; ; : : 
ae : ; é : ; hibited by the wing considered. In the following ex 
ecessary | * direction—the interference drag may still be com- : é ; Bs 
: amples, the plate will be so selected that line EF co- 


edge ol | 


ing the 


puted for the case where \ = A(x) in the region DEF 
of the wing. 

The reflecting plate thickness need not have been 
considered as zero. A distribution of sources repre- 
senting its thickness may be added but, if its drag is to 
be computed with the integral relations, the restric- 


tions on X must be satisfied. Even if the thickness 


incides with the trailing edge of the wing and region 
ABC contains sources but not sinks. In Fig. 9, the 
shapes of the three interaction regions discussed above 
are shown for the case where the region ABC is tri- 
angular. It is seen that the chord of each region is the 


same. The dimensions s and d, needed for calculations, 


nan distribution does not fulfill these restrictions, the inter- are found from the geometry to be 

late on lerence drag of the wing can still be computed with s = (h/m) (1 — V1 — m2) (29 

trailing the method discussed. The drag due to this plate | 

antl thickness could be offset by appropriately recessing the d = 2Bh(1 — V1 — n,?) (30 
Such afterwing, if convenient. Alternatively, a stream 

| surface slightly above the reflecting plate in the con- DRAG OF SOME THREE-DIMENSIONAL BIPLANES 

| figuration under discussion may be selected as the 

Seow lower face of a reflecting surface lucing the same Pentagonal Wing 

. of the | - ¢ gs é producing the same 

ata el effect on the wing as the flat plate. A suitable dis- Consider the sweptback, tapered wing of pentagonal 


tribution of singularities must first be assigned to the plan form which is shown in Fig. 10. Spanwise sec- 


exterior plate face. This consideration has been dis- tions are dissimilar and the relative thickness dimin- 
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USE OF 


ishes along the span, vanishing at the tip. The geom- 
etry of the wing is defined by the parameters X, 1, c, 
| is the relative thickness at the root 


and 7, where 7 
standard 


section. The drag may be computed by 
methods of linearized theory and is given in reference 
14. The drag parameter CpB/7? is plotted vs. RX for 
n = 0, 1 in Fig. 11 for the case of symmetric distribution 
of thickness and zero angle of attack. 

A reflecting plate for this wing, designed in accord 
with the discussion of the preceding section, would be 
triangular in plan form with rounded apex; it could, of 
course, be extended farther into the undisturbed flow 
The plate and the 
other two interaction regions, ABC and DEF, 
In this example, the plate is so 


without affecting the interaction. 
are 
indicated in Fig. 10. 
oriented to the wing that EF coincides with the wing 
trailing edge. In addition, the plate is so sized that 
DEF is just circumscribed by the triangle of the after- 
wing plan form. These two assumptions give the 
maximum interaction amenable to calculation by the 
method discussed in the last section. The shape and 
orientation of each of the three interaction elements 
The configuration is suggested 
for the 
In this case all of 


are thereby specified. 
by two-dimensional considerations applied 
limiting case of n tending to zero. 
the disturbances produced on the plate side of the 
forewing are intercepted by the plate and reflected 
back to fully close the wing and balance the fore drag 


on the plate side. From the diagram 


h = [1 — Roc + 2sk]/2B (31) 


From Eqs. (29) and (31), 


s= (1 — R) (c/2B) 1 — V1 — nn?) nV 1 — n? (32 
From Fig. 10, 
| f dS = (Vid/rk) (Re — 2sk)? (33 
J ABC 
where \ = Oz Ox is the constant inclination in region 


ABC. From Eqs. (28), (82), and (33) the interference 


drag for the symmetric wing is found to be 


Cp, B/r? = —[1/R(2 — R)] [1 — (1 — R) X 
(l1— V1 — n?)/RV1—n?/;? (34) 
where Cp, = D,/(p1 Vi17/2)S, 
So = (Re®/k) (2 — R) 
r mg / oe 


Eq. (34) holds when the quantity in the bracket is non- 
negative. The lateral spreading of disturbances as 
they pass successively from region ABC to the plate 
to region DEF (which is characterized by the param- 
eter s) can become so large that the region ABC shrinks 
to a point for a fixed value of the span of region DEF. 
This condition occurs when the values of m and K are 


such that the quantity in the bracket of Eq. (34 
vanishes, and corresponds to 
Rc — 2sk = 0 (35) 


in Fig. 10. 


INTERFERING FL 
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The interference thrust, given by the negative of Eq. 
(34 
11. The total drag of the configuration for specified 


value of is given by the difference of the two curves. 


is plotted vs. Kk, for various values of m, in Fig. 


If a plate is used on both sides of the symmetric wing, 
the total thrust is twice the value shown in Fig. 11 
In the two-dimensional case, for which n 0, the can 
cellation of the wing drag on the side of the plate is com- 
plete, as is to be expected for this configuration. 
When » is in the neighborhood of 0.2 or 0.3, the reduc 
tion in drag is still of the same order; the thrust ob- 
tained with this configuration is very close to the two 
dimensional value. 

Consider now the same configuration excepting that 
all the thickness is regarded to be distributed over 
only the plate side of the wing, the value of 7 being un 
altered. The opposite side is assumed to be flat and 
uninclined. In the computation of the isolated wing 
draz, the interaction between the two faces of the wing 
at the tips must be considered. Using the methods 
of reference 15, the contribution of the tips to the drag 
of the isolated wing may be computed and ts given in ref 
erence 14. The total drag coefficient of the asymmetric 
isolated wing is found by adding tothis tipdrag coefficient 
twice the value of the drag coefficient of the symmet-ic 
isolated wing. It is found that the contribution of the 
tip drag to the total is virtually negligible (n < R 
The interference drag coefficient is four times that value 
given by Eq. (34). In Fig. 12, the total drag of the 
asymmetric wing with plate is compared to the drag 
of the isolated symmetric wing having the same total 
thickness. This basis of comparison is a conservative 
one because, for given values of n, Kk, and 7, the drag 
of the isolated wing is a minimum when the thickness 
is symmetrically distributed. In the two-dimensional 
instance, corresponding to n 0, the cancellation of 
the drag of the asymmetric wing by the plate is com- 
plete. When » is in the neighborhood of 0.2 or 0.3 
and F# of order 0.5, the drag is still much less than that 
of the isolated symmetric wing. For example, if 
0.3 and R = 0.5, the drag of the configuration is about 
25 per cent of the drag of the isolated symmetric wing 
or 12! » per cent of the drag of the isolated asymmetric 
wing. The range 0 < n < 0.3, in which interference 
of the biplane type is seen to be most effective, is par- 
ticularly appropriate to the higher Mach Number 
regime. 

The ratio of plate area to wing area S, S 
vs. R for several values of nm in Fig. 13. It is seen that 
the plate area is only about 30 per cent of the wing area 


is plotted 


with good interference. 


Triangular Wing 


The second problem considered is that of the inter- 
action between a wing with triangular plan form and a 
suitable reflecting plate (see Fig. 14). The wing has 
diamond profiles which are geometrically similar and is 
when 


characterized by the parameters 7, m, c, and 7 


the profiles are regarded as symmetric; the case of 
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Fic. 15. Symmetric triangular wing. 


supersonic leading edges—i.e., m, < .1—is of interest 
here. Again the drag is readily calculated by linearized 
methods* ' and is shown plotted vs. (1 — r) for several 
values of m; in Fig. 15 for the case of a symmetric thick- 
ness distribution. 

The design and orientation of the reflecting plate is 
dictated by considerations similar to the preceding 
example. In the case of small 7, the size and shape of 
the afterwing region having constant inclination gov- 
erns the size of the interaction elements. For maxi- 
mum interference under this restriction, the region 
DEF will be inscribed in the triangle whose altitude is 
rc and whose base is the wing trailing edge. The size 
and orientation of the three interaction regions are 


thereby determined, their chord being (rc — 4d). 
From Fig. 14(a), and Eq. (30), 
h = [(1 — r)c + d]/2B (36) 
d=cl—rni(- V1 — r°n,") Vi — rn? (37) 
From Fig. 14, 
| tao => (Vir, wk,) (rc — d)? (3S) 
ABC 


From Eqs. (28), (37), and (38), the interference drag 
for the symmetric wing is found to be 
Cp,B/r? = —(r'1—r) {ll —-(—-r7r) xX 
/ ee had <a 
(1 — V1 — r?m?)/V 1 — r?m?]? (39) 
when the quantity in the bracket is non-negative. 
Here the drag coefficient is based on the wing area 
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and the relations 


7r2(1 — 7) Ao = —T 2r +] 


Ay _= 


have been used. Eq. (39) holds as long as the after- 


wing governs; from Fig. 14 this condition may _ be 
written 
(7c —d) < (1 — rie {2 
or 
[l1—(l1—ry7r) xX 
ie i r°ny*)/V 1 — rn] < (1 — rr) (43 


If the upstream interaction region is to contain only 
compressive disturbances, then, for values of m, and 
r not satisfying Eq. (43), the forewing governs the size 
of the three interaction elements. The region ABC is 
then selected so as to be inscribed in the triangle formed 
by the leading edge of the wing and the line x = (1 
The line EF is still selected to be coincident with 


rjc. 
the trailing edge of the wing. There is then obtained 
h = rc/2B (44 
| fdS = (Vii/rki) A — 7)?e? (45 
7 ABC 
and Cp,b/7* = —(1 — 1r)/r (46) | 
Eq. (46) applies to the symmetric wing when Eq. (43 
holds with the inequality sign reversed. 
The plate thrust, given by the negative of Eqs. (39 
i 


and (46), is plotted vs. (1 — r) for several values of 7, 
in Fig. 15 and may be compared to the drag of the 
The total drag of the biplane 


isolated symmetric wing. 
By using a reflecting plate 


configuration is also given. 
on both sides of the wing, the reduction in drag could , 
be doubled. Fig. 15 shows that the relative wave drag 
reductions are not nearly as large as in the case of the ' 
pentagon wing. That the triangular wing is not par- 
ticularly suited for use with a reflecting plate may be 
seen from an application of two-dimensional methods 
to the limiting case of m, tending to zero and (1 — 7) = 
0.5. Only one-half of the compressive disturbances 
from the forewing on the plate side are intercepted and 
reflected back to the afterwing by the plate. Thus the 
plate eliminates only one-half of the two-dimensional 
drag of the plate side of the wing. This observation 
is borne out by Fig. 15. A similar discussion applies 
when ; is non-zero. 

The ratio of reflecting plate area to wing area is shown | 
lf(l —vr) = 


in Fig. 16 and is seen to be quite small. 
0.5 and nm, = 0.3, it may be seen from Figs. 15 and 16 
that the wave drag reduction obtained with one plate 
is about 23 per cent of the drag of the isolated sym- 
metric wing, whereas the area of the reflecting plate 
is about 26 per cent of the wing area. 

Due to the asymmetry of the configuration when 
only one plate is used, a symmetric distribution of wing 
thickness does not give the minimum configuration drag 
The minimum drag} 


for given values of m, 7, and r. 
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USE OF 


corresponds to a profile having more than half of the 
total thickness on the plate side. Let 7; be the con- 
stant relative thickness on the plate side of the wing and 
(r — 71) be the relative thickness on the other side, 
where 7 represents the total relative thickness [see 
Fig. 14(b) ]. 
39) and (46) with 7 replaced by 27;. 


The interference drag is given by Eqs 
The total drag 
of the configuration, written in terms of 7; and (7 — 7;) 
with the aid of references 16 or 8, may be minimized 
with respect to variations of 7/7. The total con- 
figuration drag, for optimum 7;/7, is plotted vs. (1 — 
r) for several values of m,; in Fig. 17. These results are 
compared to the drag of the symmetric isolated wing 
in this Figure. The optimum ratio of 7; 7 is also shown. 
By comparing Figs. 15 and 17, it may be seen that the 
configuration drag is somewhat reduced by the redis- 
tribution of the thickness. If (1 — 7) = 0.5 and n, = 
0.3, the drag reduction in the asymmetric case is about 
30 per cent of the drag of the isolated symmetric wing. 
It may be pointed ovt that the size and orientation 
of the plate remains unchanged by the redistribution of 


thickness. 
Triangular Wing with Pentagonal Protuberance 


One can design the reflecting plate large enough to 
intercept a major part of the disturbances originating 
from the forewing. A smaller reflecting plate, de- 
signed to intercept only disturbances produced near 
a possibly thick root section, can also be envisaged. 
Consider a symmetric triangular wing to which there is 
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added a pentagonal protuberance on one side in the 
vicinity of the root section. There results a wing with 
a variable spanwise distribution of relative thickness. 
The protuberance might be regarded as a storage vol 
ume. Such a configuration, obtained by geometric 
superposition of the wings discussed in the first two 
examples, is shown in Fig. 18. It is assumed that 
nm, < 1. If 7; represents the relative thickness of the 
symmetric basic wing and if 7,;; represents twice the 
maximum relative thickness of the protuberance on 
one side of the wing (see Fig. 18), then the wing is geo 
metrically characterized by the parameters 7, m, ¢, 71, 
and 7, where 

T = ry,/r1 (47) 
The isolated drag of the protuberance side of the wing 
has been computed and is given in reference 14. The 
drag of the opposite side may be found in references 16 
or 8. The complete drag of the isolated wing is plotted 
vs. (1 — r) for several values of m, in Figs. 19, 20, and 
21, which correspond to values of 7 of 0.5, 1.0, and 2.0, 
respectively. Also shown in each Figure for com- 
parison is the drag of the isolated wing without the 
protuberance. 

The size, shape, and orientation of the reflecting plate 
is identical with the first example, providing that the 
chord c is the same in each case. Then the distances 
h and s are given by Eqs. (31) and (32), respectively, 
wherein & is replaced by r, k by k, and by m,. From 
Fig. 18, 

. 


| fdS = Vy(A; + As) (rc — 2ski)? wk, (48) 
ABC 
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With use of Eq. (28), there results for the interference 
drag 
Cp,B/ry? = -—U + T) [r/ — 7] X 
(1+ 71 —7r)/r] [1 -—- (1-7) X 
/ 2\ /-2/ 272 
(1 — V1 — m?)/rV 1 — n?/? (49) 


providing the squared quantity in brackets is non- 


negative. 
The interference thrust, given by the negative of Eq. 
(49), is plotted vs. (1 — r) for several values of m in 


Figs. 19, 20, and 21, these curves being drawn for 
values of 7 of 0.5, 1.0, and 2.0, respectively. Also 
shown in these Figures is the total drag of the wing 
with the plate. These curves may be compared to 
those previously cited in the Figures which describe 
the isolated wing and the isolated wing without pro- 


tuberance. For the case m = 0 it may be seen, by 
comparing the three Figures, that the drag of the bi- 
plane is independent of 7 for all values of (1 — 7). 


This result is to be expected from two-dimensional 
considerations. If m; = 0.3 and (1 — 7) = 0.5, the 
Figures show that the drag of the biplane increases only 
slightly as 7 is increased successively from 0.5 to 1.0 
and then to 2.0, whereas the drag of the isolated wing 
increases a great deal when 7° is increased in the same 
way. When ™ = 0.3, (1 — 7) = 0.5, and T = 0.5 or 
1.0, the plate thrust overbalances the drag increment 
attendant upon the addition of the protuberance to 
the wing; when 7 = 2.0, the two opposing effects vir- 
tually balance. Fig. 12, after minor adjustment of the 
ordinate and abscissa, can be regarded as the limiting 
case 7 > 1. Fig. 15 does not, however, correspond 
to the case 7’ = 0, because the reflecting plate in the 
second example is larger than in the third. 

It is apparent that the thickness of the wing in the 
vicinity of the root may be significantly increased 
with little or no penalty in drag. Moreover, the re- 
flecting plate required is quite small. The ratio of 
plate area to wing area S,/S,, is plotted vs. 7 in Fig. 22. 
When m = 0.3 andr = 0.5, S,/S, = 0.236; if r is de- 
creased to 0.4, S,/S, falls to 0.147. Figs. 19, 20, and 
21 show that the increase in biplane drag when r is 
decreased from 0.5 to 0.4, for given m, and 7’, is not 


large. 
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PRACTICALITY OF THE THREE-DIMENSIONAL BIPLANE 


The virtual absence of the thickness drag of the 
symmetric two-dimensional Busemann biplane (see 
Fig. 5) has been known since 1935,! but a number of 
undesirable features appear to limit the usefulness of 
this configuration. Among these are numbered the 
sensitivity of the favorable interference to variations 
of the Mach Number from the design value and the 
possibility of choking in the internal passage. The 
maximum practical contraction ratio between the two 
wings is limited, as in the case of inlets, by the require- 
ment that the configuration start at, say, the design 
Mach Number”; this limits the thickness of the wings. 
Furthermore, the elimination of the monoplane thick 
ness drag is offset by an increase in other forms of drag 
when the symmetric biplane is substituted for the 
monoplane, providing that the two systems are com- 
The following comparison is pre- 


pared for equal lift. 
The two-dimensional profile 


sented in reference 18. 
at an angle of attack supports wave drag due to thick- 
ness, wave drag due to lift, and skin-friction drag; each 
form of drag is normally of the same order. If the 
thickness is suitably redistributed and a second wing is 
added to form a symmetric biplane, then the thickness 
drag is eliminated, the drag due to lift is unaltered, but 
the skin friction is roughly doubled because the wetted 
area has been doubled. If the chord of each wing is 
halved and the angle of attack of each is doubled to 
preserve the same lift, then this biplane exhibits no 
thickness drag, twice the induced drag of the original 
monoplane, and roughly the same skin-friction drag 
as the monoplane. If, for example, the three forms 
of drag of the original monoplane are all equal, then 
neither biplane cited offers an aerodynamic advantage 
over the monoplane. 

This assessment of the symmetric lifting two-dimen- 
sional biplane is misleading, however, because it con- 
siders only triangular profiles. A more efficient lifting 
two-dimensional biplane results when trapezoidal pro- 
files, such as appear in Fig. 10 for R < 0.5, are adopted, 
This modification allows some lift to be produced effi- 
ciently in the internal passage of symmetric thickness 
arrangements. Furthermore, a reflecting plate or other 
smaller wing can be substituted for the second symmet- 
ric wing, the former not conforming to the central 
streamline of the symmetric thickness arrangement. 
The chord of these auxiliary surfaces need be only a frac- 
tion of the chord of the basic trapezoidal wing. The 
proposed configurations readily excel the monoplane 
on the basis of the comparison of total drags used in 
reference 18 and repeated above. 

In the case of the three-dimensional biplane con- 
figurations herein considered, the possibility of elimi- 
nating most of the thickness drag of the isolated wing 
has been shown to remain. This provides the oppor- 
tunity to increase the wing thickness and storage vol- 
ume, especially near the root section, without a large 
drag penalty. But the three-dimensional biplane also 
exhibits the same undesirable sensitivity of the favor- 


— 


an = Lo in 


IPLANE 


of the 
e (see 
iber of 
ness of 
ed the 
iations 
nd the 

The 
he two 
equire- 
design 
wings. 
thick 
ft drag 
or the 
> com- 
is pre- 
profile 
thick- 
r; each 
If the 
ving iS 
ckness 
d, but 
vetted 
ying is 
led to 
its no 
riginal 
. drag 
forms 
, then 
untage 


jimen- 
t con- 
lifting 
1 pro- 
ypted, 
d effi- 
‘kness 
other 
nimet- 
entral 
ment. 
i frac- 

The 
plane 
ed in 


con- 
elimi- 
wing 
ppor- 
Ev y]- 
large 
> also 


avor- 





ete keET 


USE OF INTERFERING FLOW FIELDS 





DRAG OF ISOLATED —-— 


WING 
DRAG OF ISOLATED WING ---- 
WITHOUT PROTUBERANCE 
INTERFERENCE THRUST——— 
DRAG OF WING WITH 
PLATE 





ha 








aan 

















3G. 2 7 ng Tr y it} » FOI r -rance - 
Fic. 19. Triangular wing with pentagonal protuberance, 7 = rue. Tl, Tiengiier wing ia tts Sor ee: 7 
0.5. ‘ 
T |] 
25 1 ng 
2 
/ 3 
4 
/ 5 
20 Y, 
H{/ , 
15 ‘dl 





Se~ 





te 
—~ 





x 


ee 





— 


Wy 




















77 



































10 
05 
=) 
(1-r) 
Fic. 20. Triangular wing with pentagonal protuberance, 7 = Fic. 22 


1.0. 


| 2 3 4 5 
r 


Triangular wing with pentagonal protuberance: ratio 
of plate area to wing area 











16 JOURNAL OF THE AERONAUTICAL SCIENCES 


able interference to variations of the stream Mach 
Number from the design value which is characteristic 
of the two-dimensional biplane. This effect is a sub- 
ject for further research. The geometric limitations 
which are imposed by the contingency of choking in 
the two-dimensional configurations are largely elimi- 
nated in the case of the three-dimensional swept bi- 
planes discussed because spanwise flow is unrestricted. 
Therefore, large thickness distributions in the vicinity 
of the wing root may be used collaterally with a small 
gap. The reflecting plate is relatively small, the ratio 
S 

discussed, so that the increase in skin friction or drag 
due to lift which attends the reduction of thickness 
drag when the biplane is substituted for the isolated 


S, falling in the range 15-30 per cent for the cases 


wing need not be severe in the case of suitable lifting 
configurations. This point is illustrated by consider- 
ing the configuration illustrated in Fig. 10 for the 
limiting case where » tends toward zero. For equal 
angles of attack, for equal chords, and for all values 
of 2, the lift and pressure drag of this configuration are 
the same as the lift and drag due to lift, respectively, 
of the isolated symmetric pentagon wing. Thus the 
thickness drag of the latter is, in effect, eliminated 
without changing the induced drag; the increase in 
wetted area and, hence, skin friction, is only 33 per 
cent when R = 0.5. When x is non-zero, it can still be 
expected that the lifting biplane will exhibit sub- 
stantially less total drag than the isolated symmetric 
wing generating the same lift, whether the comparison 
be made on the basis of equal chord or equal wetted 
area. It is believed that a consideration of plate thick- 
ness does not substantially change these conclusions. 
Lifting three-dimensional biplanes will be considered 
by the authors in a subsequent publication. 


CONCLUSIONS 


A general outline has been given of possibilities for 
favorable interference between supersonic aircraft 
components. The interference between a three-dimen- 
sional wing and a small reflecting surface in biplane 
arrangement has been discussed in terms of linearized 
flow theory. Relationships have been developed which 
permit the analysis of such biplane configurations 
without the need for determination of the flow field. 
These relations apply in the linear approximation to a 
large class of such configurations. The method has 
been applied to several biplane arrangements of prac- 
tical interest. It has been shown that the pressure 
drag of the wing can be greatly reduced or that the 
volume near the root section of the wing can be con- 
siderably increased without penalty in pressure drag. 
The biplane has been discussed from the engineering 
point of view, on the basis of which it appears that the 
configuration might be a practical one for some appli- 
cations, particularly in the higher Mach Number range. 
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APPENDIX—ANALYTIC DERIVATION OF THE INTEGRAL 
RELATIONS 


Consider an elemental area AB in the plane z = h 
whose streamwise length is c, whose elemental width 
is dn, and whose orientation with respect to the coor 
dinate system is as shown in Fig. S(a). Let this area 
support a source distribution of constant intensity 
proportional to f. The trace in the plane z = 0 of the 
downstream Mach cone with apex at point .1(0, 0, A) 


is given by 
x = BV y? +h (A-1) 


Similarly, the Mach cone with apex at the point 
B(c, 0, h) produces the trace 


(x —c) = BV y? + h? (A-2) 


These hyperbolas are indicated as OER and PFQ, 
respectively, in Fig. S(a). The perturbation potential 
at a point P(x, y, z) due to a source distribution in the 


plane z = /is given by 


(x, y, 2) = 
T(E, n) dé dn 


— (A-3) 
If. V (x — §)? — B[(y — n)? + (2 — A)? ] 


where the integration is extended over the part of the 
distribution + which is upstream of the trace of the 
fore cone from P(x, y, z) in the plane z = 4. The equa- 
¢ 


tion of this hyperbola, in terms of the coordinates & 


and 7, is 
(x — §&) = BV (y —n)?+ (s —h)? (A-+) 


The elemental perturbation potential 6¢ induced by 
the elemental source distribution considered follows 


from Eq. (A-3) as 
by(x, ¥, 3) = 
dé 


fan | (A-5) 
0 V(x — £2 - B?[y? + (s — h)?] 


where £, differs according as P(x, y, 2) is upstream or 
downstream of the downstream Mach cone from 
Bic, 0, h). When [y? + (2 — h)*| is non-zero, it is 
seen that the elemental distribution behaves like a line 
source whose constant elemental lineal intensity 1s 
proportional to fdy. Therefore, the elemental area of 
width dy may be rotated about its streamwise axis 
without changing the flow field. If P(x, y, 2) is up- 
stream of the Mach cone from B(c, 0, /), then &; is given 
by Eq. (A-4), and it follows from Eq. (A-5) that 


x 
BV y? + (¢ — hy)? 
lye+(s —h)*]>0 (A-6 


6¢(x, vy, 3) = —fdn cosh! 


In Eq. (A-6) the principal value of the inverse function 
is inferred. By differentiation of Eq. (A-6), the x. and 
z components of the elemental perturbation velocity 
in the plane z = 0 are found to be 
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—fdn Vx? - B*(y? + h?) (A-7) 


—fdnhx/(y? + h?)V x? — By? + h?), 
(y? +h?) >0 (A-8) 


respectively. 

We now undertake the derivation from Eqs. (A-7) 
and (A-S) of several integral relations that serve to im- 
plement certain previously discussed calculations and 
to exhibit some features of the flow field considered. 
It is convenient in this regard to introduce the oblique 
coordinates a and 6 in the plane z = OQ, the axes of 
which are inclined by the angle tan~'k with respect 
to the y axis as shown in Fig. S(a). The coordinates 
a and £ are related to the coordinates x and y by the 


inverse transformation 


x = (k V1 4+ R) (a + B) (A-9) 


y = (1/V 1 + R?) (8 — a) (A-10) 
By substitution of Eqs. (A-9) and (A-10) into Eqs. 
(A-7) and (A-8) there results 
é6u(a, B, O) = 
V1 + R? fdn 
BV na + 8)? — (8 — a)? — hX(1 + &2) 


(A-11) 


éw(a, 8B, O) = 
(1 + k*)nh fdn(a + B 
~ [(8 — a)? + 1 + k?)| X 
V n?(a + B)? — (8 — a)? — AX(1 + R?) 
(A-12) 


where n = k/B. 
inclined with respect to the positive y axis by the angle 


Consider two parallel lines which are 


tan ~'k, and which are tangent at points E and F, respec- 
tively, to the hyperbolas given by Eqs. (A-1) and (A-2); 
these lines are identified by EE’ and by DC in Fig. S(a). 
From the geometry, the equations of these lines are 
found to be 
ag = (V1 + k? 2k) BhV 1 — n? (A-13 
ap = (V1 + k2/2k) (Cc + BhV 1 — m2) (A-14 


respectively. 
We first evaluate the integral of 6u and of é6w along a 


generic line a = const. between the two branches of 
the hyperbola OER; this line segment is denoted by 
KJ in Fig. S(a) and lies in the region CED. That is, 
ap < a< arp (A-15) 

From Eq. (A-11), 

has Vit k fd 

6u(a, B, 0) dB = — = x 
J Bx B 

*Bs dB 

7 (A-16) 
Jee Vn*(a + B)? — (8 — a)? — h2(1 + ?) 


where 8x and #, are found by substituting Eqs. (A-9) 
and (A-10) into Eq. (A-1) and solving the result for 8 


in terms of a: 
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(1 + n?) | 
Bey = ei. = 
(1 — n°) 1 — n° 
V 4tn7a? — h2(1 + Rk) 1 — we? A-17 


After rearrangement of the radicand, the integral in 
Eq. (A-16) is readily evaluated by trigonometric sub 


stitution. There results 
Jeni 
(1/V 1 + ef 6u(a, B, O) dB = —rfdn BV | n° 
e’ DOK 
A-18 
Similarly, from Eq. (A-12), 
AJ 
| 6w(a, B, O)dB = —(1 + k*)nhfdn X 
/ BK 
[- \a@ + B)dB 
J px (8B — a)? + h2(1 + R?)] X 
V n*(a + B*) — (8B — a)? — h2(1 + R? 
(A-19) 
After making the transformation 
y = (8 — a) (a + B) (A-20 


the integral in Eq. (A-19) may be evaluated through 
use of two subsequent trigonometric transformations. 
There results 

6w (a, B, O)dB = 


(l1/V 1 + k?) — mfdn, h+#0O 


PR 


A-2] 


Since d8/V 1 + k? = dy, Eqs. (A-18) and (A-21) corre- 
spond to Eqs. (16) and (17) in the main text. It is 
significant that the value of each integral is independent 
of a. 

We next evaluate the area integrals of éu and é6w 
over the region CED [see Fig. S(a)]. Since the ele 
mental area is given by 
2?) |dadB (A-22 


dA = [2k/(1 4 


it follows that 


bd Fay eAy 
6udA = [2k/(1 + k?)] da | 6u (a, B, O) dg 
CED J at 


. ¢~ BK 

A-23 
The first integral is given by Eq. (A-18); the second 
integration is trivial. From Eq. (A-23) 


6udA = af cdn’ BV 1 — n? 
J CED 
Similarly, 


” *arF bd I 
} ébwdA = [2k (1 + k?) | da dw(a, 8, 0) d3 
7 CED 


e/ at e/ BK 
A-25 
With use of Eq. (A-21) there results 
| bwdAi = — nf cdn, hk #0 (A-26) 
7 CED 


Eqs. (A-24) and (A-26) correspond to Eqs. (12) and 


(13) in the main text. 
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We lastly discuss the line integral of induction along 
a generic line ay < a = const. between the branches of 
the hyperbola OER. 
NM in Fig. (8a). 
elemental distribution AB by a semi-infinite elemental 
source distribution beginning at A and a semi-infinite 
Then the 


This line segment is denoted by 
It is convenient to represent the 


elemental sink distribution beginning at B. 
induction due to the former in the entire region down- 
stream of hyperbola OER is given by Eqs. (A-7), 
(A-8), (A-11), and (A-12). The latter 


induces a flow in the entire region downstream of 


distribution 


the hyperbola PFQ. This flow is also described by 
Eqs. (A-7), (A-8), (A-11), and (A-12), with f replaced 
by —/f, in terms of a second coordinate system whose 
origin is shifted downstream relative to the first one 
by a distance c. The complete flow in the region down- 
stream of hyperbola PFQ is then given by superposi- 
tion. Because the values of the integrals in Eqs. 
(A-18) and (A-21) are independent of a, it follows that 
the integral of induction from N’ to M’ due to the 
semi-infinite sink distribution is equal and opposite to 
the integral of induction from N to M due to the semi- 
infinite source distribution. Then the induction due 
to the elemental distribution AB for all a > ap is such 


that 


eBM 
| 6u(a, B, 0) dg = 0, (A-27) 
e/ BN 
%BM 
| Ow (a, 6, Q) dp = 0, h A~ UV A-28) 
ed BN 
Eqs. (A-27) and (A-28) correspond to Eqs. (14) and 


(15) in the main text. 
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An Experimental Study of the Response of 


Simple 


Panels to Intense Acoustic Loading 


LESLIE W. LASSITER,* ROBERT W. HESS,* ann HARVEY H. HUBBARD” 
Langley Aeronautical Laboratory, NACA 


SUMMARY 


his paper deals with the measured response of simple panels 


to intense noise of the type encountered in the vicinity of aircraft 


power plants. Stress data are presented for panels exposed to 


noise from a jet engine at levels of approximately 160 db. as well 
is to nominal levels of discrete and random noise in controlled 


laboratory tests. Even at noise levels as high as those encoun 


tered near jet engines the data are shown to be in general agree 


ment with calculations by the power spectrum procedure sug 
gested by Miles 
Damping is shown to be a particularly significant parameter 


ind is evaluated for various test conditions. Opportunity is 


taken to discuss the techniques of testing with a high intensity 


laboratory siren and to present some noise spectra measured 


near the exhaust of a turbojet engine Results of an exploratory 


investigation of fatigue life and crack growth are also presented 
for various panels which failed under discrete frequency acoustic 


loading. Fatigue life was noted to decrease markedly for a 


nominal increase in the noise intensity level 


SYMBOLS 


frequency 


w damped natural frequency of panel 
g\lw power spectrum of random noise 
Z(w) = mechanical impedance of panel 
Yah power spectrum of panel stress respons¢ 
a” mean square stress 
Vo root mean square stress 
omar. = stress amplitude 
Ost stress under given static load 
\P static pressure differential 
6 damping ratio, ¢/¢ : 
S stress per unit static pressure 
t panel thickness 
R panel radius of curvature 
root mean square sound pressure at siren funda 
mental frequency 
D diameter of engine tailpipe 
slant distance from tailpipe (along lines 15° to jet 
centerline 
d radial distance from 15° jet boundary 


characteristic panel dimension 


INTRODUCTION 


penn DAMAGE to the skin structure of airplanes 
and missiles due to noise from the power plants is 
now a serious operating problem. Of particular 
concern are those failures due to the intense random 
noise fields near jet engines. There is a need for evaluat- 
ing the significant parameters involved in the problem 
in order to permit the design of optimum structures to 


withstand these fluctuating loads. 


Presented at the Structures Session, Twenty-Fourth Annual 
Meeting, IAS, New York, January 23-26, 1956 
Scientist 


* Aeronautical Research 


The present paper is devoted mainly to a discussion 
of the response of simple panels to intense random noise 
inputs. The concepts of generalized harmonic anal 
ysis are applied to this latter problem, and the results 
are compared with experiments for both flat and curved 
panels. As a matter of interest some sample spectra 
of the noise near a J-34 turbojet engine with and with 
out afterburner are presented. 

There are brief discussions of some exploratory tests 
of the fatigue of simple riveted panels by means of 
discrete frequency noise from a laboratory siren 
The manner in which these failures occur is illustrated 
and some indication of how the fatigue life varies with 


intensity of the input noise is given. 
(I) PANEL STRESSES 


METHOD OF CALCULATION 


The concepts involved in analysis of the problem of 
panel excitation by random loading are outlined and 
illustrated in Fig. 1. A random acoustic input or loading 
which may have a power spectrum such as that shown 
in the top part of Fig. 1 exists, and it is desired to deter 
mine the power spectrum of the stress response which 
this loading will produce in the panel. The two spec 
tra are related by the transfer characteristic of the 
panel, which is the square of the admittance or fre 
quency response to a unit-pressure sinusoidal input and 
in simplified form may be as shown in the middle part 
of Fig. 1. The power spectrum of the stress response 
then results from the product at all frequencies of the 
input spectrum ¢y(w) and the square of the admittance 

1/Z(w such that 


the half-power bandwidth of the panel admittance curve 


If the system is sharply tuned 


is small relative to that of the loading spectrum, the 
the 
spectrum takes the same form as the admittance curve. 
the 


a simple manner if the panel is as 


response ¢,(w) is essentially resonant and stress 


As Miles and others" * have shown response can 
be calculated in 


sumed to behave as a single degree of freedom system 


which is sharply tuned. Thus, the square of the 
impedance |Z(w) * is taken as 
Z(w) - ; [1 — (@/ Wo) |? ms 16 Ww Ww : h \ l ) 


is taken as the power spectral 
It should be 
noted that in this simplified presentation of the prob- 


and the input ¢y(wo 
density at the panel resonant frequency. 


lem, the spatial correlation of pressure over the panel 
area is tacitly assumed to be unity, that is, all parts of 
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the panel are assumed to be excited with pressures in 
the same phase. Where the inputs are of a scale which 
is not large relative to the panel dimension this condi- 
tion is not met, and the input loading spectrum must 
be one in which correlation effects have been included. 
Such effects, however, will not be considered in the 
present paper. 

The mean square response go”, which is a measure of 
the area under the stress power curve, becomes 


o? = 1/46 ao ¢n(wo) So? (2) 

The mean square stress -s thus seen to be dependent 
upon four parameters, three of which relate to the 
panel. wo is the damped natural frequency of the 
panel, So is the stress per unit input pressure under 
static loading, and 6 is the damping in terms of the 
critical value. The admittance curve of the panel 
affords a measure of damping, either by its bandwidth 
26a at the half-power frequencies, or by the value 
(So 26) at resonance. The input ¢y(wo) is the power 
spectral density at frequency wo. 

The mean square stress is a useful index of the panel 
but still another characteristic of the re- 
sponse is of interest. This is the time history of stress. 
Fig. 2 gives a comparison of stress time histories of a 


response, 


given panel for periodic loading and for random loading. 
For periodic loading, with one Fourier component of 
the noise near the panel natural frequency, the stress 
response is essentially sinusoidal at that frequency and 
the envelope is of constant value. 

For random loading Fig. 2 indicates that the response 
is again essentially sinusoidal, but in this case the 
envelope is not constant. In fact, it can be shown that 
the amplitude of a single cycle of the response is a 
statistical quantity having the so-called Rayleigh prob- 
ability distribution. This type of response is char- 
acteristic of a sharply tuned system randomly excited, 
and the interval between the beats indicated in Fig. 2 
is an indication of the damping of the system. 


EXPERIMENTAL RESULTS 


The remainder of the present paper will deal with (1) 
determination of the structural and loading parameters 
necessary to a calculation of the mean square stress, 
(2) a comparison of measured stresses with those calcu- 
lated by Eq. (2), and (3) results of some exploratory 
fatigue tests of panels under acoustic loading. The 
panels used in the investigation were of 2024 t3 alu- 
minum alloy, measuring 11 by 13 in. with simulated riv- 
eted-edge conditions. The panels were mounted on a 
rigid plate of 1l-in. thickness and this was subsequently 
bolted to a supporting chamber. Flat panels of thick- 
nesses (0.032 in. and 0.064 in. and curved panels of 
thickness 0.032 in. are considered. The curved panels 
consisted of flat stock rolled to a radius of 4 ft. 

Panel Ckaracteristics 

Static Response—Fig. 3 illustrates the static response 
of flat and curved panels of 0.032-in. thickness. Static 
stress is plotted as a function of pressure differential 
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across the panel for loading in both directions. For 
flat panels the response indicates that for loading in 
either direction the system is linear over a range of 
pressure and then nonlinear in the char 
acteristic manner of a system of increasing stiffness with 


becomes 


deflection. 
dependent upon the direction of loading. 
pressure on the concave surface the response is similar 
to that of a flat panel. For negative-pressure loading 
on the same surface, however, the response curve turns 
upward, indicating that the nonlinearity involves a 


For curved panels, the static response is 
For positive 


decreasing spring stiffness. 

Dynamic Response--The dynamic response of the 
various panels was determined by exciting the panel 
with periodic noise from a large siren. By holding the 
acoustic pressure output of the siren constant and 
varying the frequency, an admittance curve for that 
pressure level is obtained. Fig. 4+ gives sample results 
of the measured admittances of flat and curved panels 
of 0.032-in. thickness. Stress amplitude is plotted as a 
function of frequency for two input pressure levels to 
the two panels. Although the levels used with the 
flat and curved panels are not the same, the curves will 
illustrate the main differences between the two types. 
At the lowest pressure for each panel the curves are 
similar and tend toward symmetry about the natural 
frequency, which is observed to be appreciably higher 
for the curved panel than for the flat panel. Actually, 
the curved panel response shown here is for a mode very 
much like the third of a flat panel (node lines in 
direction) and the flat panel response is for the 
These were chosen because in each case 


long 
first 
mode. this 
particular mode was found to be the primary one ex- 
cited when the panel was forced by a random input. 
These are not necessarily typical since the modal re- 
sponse probably depends upon the radius of curvature 
of the panel and upon the spatial scale of the loading 
disturbance also, but they are of interest here because 
of the relation they have to the results obtained under 
random excitation. 

In contrast to this similarity of response of flat and 
curved panels at low input pressures, the top curves 
of the Figure illustrate the differences which occur at 
high levels of excitation. The flat panel assumes a 
response which is skewed toward higher frequencies, 
whereas the curved panel response is skewed toward 
lower frequencies. These two types of response are 
indications of two types of nonlinearity.* For the flat 
panel the nonlinearity arises from an increasing stiff 
ness with deflection, and for the curved panel at this 
particular input level the nonlinearity arises from a 
condition of decreasing stiffness with deflection. At 
intermediate levels the curved panel also indicates a 
slight stiffening effect. In practice, the double-valued 
portions of these skewed curves are nearly impossible 
hence, the half-power frequencies are some 
For that reason, the 


to obtain 
what a matter of conjecture. 
determination of damping 6 was based on the resonant 
amplification, using the relation 
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Panel Damping—Fig. 5 illustrates the nature of the 
damping results obtained from one of the test panels. 
The Figure gives the variation of damping with stress 
(or deflection) for a 0.064-in. flat panel for two condi- 
The 
top curve is for termination in the mounting chamber 
The lower 


tions of acoustic termination on the back side. 


of material of high acoustic absorptivity. 
curve is for termination with a material of appreciably 
lower absorptivity, and it is seen to indicate damping 
values of roughly half those of the other case. In neither 
case does the material actually come in contact with the 
panel—hence, the difference in damping is taken as 
evidence of the importance of radiation losses to this 
type of test mounting. For either case, the damping 
increases with stress at the higher stresses, but when 
the termination is highly absorbent the increase is at a 
much higher rate. 

Damping for the curved panel was found to be lower 
at low stresses and higher at high stresses than the 
values for the flat panel shown here, thus following a 
somewhat different variation with stress. 

The damping data of this Figure, together with the 
static stress response and natural frequency from Figs. 
3 and 4, respectively, provide the information necessary 
for a calculation of the mean square stress for a given 
input noise spectrum, and it thus remains to compare 
measured values with those calculated. Before pro- 
ceeding to this, however, it might be well to consider 
some of the features of a turbojet noise field. 


Some Characteristics of the Noise Inputs 


engine 
the 


The present tests were made with a J-54 
equipped with afterburner. This is not one of 
more powerful engines, but some of the characteristics 
of its noise field are expected to be more or less typical. 

Spectra-—Fig. 6 gives sample spectra obtained near 
the flow boundary for this engine at the 100 per cent 
r.p.m. condition and for afterburning conditions. 
Sound pressure is given as the spectrum level value in 
the units psi Yeps. At 100 per cent r.p.m. the spec- 
trum shape is more or less typical of those obtained 
heretofore in the near sound field of an engine in that 
the highest pressures occur in the frequency range 100 
cps to 1,000 eps. When the afterburner is cut in, the 
spectrum level is raised by about 8 to 10 db. as a result 
of the higher efflux velocity of the jet—-root mean square 
sound pressure varies about as the second power of jet 
velocity in the near sound field.! In addition to this 
increase, Operation of the afterburner introduces a very 
intense periodic component into the spectrum at 125 
eps. This apparently results from an internal reso 
nance of the tailpipe and is not necessarily typical of all 
afterburners. In the afterburner tests discussed herein 
the test panel had 
and the stress time histories were similar to that of 


Fig. 2 for random input, thus indicating that the load- 


a resonant frequency of 160 eps 


ing was primarily random. 
Spatial Variations—The quantity of primary inter- 


est from these spectra is the value at the panel natural 
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o . ‘ l 2 rie . . ° 
frequency, that is [¢gy(wo/2mr)]~. This characteristic 


of the problem allows a fairly simple means of charting 


the spacewise variations of sound pressure near the 


engine, since several pertinent frequencies can be 
selected and the spatial distribution of pressure at those 
gives such plots for fre- 
The left- 


hand plot is a plot of root mean square sound pressure 


frequencies plotted. Fig. 7 


quencies of 200 eps, SOO eps, and 1,600 eps. 
lyy(ao 2a)]°* as a function of slant distance z/D meas- 
ured from the tailpipe along a line inclined at the angle 
of the jet boundary (15°) and displaced radially one 
diameter from the boundary. The lower frequencies 
are seen to reach their maximum value at some distance 
downstream of the tailpipe (11 diameters for 200 cps). 
At 800 eps the maximum occurs at about 7 diameters 
and at 1,600 eps the maximum appears to be upstream 
of 2.5 diameters. 

The right-hand plot of Fig. 7 illustrates the relative 
variation of sound pressure for 200 cps and 1,600 eps, 
with radial distance d/D from the jet boundary at a 
slant distance of 8S diameters downstream of the tail 
pipe. The decrease of pressure with distance from the 
boundary is rather rapid for 1,600 eps but for 200 eps 
the decrease is at a comparatively slow rate. The 
trends illustrated by this Figure are thought to be 
more or less typical of turbojets in general, although 
the magnitudes involved are undoubtedly lower than 


exist near some of the latest model turbojets. 


COMPARISON OF THEORY AND EXPERIMENT 


The measured stresses to be compared with theory 
herein were obtained by selecting locations in the 
turbojet noise field such that a range of input pressures 
the lowest 


could be obtained. The measurements at 


pressures were obtained in the sound field of a 4-in. 
unheated-air jet... The damping used in a stress cal 
culation for random input can be obtained from Fig. 5 


by an iterative process; however, a more direct pro- 


cedure for obtaining the calculated stress without 
the necessity of determining the damping as such is 
available. This involves rewriting Eq. (2) as 

Vv o- Vv 6 = [ar t wo CN Ww So? ] . (3 


Then, based on the damping curve of Fig. 5, the quan 
tity V/ ,26 is plotted as a function of ~W>2 as shown in 
Fig. 8S. This allows entry into the ordinate scale with a 
value computed from the right side of Eq. (3), and the 
calculated stress is obtained directly from the Figure. 
Figs. 9 and 10 give comparisons of measured stresses 
with values calculated in this manner for flat and curved 
panels for a range of input pressure [¢y(wo/27) | 
Flat Panel 
in. and 0.064 in. 
stress variations for each panel; the data points repre- 


Fig. 9 is for flat panels of thickness 0.032 
The curves represent the calculated 


sent measured values. 

The indications are that the theory is in rather good 
agreement with experiment throughout the range of 
the tests for both panels, despite the fact that a mark- 
edly nonlinear response is associated with the higher 
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pressures. It will be noted, however, that the non- 
linearities are at least partly accounted for in the em- 
pirical values used for the parameters. 

It should also be noted that the nonlinearity in Fig. 9 
is partly due to the fact that the backing condition for 
the low-pressure runs differed from those of the turbo- 
jet runs, thus yielding different dampings even for the 
same input pressure. Although the 0.064-in. data in- 
dicate that the theory is somewhat unconservative at 
the higher pressures, similar comparisons for other 
thicknesses do not always indicate this. Thus the dif 
ferences appear to be due to experimental scatter or 
perhaps to a systematic error. 

Although the tests with the 0.032-in. panel were lim- 
ited to the 100 per cent engine condition because of the 
acoustical resonance of the afterburner, an evaluation 
of the effect of panel thickness can be made at the lower 
pressures. If in Eq. (2) the quantity wo 5S)? is substi- 
tuted for in terms of panel dimension a and _thick- 
ness /, there results the proportionality relationship 


Ver aa [ew(wo)]/*/8” t°”* (4) 
which indicates that root mean square stress should 
vary inversely as thickness to the 3.2 power. Thus it 
would be expected that a doubling of thickness, as in 
Fig. 9, would result in a stress reduction of 0.37, as is 
shown by the two theory curves. The general agree- 
ment of experiment and theory in each case appears 
to confirm the calculated effect of thickness. 

Curved Panels Fig. 10 gives a comparison of stress 
for panels of given size and thickness (0.052 in.) having 
different radii of curvature; and it allows evaluation of 
the theory for the curved panels. The top curve is 
taken from the preceding Figure for a flat panel (R = 
o); the lower curve is the calculated stress for a 
curved panel of 4-ft. radius. The measured stresses 
in the latter case are for excitation by the 4-in. air jet 
only, hence the more limited pressure range. Over 
this range the theory and experiment are again in 
rather good agreement. It may be that the agree- 
ment for this panel is somewhat fortuitous, since, as 
previously discussed, the panel responds in a mode 
much like the third, and the value of Sy is that for uni- 
form static loading in the nature of a first mode re- 


sponse. 


(II) PANEL FATIGUE CHARACTERISTICS 


From the results of Figs. 9 and 10, it appears that a 
power spectrum approach to the stress problem of 
excitation affords a_ usable 
However, it still remains 


panels under random 
means of stress prediction. 
to investigate the fatigue properties of the panels under 
acoustic excitation. 

In view of the differences indicated by Fig. 2 for 
panel response to periodic and random inputs, it would 
be expected that for a given panel and mean square 
stress there would be a difference in fatigue life in the 
two cases. Thus the specific fatigue properties of a 
given type of panel should be determined with random 
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loading. This in itself presents serious problems of 
testing, since the generation of random noise is by 
nature an inefficient process and to date the best gen- 
erators are the turbojets themselves. These are not 
generally available for an extended program of fatigue 
testing. However, the generation of intense periodic 
sound with a siren is economical, and it appears that 
considerable information can be obtained from fatigue 
testing with a siren. This is particularly true in the 
investigation of trends, such as was undertaken in the 


present tests. 


TYPE OF FAILURES 


For 0.032-in. flat panels with simulated riveted- 
edge conditions, Fig. 11 illustrates some characteristics 
of the failures obtained under periodic acoustic ex- 
citation. The regions identified as (1) in the photo- 
graph are typical of the first stage of the failures, 
wherein small cracks originate around rivets where 
stresses are highest. The time of appearance of such 
cracks is taken as the fatigue life of the specimen. 
Further excitation of the panel usually, in a rather short 
time, results in growth of the cracks around individual 
rivets until the cracks from adjacent rivets meet and 
form a split running nearly the length of a side as shown 
by the region denoted as (2). Once this occurs, the 
subsequent damage usually takes the form of new 
cracks extending toward the center of the panel and 
often growing in such a way as to leave holes in the 


panel [see region (:3) ]. 


EFFECT OF SOUND LEVEL ON FATIGUE LIFE 


Fig. 12 illustrates the importance of the sound level 
to the problem of fatigue life of a panel. This Figure 
presents a plot of the sound level in decibels as a func 
tion of fatigue life in minutes for flat panels of 0.032-in. 
At 140 db. a single sample was tested for 
At 154 db. a family of 


thickness. 
over 500 min. without failure. 
seven panels yielded a mean fatigue life of 14 min. and 
at 160 db. a family of 5 panels had a mean life of 3 min. 
At 166 db. two specimens indicated that the life was 
only of the order of 30 or 40 sec. Thus, it is observed 
that an increase of the decibel level from 140 db. to 
166 db. decreased the fatigue life of these particular 
panels from a matter of hours to a matter of only a 
few seconds, or by a factor of roughly 500. Although 
the data shown in Fig. 12 apply directly only to the 
particular configuration tested, the trend illustrated by 
the curve appears to be more or less general. 


CONCLUDING REMARKS 


Consideration has been given to the problem of the 
response of simple panels to intense acoustic excitation. 
Application is made of the method of generalized har- 
monic analysis to the problem of prediction of stresses 
due to random loading and the results are shown to be 
in good agreement with measured data. Exploratory 
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Stagnation Point of a Blunt Body in 
Hypersonic Flow 


TING-YI LI* ann RICHARD E. GEIGER** 
Rensselaer Polytechnic Institute and General Electric Company 


SUMMARY 


Phe of this paper is to present a method of calculation 


devised to yield all the important information on the symmetric 


purpose 


inviscid hypersonic flow in the stagnation point region of a blunt 
body The problem is the same as that considered by Hayes? 
who used a slightly different approach. It is demonstrated that 
Hayes’ results are valid in the stagnation point region and can 
hence be considered a basis for constructing less restricted solu 
tions 

Equations are presented giving velocity, pressure, 


The 


pressure c¢ efficient 


detachment 
distance, and vorticity values of shock detachment dis 


tance and body are compared with experi- 


mental data for spheres. The pressure comparison shows that 
the results of Hayes and the theory presented herein represent a 
better approximation than the Newtonian impact theory for 
hypersonic Mach Numbers 

In conclusion, the possibility of refinements to this analysis is 


discussed. 


SYMBOLS 


1,, 4), 49 = constants of integration, two-dimensional prob 
lem 
Co, Ci, ¢ = constants of integration, three-dimensional prob 
lem 
( = pressure coefficient 
= velocity function, defined by u = xf’(y 
A = free-stream density divided by density behind 
bow shock at body axis 
K = body curvature, positive for convex body 
= pressure 
= distance between body axis and arbitrary point 
vy 
R = radius of curvature, negative in the present 
analysis 
= velocity component in x-direction 
= velocity component in y-direction 
\ = curvilinear distance along body surface 
= distance normal to body surface 
a = flow deflection angle immediately behind shock, 
relative to free stream 
= shock angle relative to free stream 
= shock layer thickness 
5* = shock layer thickness on the axis of symmetry 
kinematic viscosity 
vorticity number, (w/ug) V vx /upg 
= density 
= angular position on sphere surface, relative to 
free stream 
A = body surface slope relative to free stream 


OW / OX (Ou/OY) 


= vorticity, 


Presented at the Hypersonics Session, Twenty-Fourth Annual 
Meeting, IAS, New York, January 23-26, 1956 

* Consultant to Special Defense Projects Department, Gen- 
eral Electric Company, and Associate Professor, Department of 
\eronautical Engineering, Rensselaer Polytechnic Institute. 
Defense Projects Department, General 


** Engineer, Special 


Electric Company. 


Subscripts 
B = body 


S = shock 
= free stream 
2 immediately behind normal shock 


INTRODUCTION 


bn PROBLEM of the high-speed flow past a blunt 
nosed body has been a subject of considerable 
research. In this problem the bow shock wave is 
detached from the body surface. Heybey' recently 
proposed a method of analysis designed particularly 
for an estimation of the shock detachment distance in 
the stagnation point region. In his paper Heybey 
cited many interesting references. In particular, he 
mentioned that Schwartz suggested that for hyper- 
sonic flows the detachment distance should be uniquely 
dependent upon the density ratio across the shock 
wave. The important role of the density ratio param 
eter is even more clearly brought out by the recent 
investigation of Hayes" Hayes has calculated’ 
the symmetric hypersonic flow past a 


body, and his results, presumably, can be considered 


blunt-nosed 


as a first step towards the understanding of the charac- 
This paper approaches the 
The 


present analysis shows that Hayes’ results are valid 


teristics of such flows. 
Hayes’ problem in a slightly different manner 


in the stagnation point region and, hence, can be 
considered a basis for constructing solutions for other 
less restricted regions. 

An understanding of the inviscid hypersonic flows 
past blunt bodies is important from a practical point of 
view. Presumably the inviscid analysis can provide 
an accurate description of the main stream conditions 
outside the viscous boundary layer. Analysis of the 
boundary-layer behavior, then, will provide important 
design information on skin friction and heat transfer 
The heat-transfer rate near the stagnation point of a 
blunt body, in particular, is important for designers. 
It is known that this heat-transfer rate depends on 
the magnitude of the velocity gradient along the body 
surface.*:> It has also been shown recently that the 
hypersonic boundary-layer behavior in the presence of 
a curved bow shock depends upon a “vorticity num- 
ber." 

These considerations clearly point out the impor- 
tance of an inviscid theory that provides a more or 
less complete description of the disturbance flow field. 
This is the main objective of the present study which 
presents a method of calculation devised to yield all 
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the important information on the inviscid hypersonic 
flow in the stagnation point region. 

Hayes’ method of small density ratio represents a 
better approximation than the Newtonian impact 
theory for finite hypersonic Mach Numbers. This is 
further substantiated in the present computations as 
is shown by a comparison of the pressure coefficient 
on the body surface according to these results with the 
classical Newtonian value. 

Results of the present study are compared with avail- 
able experimental data, and generally speaking, the 
comparisons show good agreement. Further refine- 
ment of the present analysis is discussed. 


ANALYSIS 


General Theory 
The purpose of the present investigation is to devise 
a method of computing the steady symmetric inviscid 
flow field about a body with a blunt nose at hyper- 
sonic free-stream Mach Numbers (\/.. > 1). 
The basic equations of motion governing the flow 
of an inviscid fluid are as follows: 
[1/(1 + Ky)] u(Ou/ox) + 
v(Ou/Ov) + [A/(1 + Ay)] ww + 
[1/(1 + Ky)] (1/p) (Op/ox) = 0 (1) 
[1/1 + AKy)] u(dv/Ox) + 
v(dv/Oy) — [A/(1 + Ay)] wv? + 
(1/p) (Op/oy) = 0 (2) 
The continuity equation is 


[1/11 + Ky)] (0/dx) (pu) + (0/Oy) (pv) + 
[A/(1 + Ky)]pv + (7/r) } [ou/(1 + Ky)] (Or/Ox) + 
pv(Or/dy)f = 0 (3) 
\0 two-dimensional flow 
1 = 
F axially symmetrical flow 


In these equations x and y are the axes of an orthogonal 
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curvilinear coordinate system, in which the origin is 
taken at the nose of the body, x is measured parallel 
to the body surface, and y is measured normal to the 
body surface. wu and v represent the velocity com 
ponents in the x and y directions, respectively Sec 
Fig. 1.) 

A (in general, a function of x) represents the local 
curvature of the body surface: 


K = —(d0/dx) = — (1/R, 


6 is the local inclination of the body surface with respect 
to free stream, and R, is the local radius of curvature 
p and p denote, respectively, the 


ry is the distance from an arbi- 


of the body surface. 
density and pressure. 
trary point x, y measured normal to the axis of sym 
metry which is parallel to the free stream. 

Eqs. (1) to (3) in their present forms are not amen- 
able to simple mathematical analysis. However, in 
the nose region they are susceptible to great simpli 
fications, as follows: 

(1) Near the stagnation point u u_, v/u, = O [R|I. 
For hypersonic flow, k—the density ratio across the 
shock wave—is very small, so that higher order terms 
like Auv, Au? can be neglected. 

(2) In the disturbed flow region between the shock 
wave and the body surface—i.e., 0 < y < 6, where 6 
represents the shock layer thickness—the quantity 
Ky = O[6/R,] can be neglected if 6 R, < 1. For 
the hypersonic flow problem, 6 Ry < 1 is, of course, 
a reasonable estimate. 

(3) In the vicinity of the stagnation point, 7 > x. 
Therefore, for small x, ry can be replaced by x in Eq. (3). 

(4) As the free-stream Mach Number approaches 
infinity, the Mach Number immediately behind a nor- 
quantity 


mal shock wave, j/2, approaches the 


V(y — 1)/2y- Therefore, the flow field near the 
stagnation point can be regarded essentially as in- 


compressible. 
p™ pp. = p kR (5 


For strong shock waves produced by blunt bodies at 
hypersonic speeds, k < 1. With these simplifications, 
the basic equations become 

u(Ou/Ox) + v(Ov/Ov) = (—kR/p,.) (Op/Ox) (la) 
u(Ov/Ox) + v(Ov Oy) = (—k/p.) (OP Oy) (2 


1=0 (Ou/Ox) + (Ov/dOy) = O (3a 


[O(ux) Ox] + [O(vx)/Ov] = 0 (3b 


These differential equations must now be solved to 
yield solutions satisfying appropriate boundary con- 
ditions. These are as follows: 

The tangency condition at the body surface is 


7 = 0 v= O (6 
The shock wave conditions at a point A are: 
Us = G4 cos (6 — a) (7) 


V4 = —q, sin (0 — a) 
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STAGNATION POINT OF A BLU 
1, = V (ku, sin B)* + (u, cos 8B (9 
Db,= dD. + pu (1 — Rk) sin? B (10 
Fig. 2 identifies the quantities appearing in these 


equations as follows. 

a is flow deflection across the shock wave at A, @ is 
the local shock wave angle at A and g, is the magnitude 
of the velocity vector after the shock wave at A. 


Eqs. (7) and (S) result from considerations of the 
geometry of the present coordinate system. 

Eq. (9) is derived from the conditions of conservation 
of mass and conservation of momentum along the 
shock wave direction, while Eq. (10) represents the 
conservation of momentum normal to the shock wave. 
Another condition namely, the stagnation point 
condition—is simply 

r= a y= @ u=v=0 (11 


Thus, the problem is to seek the solutions of the differ- 
the above 
for 


ential Eqs. (la), (2a), (3a), or (3b) with 


boundary conditions. This will be carried out 


both the two-dimensional and the axially symmetrical 
flows. 


Two-Dimensional Flow 


Near the stagnation point, it can be assumed that 


u = xf'(y) = x(df/dy) (12) 
Eq. (8a), the continuity equation, is satisfied if 
v= —f(y) (13) 


rhe choice of the particular « form satisfies Eq. (11). 
It can be justified by the following argument. In 
general, u (x, y) may be expressed for small x as an 
ascending series of the form 


u(x, ¥) = do + QQx + 2x7 + 


where the a's are functions of y. Clearly the sym- 
metry of the problem permits only terms in odd powers 
near the stagnation point 

There- 


of x, and for small x’s—.e., 
third and higher powers of x are unimportant. 
fore, Eq. (12) represents a satisfactory approximation. 

‘The quantity p can be eliminated from Eqs. (la) and 


2a) by cross differentiation. The result is 


u(Ow Ox) + v(Ow OV) = O (14) 
for two-dimensional flow. The vorticity, w, is ex- 
pressed by 

w = (Ov Ox) — (Onu/Oy) 


Physically Eq. (14) means that w remains constant on 


a streamline. Eqs. (12), (13), and (14) can be com- 
bined to yield the total differential equation for f (y) 
namely, 


(15) 


Phis third-order ordinary differential equation has the 
general solution 


f(y) = Ayexp (Avy) + Aiexp (—Azy) (16) 


N 
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where Ao, A; and A» are integration constants to be 


determined from the boundary conditions. The tan- 
gency condition on the body surface, y = 0, yields 
Ay = —A,. Thus, Eq. (16) becomes 

—v = f(y) = 2Ao sinh Aoy (17 
From Eq. (12), 

u = xf’ = 2ApAox cosh Aoy (18) 
The expression for w is 

w = —xf” —2A,A,*x sinh Asy (19) 


The boundary conditions at the shock wave are 


2A, sinh 4.6 = V (ku. sin 8B)? + (u_ cos B)? X 

sin (9 — a) (20 
2A,Aox cosh Ax = V (ku_ sin 8)? + (u_ cos 8B)? X 

cos (6 — a) (2] 
QA ,Ao*x sinh Ad = 

—u_*(1 — k)? sin B cos B(de dx 
(22 
V (ku_ sin 8B)? + (u. cos 8)" sin (@ — a 


Eqs. (20) and (21) follow directly from Eqs. (7) and 


(8S). Eg. (22) can be established easily by using Eqs. 
(la) and (2a). Near the stagnation point consider 


x, — 0, then A — C where C is the point behind the 
normal shock on the axis of symmetry (see Fig. 2). 
Obviously, when the shock wave conditions are applied 


near C, the following limiting values must be taken: 


x20, av2aO, 6~7~ 272, Bor 2, 
cos (@ — a) ~>0, cos8—~>O0, sinBml, 
sin (@ — a)—~>1, 6—~ 6* 


It can easily be shown that to the present order of 


accuracy 
. cos®s dg l 
lim =|{—- ee ve 
er dx] .-+0 Rs 
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i cos (@ — a) ( d6 L **) l 4 k— 1] 

1m = = ss = 

x0 x dx dx x0 R, kRs 


where Rs and R, are the radii of curvature of shock 
and body at the body axis.—| When these considerations 
are applied to the conditions in Eqs. (20)—(22), it is 


found that 


2A, sinh A.d* = ku 20a) 
2A, Az cosh A.d* = — (u,/Rs) (21a) 
; k 
2A, Ap” sinh A.d* = (1 — k)? a — 294) 
( s)” 


These equations can be solved simultaneously for 
Ao, A, and 6*. The results are: 
Ay = ku, VR(2 — k)/2(1 — k) (23) 
A, = —(1 — k)/kRs (24) 
6* = [(—kRs)/(1 — k)] sinh—'! [(1 — k)/V R(2 — k)] 
(25) 
Substitution of the results in Eqs. (17) and (18) yields 
the following velocity profiles in the shock layer: 
u/u., = (—x/Rs)Vk(2 — k) cosh ) (y/6*) sinh! 
(1 — k)/VR(2 — k)]} (26) 
v/u. = —k[VR(2 — k)/(1 — k)] sinh 
) (y/6*) sinh-' [(1 — kR) Vk(2 — k)]} (27) 
Hayes’ obtained the results below by a different ap- 
proach: 
6* = —(kR/V1 — 2k) sinh! V(1 — 2k) /2k > 
kRR/(1 — k)] sinh [1 — k)/V 2k] (28) 





u/u,, = —(x/R) V 2k cosh [(y/6*) X 
cosh! (1/2k)] ~ —(x/R) V 2k X 
cosh } (y/6*) sinh—! [(1 — k)/VW2k]} = (29) 
v/u,, = —k V 2k/(1 — 2k) sinh [(y/6*) sinh! 
V (1 — 2k)/2k] ~ —[kV 2k/(1 — k)] X 
sinh } (y/6*) sinh~! [(1 — k) V 2k]! (30) 


Obviously, these results are in agreement with the 
present results as far as the leading terms containing k 
are concerned. 
The expression for surface pressure coefficient is 
found to be 
C, = (2 — k) [1 — (x/Rs)?] (31) 


The general expression for vorticity is 


Rs) Vk(2 — k) sinh 


((1 — k) VkR(2 — k)] sinh } (y/6*) sinh~! 


(1 — k)/VR(2 — k)]} (32) 


w= (%/6*) (1 


which shows that the vorticity vanishes on the surface 
of the body. 


) 


+ In the present computations, take Rp ~ Rg = R 
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Axially Symmetric Flow 
Near the stagnation point the w variation is assumed 
to be as in Eq. (12). The continuity equation, Eq 


(3b), is satisfied if 


v= —2f(y) 


When / is eliminated from Eqs. (la) and (2a) for 


axially symmetric flow, the resulting equation is 
u (Ow/Ox) + v(Ow/Oy) = (u/x)w 34 


This equation is the well known law of variation of 
w along stream tubes. Combination of Eqs. (12), 
(31), and (34) yields simply 


Eq. (35) in turn yields the general solution 

fly) = G+ Gy t+ Gy? 36 
in which the integration constants Co, C), Co, must be 
determined from the boundary conditions. The tan 
gency condition on the body surface implies that Co = 0 
The shock wave conditions considered near the stag 
nation point—1.e., x — O-—yield the following equa 


tions: 
C,6* + C.6* = (ku, /2 37 
Ci + 20.6* = —(u,/Rs 38 
C. = [(1 — k)?/2] [ku /(RRs)?] (39 





By substituting C, from Eq. (39) into Eqs. (37) and 
(38) and by solving for C; and 6*, the following equa- 


tions result: 
CQ, = * V1 — (1 — Rk)? (u. Rs (40 
6* = {[-1 + V1 — (1 — &)2]/[(1 — 2] RRs (41) 


uw, vcan now be expressed as follows: 





Uu Stas Bek Dee —{1 — (1 — k)?] 
u 5* (} ~— ») 
2 F2V 1 — (1 — &)? — (1 — &)* 3 19 
(1 — k)? 0” x 
% Z| #VI- Cl ky? — [1 — (1 — k)*] 
ul 6* it ae 
34901 — (1-8) - (1-8 ¥] 
2(1 — k)? 6* | , 


There is the question of the choice of the sign before 
the radicals in Eqs. (40) to (43). It can be seen that if 
the ‘“‘bottom”’ signs are chosen, Eq. (43) will give a 
value of v which undergoes a reversal in sign in the 
range 0 < (y/6*) <1. The physical meaning of such 
a solution is that there exists a flow reversal in the 
shock layer near the stagnation point. Since an in 
viscid theory should not permit a solution with flow 
reversal, the ‘‘top’’ sign in Eqs. (38) to (43) should be 
adopted. A comparison with Hayes’ work shows 
again that for k < 1, the results are in agreement. 
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STAGNATION POINT OF A 
The surface pressure coefficient is 
C, = (2 — k) [1 — (x/Rs)?] (31) 


which is the same result as obtained in the two-dimen- 
The vorticity on the body surface does not 
It can be 


sional case. 
vanish as in the two-dimensional problem. 


expressed as 


w = [(—u,)/Rs] [-—1 + Vk(2 — k)] (x/6*) (44) 


A comparison of the two-dimensional and_three- 
dimensional shock detachment distances according to 


the present study is presented in Fig. 3. 


RESULTS 


Since primarily three-dimensional shapes are of 
practical interest, the ensuing discussion and com- 
parisons with experiment treat only the three-dimen- 
sional case. 

A comparison of experimental values of pressure 
coefficients on a sphere with those computed in this 
analysis is presented in Fig. 4 for Mach Numbers 
1.97 and 5.8. The agreement is surprisingly good in 
view of the simplifying assumptions made in the theory. 
In particular, the values of the density ratio, k, for the 
test data—0.3881 and 0.192, 
tainly not very small compared to 1.0 as was assumed 
in the theory, yet the correlation is quite good. More- 
over, the agreement persists well beyond the theoretical 
As would 


respectively—are cer- 


limits of applicabilityt of the present theory. 
be expected, the agreement is better at the higher 
Mach Number (smaller k). 

Figs. 5a and 5b present a comparison of theoretical 
and experimental values of shock detachment dis- 
tance for spheres in air and in monatomic gases, respec- 
tively. In both 
theory and test improves with increasing Mach Num- 
The density ratio, k, in the tests 


cases the agreement between the 
ber (decreasing k). 
of Figs. 5a and 5b ranges from 0.24 to 0.54 so that 
the small density ratio assumption of the theory does 
not, strictly speaking, apply for these data. In view 
of this, the difference between theory and the experi 
mental measurements is not considered excessive. 

Recent experiments on hemispheres reported by 
Korobkin’ support the assumption that the uw velocity 
component is directly proportional to distance along 
the surface of the body. Fig. 6 shows the nondimen- 
sionalized velocities measured by Korobkin at four 
different Mach Numbers in the NOL Aeroballistics 
Tunnel No. 1. Within the limits of the data, it is 
seen that the local velocity on the surface is indeed 
linear up to 45 degrees and deviates from linearity only 
slightly thereafter. 

The velocity gradient at the stagnation point is a 
quantity of considerable importance in the blunt body 
heating problem. According to the present analysis 
the quantity is given by 


t+ See Discussion, page 31, on the limits of applicability of 


the theory. 
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This estimate is compared with Korobkin’s experimental 
measurements in Fig. 7. The agreement is good, the 
theory being consistently high by about nine per cent 
for Mach Numbers ranging from 2.0 to 5.0. Also shown 


in Fig. 7 are values of 
Ou ; 
has for y = 1.4) 
Ox =0 
at Mach Number infinity as predicted by the present 
analysis and by Lees.’ The latter two values are 
practically identical. 

The theoretical distribution of velocity and_pres- 
sures along the axial streamline between the shock and 
the stagnation point are presented in Figs. S and 9, 
respectively. No experimental data are available for 
comparison. 

In an analysis of the simple shear flow past a flat 
plate in a viscous fluid,® the senior author of the present 
paper has shown that ‘“‘vorticity number,” rather than 
the vorticity itself, is an important parameter in study- 
ing the effect of vorticity upon the boundary-layer 
behavior. The vorticity number is essentially the 
ratio of the vorticity in the free stream to the average 


vorticity in the boundary layer. It is expressed as 
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STAGNATION 


WwW = ( Ov Ox } 
A= 


— (Ou/Oy) 
boundary-layer thickness near the stagnation 
» 2 > 
point 


The vorticity number for spheres has been calculated 
by the present theory (Fig. 10). Because local velocity 
on the surface of the spheres varies linearly with x, the 
The 
and 


vorticity number is constant along the surface. 

vorticity number is a function of density ratio 
Reynolds Number (based on free-stream velocity, vis- 
It is 


seen that the vorticity effects increase as Mach Number 


cosity in the shock layer, and sphere diameter). 
increases. However, on spheres of practical sizes, ex- 
cept in extremely low density flows, the magnitude of 
£ remains small. It can therefore be anticipated that 
the vorticity effect upon the boundary-layer behavior 
in a blunt-nosed body may also be small. 

An estimate of the ratio of boundary-layer thickness 
to shock layer thickness near the stagnation point 
expressed in terms of the density ratio and the vorticity 
number may be obtained from Eq. (46). This ratio is 


given by 


A/é = Of& V 2k 





It is less than the vorticity number, which was shown 
above to be small for spheres of practical sizes in all 
but extremely low density hypersonic flows. Hence 
it appears that the boundary-layer thickness is neg- 
ligible near the nose of the sphere. 


DISCUSSION 


Two restrictive assumptions made in the present 
inviscid analysis are (1) rx ~ 1, and (2) k = p/p, 
constant. 

The assumption 7 x ~ 1 affects only the continuity 
equation which, for constant curvature and constant 


density, 1s 


uor Ox) + (rou/Ox) + (1 + Ay) X 


| (vor Oy) + (rdv Ov)} + Avr = 0 
Now r= [1 + Ay] Ry sin (x/ R,) 17 ) 
hence or Ox = [1 + Ay] cos ¢ (48) 
and Or Oy = sin ¢ 19) 


To say r = x implies that Ov Oy ~ 0 or sin ¢ < 1 and 
that Or Ox ~ 1 or Ay = O[6 R, | = OR] <1. 

Thus, the assumption 7 ~ x requires k < 1 (which is 
consistent with the explicit use of that assumption else- 
where in the theory) except on or very near the body 
surface. Moreover the assumption restricts applica- 
tion of the theory to the immediate vicinity of the 
stagnation point. 

Sufficient information concerning high-speed flows 
about spheres is available to permit an estimate of the 
extent of the subsonic flow region and, hence, an assess- 
ment of the validity of the constant density assump- 
tion. 

The between the modified 


excellent agreement 
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Newtonian hypersonic approximationt and experi 


mental pressure measurements on spheres has been 
pointed out by a number of authors—e.g., Lees,’ 
compressible flow, application of 


readily yields the Mach Number 


Oliver.» Assuming 
the 


distribution ou the sphere surface. 


Bernoulli's law 
The locations of 


modified Newtonian approximation referred to here 


+ The 


alters the simple Newtonian relation for pressure coefficient by 
replacing the approximate stagnation pressure coefficient of 2.0 
New 


tonian impact theory for application to hypersonic blunt bodies 


with its exact value. This type of modification of the 


has been suggested by a number of investigators on the basis of 


their experimental measurements 
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the local Mach Number 0.30 and 1.00 points on the 
sphere were thus calculated for the free-stream Mach 
Number ranging from 1.5 to 10.0 (y = 1.4). The re- 
sults of the computations are shown in Fig. 11. 


Fig. 12 shows a representative contour of constant 


total velocity for k = 0.20. This corresponds to flow 
past a sphere in ideal air [y = 1.4] with a free-stream 
Mach Number of 5.1. If one considers the three- 
dimensional Prandtl-Glauert factor, 1/(1 — Jf), a 


fair criterion of the significance of compressibility 
effects and permits about 9 per cent error, the local 
Mach Number 0.30 line represents the limit of validity 
of the constant density assumption for the calculated 
example. 

In view of the relatively severe limitations imposed 
upon the applicability of the present theory by the 
assumptions r ~ x and k constant the excellent agree- 
ment with experimental pressure coefficients well be- 
yond the theoretical limits of validity of the theory 
would seem fortuitous. It is conjectured that perhaps 
the errors incurred by the various simplifying assump- 
tionst are compensating on the sphere’s surface. 

It is interesting that the present theory gives pressure 
coefficients very similar to those given by the Newtonian 
hypersonic approximation. The present theory gives 
C, = (2 — k) (1 — ¢°), while the simple Newtonian 
value is C, = 2 cos’ ¢. Within the theoretical range 
of applicability of the present theory ¢ = sin @ and 
hence the pressure coefficient may be written C, = 
(2 — k) cos? @ The quantity (2 — k), which is the 
pressure coefficient at the stagnation point as given by 
the present theory, is a density ratio correction to the 
Newtonian approximation? resulting from consider- 
ation of boundary conditions immediately behind the 
shock. The present theory has also the apparent 


1 r~x,k = constant, and k< | 


~ See footnote on page 31. 
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advantage that it describes the flow field between the 
shock and body near the stagnation point while the 
Newtonian theory yields only surface pressures. 

It is noted that the development of the theory of this 
paper does not depend upon an explicit energy rela 
tion. This is an advantage that assumes considerable 
importance for the hypersonic flow of a real gas about a 
blunt body. In the present formulation, the theory 
can be applied to the case of a dissociated or 10n1zed 
gas when the value of k is correctly estimated, pro 
vided the constant density assumption is not invali 
dated by real gas phenomena. It should be noted that 
the theory assumes thermodynamic equilibrium; hence 
time 


some real gas phenomena, such as relaxation 


effects, are not accounted for. 

The hypersonic blunt body theory developed herein 
has the advantage of simplicity but its extent of appli 
cability is limited to the vicinity of the stagnation point. 
Steps have been initiated to extend the range of ap 
plicability of the theory by removing some of the re 
strictive assumptions. 

To provide an understanding of the curvature effects, 
an attempt is being made to treat the problem with 
the terms containing A retained in the basic differential 
equations. The task is being carried out for the case 
of axially symmetric flows for bodies with a constant 
K. In this approach it is assumed thatr ~ (1 + Ay)a. 
A preliminary examination of the results obtained to 
date appears to indicate that the effects of curvature 
termis are: 

(1) To increase the 
Hence, it is anticipated that the agreement between 
theory and experiment will be improved. 

(2) To provide a pressure relief on the body surface. 


shock detachment distance. 


This is in qualitative agreement with the centrifugal 


force correction to the classical Newtonian impact 
theory. 
The details of this analysis are presently being 


worked upon. It is anticipated that specific numeri 
cal calculations will be made to illustrate quantitatively 


the effects of curvature. 
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Unsteady Flow Around Thin Wings at High 
Mach Numbers’ 


MARTEN T. LANDAHL* 
Massachusetts Institute of Technology 


SUMMARY 


Che perturbation velocity potential for unsteady flow around 


thin wings is expanded in a series of .\/~?, taking into account 
linear and second-order terms in perturbation velocities. The 
first terms for the pressure are seen to be in agreement with 
those predicted by the Hayes-Lighthill piston theory. Addi 
tional terms provided by the present theory give improved ac 
curacy for moderately high Mach Numbers and also take into 
account some effects of three-dimensionality, such as moderate 
leading-edge sweep and spanwise downwash variations. Com- 
parison with previous solutions suggests that the present theory 
is sufficiently accurate for practical purposes down to J = 2 \ 
calculation performed for the simple case of two-dimensiona! 
bending-torsion flutter shows good agreement with results ob 


tained from the exact linearized theory for Mach Numbers 


abov € ] Oo 
SYMBOLS 
= distance from mid-chord line to elastic 
axis (positive aft), divided by > 
= semichord at wing root 
( = pressure coefficient, (p — p 1/2) pl? 
Cc; = lift coefficient, lift/(1/2)pl2S 
rv, ¥) function giving the shape of the airfoil 
C9, & = function giving the instantaneous location 
of the (upper) wing surface 
=v¥ 1 
k = reduced frequency of oscillation wh / l 
L», , WJ, = flutter coefficients as defined in reference 5 
7 = free-stream Mach Number 
= local static pressure 
= free-stream static pressure 
” = radius of gyration in pitch, divided by > 
= Laplace transform variable 
S = wing area 
= dimensionless time obtained by multi- 
plying by U’/b 
free-stream velocity 
= downwash on wing in linearized form 
oe = dimensionless Cartesian coordinates ob- 
tained by dividing by } 
= XH = location of pitching axis 
7 = distance from elastic axis to center of 
gravity (positive aft), divided by } 
a = angle of attack 
= ratio of specific heats 
6 = the larger of the quantities « and A 
€ = thickness ratio of wing 
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r = «imensionl im} litude of time-deper 
ent wing motion 

m = ratio of wing mass per unit 5] to th 
quantity 4o/ 


= free-stream density 


p 

Q?=¢rty dimensionless perturbation velocity 
tential 

¢ = first-order potential (in 6 

¥ = second-order potential 

w = circular frequency of oscillation 

Why Wer = circular frequencies of uncoupled vibration 


in bending and torsion, respectively 


INTRODUCTION 


PAPER by Lighthill,! Hayes’’ piston theory 


ly A 
concept was used for calculating the unsteady flow 


around thin airfoils at high Mach Numbers. Very 
simple results were obtained which included the effect 
of nonlinearity. The second-order term was given 


exactly in the limit of high Mach Numbers, whereas 
the third-order term, which is influenced by the change 
in entropy, was only approximately given. 

According to Goldsworthy* the relative error in 
piston theory is of order \/ In Lighthill's result the 
linear term for the pressure coefficient is of order 6.\/ 
and the second-order term of order 6°. Consequently, 
the dominating error term must be of order 6.1/~*, since 
6. is considered small. Naturally, there is no point 
in retaining the second-order term in Lighthill's ex 
pression unless the error is substantially smaller than 
this term. It follows from the above that this occurs 
when M/* > 6—'. 
M* > 20, which clearly indicates the limited accuracy 
In the 


present paper, the solution is therefore expanded in 


Thus, a 6 of 5 per cent requires 


of piston theory at moderate Mach Numbers. 
2, and a few more terms are included than 


powers of \/ 
the leading terms constituting piston theory 


METHOD OF SOLUTION 


We consider a thin wing, situated mainly along the 
xy-plane. The undisturbed stream flows along the 
positive x-axis. Nondimensional coordinates x, y, 2, ¢ 
All lengths are made dimensionless by divid- 
The time is 


are used. 
ing by the semiroot chord 6 of the wing. 
made dimensionless by multiplying by UL’), where 
U is the free-stream velocity. Then we can set U’ = 
1 without loss of generality, and, consequently, the 
velocity of sound is 1/7. To second order in per- 
turbation velocities the flow is isentropic and irrota- 


tional if viscosity is neglected. Then the perturbation 
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velocity potential exists and _ satisfies the following 
differential equation, omitting third-order terms: 
(1 — M?)d,, + dy + ¢o22 — 2M°b,, — MP?bn = 
M?(y — 1) (bp + oe) (Ger + yy + O22) + 
M?[(0/Ox) + (0/Ot)] (6,7 + 6,7 + ¢:7) (1) 
The position of the upper wing surface is, at every 
instant, given by 


z = bh(x, y, t) (2) 


Then we get the following boundary condition on the 
upper wing surface: 
o.(x, v, dh, t) = (1 + ¢,)6h, + $,6h, + 6h, (3) 
It is convenient to impose a condition at s = +0 in- 
stead. Using a Taylor series expansion and neglecting 
terms of third or higher order, we get 
o.(+0) = 6(h, + h,) + b[h,6, — hd.: + h,d,]|,= +0 (4) 
We also require ¢ and its derivatives to vanish far 
upstream. For the pressure we have, in isentropic 
flow, 
Cp = (2/yM*) G1 — [vy — 1)/2)IP(2¢, + 
26, + $:7 + by? + o7)$"7" — 1) (5) 
or, expanding and retaining linear and quadratic terms, 
C, = —2(¢, + o:) + Mo. + ¢,)? — 
(o,° + o,° + $:*) (6) 
On the wing surface the pressure can be expressed 
using the values of ¢ and its derivatives on z = +0, as 
follows: 
Cp = —([2(¢, + o) — M*(¢, + 6,)? + 
or? + o,? + o.7 + 25h(G,. + G:2)|,=40 (7) 
or, making partial use of Eq. (4), 
C, = —([2(¢, + o.) — M*(¢, + ¢,)? + 
o;? + ¢,7|,-+40 — 267h[(0/Ox) + (0/OF) 7h — 
6-(h, + h,)? (8) 
To solve the nonlinear equation for ¢ we use an iter- 
ation procedure. Putting 
oO ¢ + y (9) 
where ¢ = 0(6) and y = 0(6 
we get from the above equations 
(l — M?) ¢,, + Suy + g -_ 
2M*y,, — M*e, = 0 (10) 
(1 — M2)bre + Wy + vee — 2M, — Mn = 
M?(¥ — 1) (¢, +- C1) (rz + Cyy + 22) + 
M2[(0/dx) + (0/df)] (¢2 + ¢,2 + ¢2) (11) 
¢:(+0) = 6(h, + h,) = dwo(x, y, t) (12) 


¥.(+0) = 6(h-¢, — he. + hyey)e=+0 (13) 


Cy = [2(¢, ae 2) a 2(, + v1) oat 
M?(¢, o ¢1)? +- ¢,” + oy ls ate 
26*h[(0/Ox) + (0/Ot) |*h — 62(h, + h,)? (14) 
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By use of an iteration the first term on the right-hand 
side of Eq. (11) can be somewhat simplified, so that 


C~ Pk, +4, +4, ~ 20. ~ Wa © 
MA(y wa 1) (Or ¢,) [(O Ox) + (Oo Ot) |"¢ + 
M?[(0/Ox) + (0/Ot)] (¢2? + gy? + ¢22) (15 


To obtain an expansion for high Mach Number we put 


= Me (16) 


¢ = (6 M) lea + A/ IP) go + / MA gs] + 
0(6 AT") (1S 





Y = 6" IW + (| AI?) Ys | 1+ ((6? 1/4) (19) 


This gives the following sequence of differential equa 
tions to solve, with their respective boundary condi 
tions on zs; = +0: 


a T Yj = J (20 
1 = ray (21 
— 2 + Snaz = —}[(Oo Ox}) 
(0/0f;) |? + (07/0917) ; ¢, (22 
Ynz, = O n> 1 (23 
—Y + Visas, = (¥ — Deisn¢: = 
2¢ Fi, (24 
Vie = —hyis,- (25) 
— Por e Y20,2 ~ (0 Ox) [(y ~~, I Fla F 24 T 
2¢12¢221 — 1 [(0/Ox1) — (0/04) |? + 
(O- Ov?) } i + (0/0Ox) (en, — 1,,)° + C1, -] (26 
Yo. — —hgo.- + (h, = h, (1, — OY + 
hy oe, (27 


For the solution Van Dyke's method! can be directly 
applied. First, a smoothing is carried out which re 
moves the jump in velocities across the shock. This 
is achieved by considering a small flexible part attached 
to the nose of the airfoil so that the relative wind is 
always parallel to the surface of the nose. The smooth 
ing technique is discussed in detail in reference +. Then 
we apply a Laplace transformation in the .\-direction. 


We set 


Llgn) = on = ce" go Gk; (28 

~~“) 

. 
L(y,) = J, = | ey, dx (29 

~/Z 
The solution for ¢, is easily found to be, for z, > 0, 
2} = —(1/s)®(s, the” (30 

where 
@o(s, fy) = [ e W(X, Vi, ty + X)) dx, (31 
~~ 
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Inversion then gives, by use of standard formulas, 


a 


Che lower limit can, if desired, be replaced by x,, 
where x,(y) gives the position of the leading edge. 


For 2: we get the following differential equation: 
be + $os:: ils — (00h) |? + 
(0? Oy,2){ Zils, #O)e (33 
The solution is obtained by putting 
= Ax + Biz" 34 


Inserting this into Eq. (33) and using the boundary 


condition, Eq. (25), gives 


1, = (1/2s*) | [s — (0/0h) }? + (0? Om"); Zils, +0 
B, = (1/2s) | [s — (00h) ]? + (0?/O917)} Zils, +0) 
O Q\* oO? | 
o~ Fea eek 
z LOX) Of; OV: 
(Xy Z fi )orl£1, Vi, +O, ft) d&) + 
0 oO\*? oO ; 
ta[(2-2) +h) e-2- 
2 Ox, = Of; Ovi" Jj « 
wu(E, Mi, + i) dh (35) 
3y using the result from Eq. (32 
O oY. 
g@=- 2 
tL \Ox, Of; Ovi" 
} (Vy 21 &) ty (Ey \ y+? & d§ 
( 0 7] ) | [ 
~ _ + vy S £1) X 
? , OX) Of; OV; jv x, 
Wi (1, Vi, + & d& (36) 


Laplace-transforming the iteration equation for ¢3 gives 


3 + Sane, = —tls — (00h) |? + 
‘ 


0? Ow*)fe "Gels, +0) — sai} ls — (0/04) ]? + 
(0? Oy,7),e "Ze(s, +0) (37) 
The solution is found by putting 
C3 = Ax a a B3z\e + ¢ 1 7¢ (3S 


By substitution, the coefficients are easily found to be 


1; = (3/48) | [s — (0/04) |? + 
(02 dy,2)!o(s, +0) (39 


B, = (34s) | [s — (0 Oh) |? + (02/dy,2)! go(s, +0) (40 
C 1/4) j [s — (0/0) ]? + 0? On?} Go(s, +0) (41) 


In particular, we find by inverting and making use of 
the above results for ¢, 
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For the nonlinear part of the solution we proceed in 
about the same way and get, after Laplace-transform 
ing Eq. (24), 


oF 7 » 
yl Y loz - - 


The solution is obtained, as before, by a proper com 


bination of s;e °” and ¢ rhe final result, taking 


) 


into account the boundary condition, Eq. (25), reads 


For yo, the Laplace-transformed right-hand side of 
Eq. (26) will include terms proportional to e~“, as 
well as terms containing 2, rhe solution is then 
a combination of 31 ~¢ » Sy and After some 
manipulation we get, for the x)-derivative on the plane 


s; = +0, 


‘(3 >) a) lf asd ; 
t| \ox, Or, dv:2 |. eS ee ree es 


2hes.., — 2(h h ¢ C1 - 2h, ¢, (45) 
Having thus obtained the approximate solutions for the 
velocity potential, we expand the expression for C, in 
the same way, 

C, = (6 WM) [Cr + AIRC, + 
(1/M‘)C, + d(C, + Ci? ] 4 
O(6 17) + O16? AP 16) 


~ 


From Eqs. (8), (18), and (19) we get, using the vari 


/ 


ables /;, x), and 2), 


C, —2¢ (47) 
Cn” = 2fi.,+ ¢ u - 2hwo, (48) 
- = 2Yo, win 2¢ 2 ” ae | - ¢£ 
C1, . (49) 
where all derivatives are to be evaluated at 2, +0. 


Inserting the results for the expanded potential gives, 
finally, after a return to the original coordinates, 


CaM = 2w (50) 
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Cr = pf (x — £&)wo(&) dé (51) 


I 


l . ‘es 
5 (37 — 5p | (x — E)wy?(E) dé + 
F *L 

| - . — | o 

i (y + 1) wy(E)D h(E) dé — i (y + 1) > x 


ad 


ey | 
{ wy(£) dé + 5 (y + l)w? — h,? — h,? - 
J Xy . 
1 I Pe 
i (y + 1) [A — A(xz) |Dih(x,) — 5 Dih?(x1) + 


hDyh(xz,) — 2hyh(xzp) + 2h,h,(xy) (54) 


where D = (0?/Ox?) + (07/dy") 
dD, = (07/Of?) + (07 Oy’) 
W(E) = Wl(E, v, f — xX + &) 
h(é) = h(é, y,t — x + &) 
h(x,) = h(xz, y, t — x + Xz) 


and the argument of wy) and / is (x, y, ¢) if not otherwise 
stated. After the integrations are carried through the 


smoothing can be dropped. 


COMPARISONS WITH OTHER RESULTS 


The above result can be subject to a number of 
checks. First, it is seen that the first linear term 
C,: and the first second-order term C,,‘”) are, indeed, 
identical with those obtained from piston theory. 

For an infinite swept flat plate, the present theory 


gives the following lift curve slope: 


Cra = (4/M) + (2/M*) (1 + tan? A) + 
(3/2M*) (1 + tan? A)? (55) 


where \ is the sweep angle. Eq. (55) is then seen to 
produce the correct expansion in powers of \/~? of 
linearized swept-wing theory, which gives, for super- 
sonic edges, 

Cr, = 4/V M* — (1 + tan? A) (56) 
At JJ = 2 and A = 40° the error in the approximate 
formula, Eq. (55), is only about 3 per cent. 

For a two-dimensional flat plate oscillating in ver- 
tical translation, the force coefficients L; and Ls, in the 
standard American notations,* become 
L, = [(1/2M*) + (3/8M*)] [(sin 2k)/Rk] + 

(3/4M*) cos 2k (57) 


(1/RM) + [(1/2M*) + (3/8M*)] X 
[(cos 2k)/k] — (3/4.M*) sin 2k (58) 


isi 
t 
II 





In Fig. 1 these results are compared with the exact 
linearized results by Garrick and Rubinow.’ At J/ = 


JANUARY, 1957 
2, the agreement is quite good for k < 1. Even at 
M = 5/3, which cannot be considered a very high 
Mach Number, the deviation from the exact theory is 
not too large for, say, k < 0.7. The present theory does 
not work for high reduced frequencies since the iter- 
ation introduces on the right-hand side of Eqs. (22) and 
(26) terms of the relative order k?.\/~*, which must be 
small for the iteration to work. This is, however, not 
a serious restriction from a flutter point of view since in 
most practical cases k < 1. 

For the nonlinear case, a symmetrical airfoil of thick- 
ness ratio ¢ oscillating in pitch about the axis x = x is 
considered. In the above formulas, then, 


bt 


6h = eg(x) + A(x — x)e™ (59 


The calculation of the pressure is quite straightforward 
(although somewhat lengthy for C,»), involving only 


relatively simple integrals like 


x 
/ ao ¢-ajc" “se 
70 


in addition to those like 


appearing in the linearized part. Terms proportional 
to e« and \* drop out when taking the pressure difference 
across the wing and can, therefore, be omitted in the 
calculation. 

Of particular interest is the coefficient k.1/;, which 
represents the damping moment in pitch. For low 
values of k, the result simplifies further, to give 


I Xo | 2 2 
ee 
3 f MoM 2M* 


l 
se(v tit 


Using this expression, the case of a 6 per cent thick 
biconvex airfoil with g = x(2 — x) is computed for two 
locations of the pitching axis—namely, at the leading 
edge (4) = 0) and at the mid-chord (x) = 1). Fig. 2 
shows that good agreement is obtained with Van Dyke's 
results’ for M > 2. 
the leading edge, one notices that the second-order 


For the pitching axis located at 


piston theory happens to give results closer to the exact 
second-order theory than the linearized theory does. 

For higher values of k, a single wedge oscillating about 
its leading edge is treated. For this case, reference 4 
gives the second-order lift and moment coefficients for 
arbitrary frequencies, the expressions being quite com- 
plicated. The present theory gives for the force co- 
efficients 
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k2(L3; + iL3) = (1/M) (1 + tk) + (2/4Rh.173) X 
e~** — 1) + (3/161) [5e7""" — 1 + (ik) X 
(l—e “")| + (e 4) (vy + 1) (A + ck) + 
(e SV?) [y — 3 + ikR(y + 1) + 

(i 2k) (By — 1) (e7*"* - 1 (61 

As seen in Fig. 3, where JJ = 2 and e = 0.10 are chosen, 

the agreement with the exact second-order result is 


quite good in the range of reduced frequencies shown 
(Rk = 0-0.5). 


The present method becomes particularly advan- 
tageous for flutter calculations. This is exemplified 
below for the case of supersonic bending-torsion flutter 
of a thin wing section. In the exact linearized treat 
ment by Garrick and Rubinow’ the flutter equations 
had to be solved by a trial-and-error method, since 
the reduced frequency and, hence, the aerodynamic 
forces, were not known beforehand. With the present 
expansion of the aerodynamic forces, however, the 
flutter speed can be obtained by iteration in essentially 
the same manner as the aerodynamic forces. First a 
solution correct within O(.1/ is sought, so that the 


first-order result in .1/ i.e., piston theory—is used 


for the aerodynamic forces. In the particular case 
treated, the following closed form expression is obtained 


for the flutter speed :' 


(uM)? jra?[1 + (@n/we)?] — 


U, X (wy/q)? + (ra? X) (Xa? = 7a)! ay 

— (Ol 
bw. \ pM(X, + a) — 

(uMaX/r,.") (@n/Wq)? — (1/3 
where 
y Yo + 2x,4 + a? + (1/3) 
; 1 + (1/77) (wn Wa)? [a? + (1/3) 
(bw. U; alll | Rk? (63) 


Using the value of & thus obtained, the correction terms 
of order .\/~* in the aerodynamic coefficients are calcu- 
lated and thence new improved values of Uy and Rk. 
This calculation, which is not reproduced here, only 


involves evaluation of a single radical as in Eqs. (62 


and (63) above. 

The iteration can then be continued by using the im- 
proved value of k and adding the correction term of 
order \/ No tables for the aerodynamic coeffi- 
cients are required, since these are calculated simul- 
taneously with the flutter speed. For the particular 
case chosen in Fig. +, it is seen that one iteration already 
gives results very close to those of Garrick and Rubi- 
now,’ for J > 1.5. For other combinations of the 
involved parameters, however, two iterations might be 
required for the lower Mach Number region. 


r< 
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CONCLUSIONS 


The method presented here for treating the steady or 
unsteady flow around thin wings for high Mach Num 
characterized as an inverse Rayleigh 
Janzen method. The leading terms of the expression 
in .\/~* are identical with those given by Lighthill's 
piston theory. Additional terms in the present solu 


tion provide extra accuracy, so that it can be used for 


bers can be 


engineering purposes roughly down to ./ = 2. For {| 
very high Mach Numbers the solution is only valid 
entire 


when 6°.)/*7 < 1, i.e., it cannot be used in the 


hypersonic region. Because of the way the solution is 
obtained, it is also restricted to low reduced frequencies 
for which k?.\/ 1, although the linear piston theory 
happens to give the correct result for k — for all 
supersonic Mach Numbers, according to Hjelte.’ 

The three-dimensionality of the flow is taken into 
The method 


does, however, predict the effect of moderate leading 


account in a fairly primitive manner. 


edge sweep and smooth spanwise variation of down 
Wing tip effects or the effect of a discontinuous 


leading-edge slope must, however, be treated by a more 


wash. 


complicated method, as shown for the linearized case | 
in reference 6. Fortunately, the influence regions be 
come smaller with increasing Mach Number and, there 
fore, probably of lesser importance. ' 

The new procedure suggested for flutter calculations, 
which logically follows from the aerodynamic calcula 
tion procedure, represents a substantial saving of labor 
since both the usual guessing of reduced frequency 
and the need for tables of the aerodynamic coefficients — 
are eliminated. 

A particular flutter example shows the iteration to 
converge very rapidly so that very few terms usually 
need be retained in the expressions for the aerodynamic 
coefficients in order to get good accuracy. Actually, 
piston theory alone already seems to give fairly good 


results for .J/ > 2. 
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end BERTRAM KLEIN 
ed for | North American Aviation, Ince. 
For ' 
valid 
entire SUMMARY tually would be replaced by a matric solution done 
ion is A simple method of matric structural analysis is presented automatically on some kind of calculating machine 
encies which is believed to have certain distinct advantages over existing In 1947, Levy laid the groundwork for a matric method 
heorv methods as explained in the text. The possible disadvantages based on an energy si ution of a redundant structure. 
or all ae auaeeee were my ne mguengierte due: and scope This approach has been extended and streamlined by 
f automatic computing equipment are made larger. Certain : ; ; , . agua se 
‘ ; é various investigators. wee rererences ZZ, 2/4, 2Y, 32 
two- and three-dimensional problems are worked to illustrate the : . . ) ; i : ; 
1 into — ease and simplicity of the method 39, and 46.) An interesting discussion of this type of 
ethod | method and the Cross-Southwell (spring matrix) type 
ding SYMBOLS of method may be found in a recent series of articles by 
i — : - . , — i 
} ae , ' Argyris.** In 1950, Heldenfel: ‘t up in matrix form 
lown ' A = effective cross-sectional area of uniform axial ta : c ‘ - agate 
1u0Us ' dlineneit a numerical iterative procedure developed by Duberg. 
more A, B, H = submatrix Levy later advocated a solution by a growing unit spring 
case = shear panel dimensions matrix type of method.** Cox*’ and Williams‘! de 
s be ~ ee ee veloped finite difference methods. Various other 
E, G = Young's and shear moduli : . eee 
here allan tials matric methods have been and are being devised. (This 
/ = g dept! 
| } = twice the stiffness of an axial element short historical survey 1S by no means intended to be 
10ns, / = axial load in the direction of an axial element all-inclusive. 
cula = shear flow For several years now the author has been thinking 
labor = shear panel aspect ratio about some way of simplifying the setup procedure in 
= shear panel gage ‘ ras a : 
encv a : matric structural analysis without subtracting from the 
a ; (, £ = joint displacements . . : . 
ients = average shear angle of shear panel generality of the approach. As it turns out, there is 
= axial displacement in the direction of an axial such a possibility; but with the simplification there is a 
n to element penalty in that the size of the matrix is larger than in 
tally Subscripts other methods. Nevertheless, the author feels that as 
ns 2 2. ,?, = numbering of joints, loads, displacements, the field of automation expands and progresses, man 
ally, and equations ageable matric size also should grow with time. 
= = secant 
700d 
THE METHOD 
INTRODUCTION 
ib The first step in any method is the idealization of the 
ITH THE ADVENT of high-speed computing ma- : ' 
; : structure to be analyzed. Only certain semimonocoque 
F chines, the door has been opened for the prac- ; aie 
mover, ; ‘ : structures will be considered here, where it is assumed 
1()9 tical use of matric methods in structural analvsis. ‘ i aa 
S é i : that the structure is composed of uniform (untapered 
Heretofore, the amount of computational work asso- 
; ae structural elements, some taking only axial load and 
ath., ciated with these methods was considered prohibitive ; rs 
. : ; ‘ ; others taking only shear loads in the form of constant 
for hand solution. As is true in many fields, the im- ee 
a ; shear flows (except for triangular panels However 
petus has been given as soon as the calculation part pce 
\ppl : : other structures and idealizations where members are 
has been found to be realizable. However, the use of j 
: f tapered or take combined loads also may be treated 
1 matrices in structures is not new. Apparently, one of ; 
a , ; = »v the method. should be note at sometimes 
; “ by th thod. It should | ted that t 
94. the first to try these methods was Baver, back in 1927, ep: : ; 
ting :; ane ass ae . ; ac a joint of a structure must be considered to be a struc 
ome in connection with reinforced concrete design. There : 
[per- me : tural element. 
seems to have been a lull in the Thirties, although 
_ matric methods found their way into dynamic and elec- 
DC trical fields.'~* The aircraft companies and civil —§. 
« e ° @ c 
engineers supplied the interest of the structural people p™~ 
on for matrices in the Forties.’~*! The author began t / 


ae | applying matrices to structural problems in 1943. ata 


The spring matrix approach was pursued. He believed 


that Southwell's systematic relaxation procedure even- ito tr nat 
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The correct solution to a stress distribution problem 
for an idealized structure is the one in which all struc- 
tural elements are in equilibrium, all displacements 
are compatible with the corresponding forces, and all 
boundary conditions in terms of forces and/or dis- 
placements are satisfied. These conditions alone are 
sufficient for the 
sufficient only to write equations (a) of equilibrium for 
each structural element (when necessary), (b) of dis- 
placement-force compatibility for each structural ele- 


correct solution. Therefore it is 


ment (when necessary), and (c) for boundary cond1- 
tions (when necessary). 

For the problems considered here, the unknown 
quantities are: (a) axial loads acting at the ends of 
certain axial elements, (b) shear flows (or forces) acting 
along the edges of certain shear panels, and (c) certain 
displacements at the corners of these panels. The 
number of these unknowns must be 7dentica/ with the 
number of available equations. 

For the structures considered here, the nonvanishing 
numbers in the section of the matrix corresponding to 
the equilibrium equations are either 1 or —1, since all 
equations are colinear. (See 
The force-displacement 


forces entering such 
illustrative problems below.) 
equation for an axial element, such as shown in Fig. 1, 


is of the form 


(1/2) (P; + Pin) = 641 — 6,) (2,--4/a) (1) 


or 


(P; + Piss) = (6:41 — 6,)2(E,,.,A /a) (la) 
which states that the average force in the element is 
proportional to the elongation of the element. The 
force-displacement equation for a uniform rectangular 
shear panel, such as shown in Fig. 2(a), is written 


Goee t} 1/2)r[ (3 + us) — 


qa = Gi4t70 = 
(v1 + v3) — (v2 + )]f (2) 


(uw; + u2)] + (1/2) 





where the term in the braces represents the average 
shear angle in the panel (times a), and where r = a/b. 
For the right isosceles triangle shown in Fig. 2(b), the 
average shear angle (times @) is 


(2a) 


ay = —&& + U3 + % — V 
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As far as boundary considerations are concerned, in 
the case of flexible supports equations similar to Ey 
(1) may be written; for rigid supports it is a matter of 
omitting equations and unknowns. Wherever an ex 
ternal load is applied to a joint, that joint must be in 
equilibrium. 

Instead of trying to eliminate unknowns by combin 
ing Eqs. (1) and (2) with the equilibrium equations, the 
method proceeds by entering the original equations 
into a single square matrix and a RHS (right-hand sid 
matrix as shown schematically in Fig. 3. 

As mentioned previously, submatrix A, in the left 
hand corner of Fig. 3, contains only the numbers 1 or 
—1. The submatrices B and C are weakly populated, 
the submatrix D is a diagonal matrix, and the sub 
matrix / is heavily populated. The size of the RHS 
matrix depends on the number of loading conditions 
that are treated simultaneously. It will be noticed 
(see illustrative problems below) that each entry in the 
matrix is a single entry—i.e., it does not consist of sums 
of numbers. Furthermore, because of the simplicity 
of the equations, the entries can be written down very 
rapidly. The simultaneous solution for the unknown 
forces and displacements may be carried out by any 
convenient method such as matrix inversion or elimi 


nation on an automatic computing machine. Some 
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A SIMPLE METHOD 


times it may be of value to rearrange the entries in the 
matrix for purposes of solution. However, the arrange- 
ment shown in Fig. 3 is simplest for formulating the 
problem. To obtain influence coefficients, it is merely 
necessary to incorporate the correct RHS and include 


the necessary equations in the matrix proper. One 
feature of the method is the ease of treating changes 
in the structure. If an element is removed from the 
structure, or one added, it is merely necessary to delete 
or add certain equations and, consequently, numbers 
in the matrix; if the size of an element is altered, only 
one or two numbers need be changed. Since we are 
working with deflections and, consequently, strains, 
thermal load terms may be incorporated directly into 
the RHS matrix. Also, 


with forces and deflections, mixed boundary condi- 


since we are working both 


tions can be treated. 


In recapitulating, it may be stated, as seen by the 
foregoing development, that nowhere in the method is 
it necessary to classify certain forces as being redundant 
and make subsequent ‘“‘cuts’’ in the structure; nor is it 
necessary to write quadratic energy expressions in 
applying the method. No joints must be “locked,” nor 
It is not necessary to 


the 


elements combined in any way. 


free bodies of whole cross-sections of 


It is not necessary to employ the artifice of 


consider 
structure. 
dummy loads to obtain deflections at unloaded points. 


Several problems are worked below to illustrate the 
mechanics of the method. The structures considered 


are of the semimonocoque type discussed earlier. 


1) Two-Dimensional Shear Lag Problem 


Consider first the plane sheet-stringer-rib panel de- 
picted in F*; 4. For purposes of simplicity, the struc- 
ture is assumed to have an axis of symmetry, all axial 
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Solution: The unknowns listed below were obtained by 


matrix inversion on the IBM 701 electronic digital 
represent the 


computer. These results appear to 


physical behavior of this type of structure for the given 


’F MATR 
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members are assumed to be identical, all panels to be 
identical, the loading is taken to be symmetrical, and the 
structure is assumed to be completely fixed at the sup 
found to exist 
; and 8 de 


ports. The following unknowns are 


6 axial loads, P;—Ps; 2 shear flows, g; and 


flections, “;—-u,; and 7v;-a, one in each direction at the 
four corners of the lower left-hand shear panel. Since¢ 
these unknowns number 16, there must exist 16 equa 


tions These equations consist of 6 axial element 
equilibrium, S axial force-displacement, and 2 shear 
force-displacement, which may be written as 
axial element equilibrium: 
1) Py —- Pi + qa =0 D # P 0 
2) P,} — Px + ga =0 5) P ya 0) 
(3) -P; — qa 0 6) Pe + qia Jol 0 
axial force-displacement: 
(1 Po ig Ps = (4 — Us k 2 pP i i k 
2) Py + Po = urok 6) Pe = (3 v2)R 
3) P; = (us — us)R ri a” vk, 

k 2EA 
(4) Ps + Ps u4k S) P 6= vk 


and shear force-displacement 


1) (2,Gt)qia = (us + us) — (41 + 


(2) (2/Gt goa 


No boundary condition equations need be written in 


this problem, since such conditions have been satisfied 
automatically when the above equations were written 
Now let k = 2/A/a = 1, and choose A at so that 
2/Gt = 4(E£/G) (A/at) = 10. Let P l Then, 
according to the scheme presented in Fig. 3, the matric 
formulation of this problem may be written as 
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0 0 0 0 0 0 0 
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loading. Also shown are the check RHS terms calcu- 
lated by substituting the unknowns back into the 16 
equations. Total computing time was less than one 
minute. 
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P,; = 0.77377893 Wu, = 3.2030849 ¢, = —1.0000001 cy = —0.000000016 
P, = 0.65552699 ie = 1 .4293059 oc =  0.000000008 Cw = —0.000000003 
P; = 0.22622108 uz = 0.79691518 Cc; =  0.000000009 Cu =  0.000000002 
P, = 0.34447301 us =  0.57069409 Cs =  0.000000003 Co = 0.000000002 
P; = —0.22622109 “1 = 0.45244218 C; = —0.000000002 C3 = 0. 0O00000005 
Ps = 0.10796915 v% = —0.21593830 cs =  0.000000003 Cu =  0.000000001 
qa = —0.22622108 v; = 0.22622109 c; = —1.0000002 Gs =  0.000000027 
qa = —0.11825193 vy = —0.10796915 cs = —0.000000016 Cie = —O0.000000057 


Eight significant figures were carried in the computa- 
tion. The worst check, percentagewise, is that for Eq. 
7, which shows an error of 2 in 107. 

Actual values of deflections may be determined by 
multiplying the above values by the actual value of 


go =a 2EA. 


2) Two Spar Swept-Wing Analysis 

The next problem considered is the analysis of a 
three-dimensional structure—i.e., a two spar 45° 
swept wing with carry-through bay and loaded under 
symmetric bending as shown in Fig. 5. The wing is 
assumed to have two planes of symmetry, rib webs 
as spar webs, identical outer axial members, and iden- 
tical values of Gf for all shear carrying members. The 
attachment bulkhead at the wing-fuselage intersection 
is assumed to be infinitely rigid in its own plane, but 
free to warp out of its plane, which means that the 
only motion at this section is axial displacements of 
the carry-through spar caps. The applied bending 
moment causes an antisymmetrical state of stress 
about the horizontal plane of symmetry, so that only 
the forces and displacement in the top half need be 
considered as unknown. 

In this problem, the unknown quantities are 
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Fic. 5. Two spar swept wing. 


ay: ) 5 spar cap, Pi-P; 
6 axl ads - : 
} axial fonds ( 1 rib cap, Ps 


3 cover panel, gi-q; 
11 shear flows <, 5 spar web, gi-gs 


3 rib web, g9-gu1 


(The positive sign convention for shear flows are shown 
in Figs. 2 and 5.) 
3 each at 5 regular panel points, w#;—1;, 
o P V)}-Vs5, and W1-W5 
17 deflections mig ; : 
| 1 each at the 2 fuselage intersection 
points, 6, and 6;. 


(The positive sign convention for displacements is 
shown in Fig. 5.) 
The necessary 34 equations are comprised of 


© spar cap 
13 axial element equilibrium 3 rib cap 
5 web cap 
mae 2 j 7 spar cap 
10 axial force displacement SS |! 
| 3 rib cap 
3 cover panel (one 
, triangular) 
11 shear force displacement . - 
| 5) spar web 
3 rib web 


The triangular panel is assumed to have loads as 
shown in Fig. 2(b), and a shear-displacement equation 
as defined by Eq. 2(a).. The forces in the carry-through 
spar caps are determined from joint equilibrium to be 

joint 4: (1/72) (P3 — qa) 
joint 7: (1/Y 2) (Ps + Pe — gsa) 

For the given structure and loading, no shear exists 
in the carry-through bay cover panels or spar webs. 
Assume 2/A a = 1| and 2/Gt = 10 for the outer struc- 
ture; a = 0.2. The values of 2A 2a are taken to 
be the following for the carry-through spar caps: 

k, = 0.73858450 and ks = 0.59969722 


Also let Py = 1. 


The matric arrangement of the problem becomes that 


for the front and rear, respectively. 


given in Fig. 6. 

The matric inversion solution is listed below to 
gether with the check RHS terms. 
computed by the IBM 701 electronic digital computer 


These results were 
in about 6 minutes. Comparison of these values with 
known solutions of swept-wing problems of the type 
treated indicates agreement as to the nature of stresses 


and deflections. 
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Py 
P, 
P; 


P; 
P, 
quia 
qra 
q3a 
gaa 
G5 
Gea 
qia 
Gs 
gga 
qiod 
quad 


—O 
—(). 
—U 
—Q. 
0 
0 
0 
—. 
—0 
—Y 
—Q. 
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96860242 
76831339 
16004218 
0313978 
2316862 
13076514 
015699017 
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Matric tormulation of two span swept-wing proble1 


— 66 


—] 


a. 


3939165 
260453 1 
383929 
$125313S 
2118685 


9.366273 


OSS39SO 
97015668 
298437 
65059807 
3954034 
260453 1 
700904 
3640058 


2118685 


825361 
5569381 


TH 


; 


cq = —1.0000004 
c = —0.00000003 1 
oa = 0 . 000000040 
C; — —0.99999968 
% = 0). 000000277 
c& = 0. 000000032 
“a = 0. 000000015 
cz = —0.000000109 
cy = —0.000000026 
Cw = —0.000000494 
cu = 0. 000000396 
Cy = 0 . 000000057 
3 = 0. 000000177 
4 = 1 .OO000003 
Qs = —0.000000102 
Cs = 0. 000000036 
cy = 0.99999985 
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000001088 
QOO0000630 
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QOO0007 26 
000000459 
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t 
Eight significant figures were carried in the calcu- sidered in the previous structure, except for the values ~ 
lations, with the worst error in the check RHS being of 2FA_ 2a for the carry-through spar caps, which are, 
for Eq. 27, 1.222 in 10°. However, since w; and w», from forward to aft, t 
F 
which appear in this equation, can be accurate only to k, = 0.80109169 k» = 0.84206305 &, = 0.59764379 t 
6 decimal places, no greater accuracy can be expected it 
when eight significant figures are carried. The matric formulation of this problem, set up ac- T 
cording to the same simple procedure used to solve the 1 
3) Three Spar Swept-Wing Analysis previous two problems, is shown in Fig. 8. The solu- | 
tion for the 67 unknown quantities is listed below. tT 
The swept wing of the previous problem is enlarged Calculations were carried out on the IBM 701 machine 
by adding another spar and adjoining structure, as in about 12 minutes computing time, after some hand - 
shown in Fig. 7. All of the basic geometry and gages, computation of an elementary nature had adjusted the i 
and the basic assumptions are the same as those con- form of the matrix to yield a rapid machine solution. t 
+ 
P, = 0.9193166 ga = —0.04833938 gqya = —O0.0477282. w. = —54.982343 us = — 2.0296676 t 
P, = 0.7604357 qa = —0.108524381 gia = —0.1SS62S6 uu, = 2.3617836 vs = —1.1147309 + 
P; = 0.5251984 ga = —0.17302184 qua = — (0.032459 i = —1.0620915 ws = —12.274127 
Py = 0.2754757 qa = —0.16589635 qsa = —O0.0180125  w, = —23.217530 59 = 0. 86083738 E 
P, = 1.0012246 q;sa = —0.03173287 qiya = —0.0118582 uy = 1.0761495 uy = 0». 84043644 t 
Ps = O.S410177 =gea = —0.13046165 goa = —0.02162213 ov, = —0.63746158 vo = —1.5171785 +- 
P; = 0.5799460 gia = —0.25054666 qua = —(Q.04772S1 @y = —6.3389383 ww = —74.408520 f 
P, = 1.0794616 qsa = 0 .0323480 G2 = —0. 1409010 6, = 0.38S958S100 wy, = 3. 7009748 LC 
Py = 1.3985520) qga = 0 .0503569 G23 = —0.05713325 us = 9.8731345 on. = —1.4434195 : C 
Py = —0.0159951 qua = 0.06221476 uy, = 5. 9608525 vt = —1.5011842 wy = —26.162823 | | 
Py, = —0.0421719 gna = 0.08383721 2, = —1.4851891 ws = —89.738974 bp = 1.8143812 § £ 
Py» = —0.0526394 gra = 0.01538240 w, = —104.83257 u; = 3.8719099 : eo 
P;; = 0.02875763 qua = 0.1382701 m = t.0415359 2 = —1.4012477 { 
Pi, = —0.1155910 gua = 0.18354615 m = —1.3590758S w; = —41.090664 r 
a 
To check the accuracy of these results, one may sub- elements, (e) free bodies of whole cross sections, and 
stitute them back into the original matrix. When this (f) dummy load artifice. 
check is made, very small errors are found to exist. (3) Only one combined matric equation is required, | 1" 
Perhaps an even stronger check of the accuracy of the and the entries are but ove number for each element in | 1 
above values is obtained by writing equations of equi- the matrix—i.e., no sums of numbers are required. BE 
librium of whole cross sections of the structure. For (4) Removal, addition, or size change of structural a 
example, consider the equation of equilibrium of a sec- elements handled easily. a 
tion taken at the rigid rib: (5) Thermal loads handled easily. tT 
Po + Pis + Ps + Pu + Po — 2(qua + qa) = (6) Mixed boundary conditions handled easily. 7 
3.0000063-4 ee een ze 
Since the correct result should be 3, the error here is a is 
about 2 in 10°. Since only eight significant figures 
were carried in the calculations, this error is entirely | ¢ | 
satisfactory. The check of vertical equilibrium at this 
section is: a Dh 
a(gu + goo + gis + ges + gis) = —0.00000062 » in 
which is an error of 6.2 in 10°. 
9 
CONCLUSIONS - 
The matric method of structural analysis here pre- sui 
sented is believed to have the following advantages 
and disadvantages: 
nes 
Advantages 
(1) simplicity and speed of setup. oti 
(2) yields simultaneously all deflections (and, hence, ani 
influence coefficients) and forces without the use of (a) est 
quadratic energy expressions, (b) redundant designation an 





and “‘cutting,’”’ (c) locking of joints, (d) combining of Fic. 7. Three span swept wing. ing 
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Fic. 8. Matric formulation of three span swept-wing problem. 
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An Experimental Investigation of Leading- 


“dge Shock-Wave - Boundary-Laver 
Interaction at Mach 5.8' 


JAMES M. KENDALL, JR.* 
California Institute of Technology 


SUMMARY 


lhe boundary layer on a slender body tends to be very thick at 


interacts with the external flow by pro 


hypersonic speeds. It 
jucing larger flow deflections near the leading edge than those 
The increased shock strength affects the 
boundary-layer growth. The flow 
s rise to an induced pressure with a negative gradient which 


jue to the body alone 
iround the boundary laver 
give 
thins the boundary layer and increases the skin friction with re 
spect to the zero pressure gradient value 


Experiments on a flat plate with a sharp leading edge (Re 
100) have been performed in the GALCIT 5- by 5-in. Mach 5.8 
hypersonic wind tunnel. The induced pressure was measured 


Profiles of 


by means of orifices in Mach 


Number, velocity, mass flow, pressure, and momentum deficiency 


the plate surface 


were calculated from impact pressure surveys normal to the plate 
surface made at various distances from the leading edge 

The results are as The 
25 per cent higher than the weak interaction theory. (2) The 
are distinctly separate for 


follows. (1) induced pressures are 


boundary laver and the external flow 


Re, as low as 6,000. (3) The shock wave location is in good 
igreement with that predicted by the Friedrichs theory for a 
body shape equivalent to the observed boundary-layer displace 


ment thickness. (4) Expansion waves reflected from the shock 


ire weak. (5) The average skin-friction coefficient tends toward 
ind nearly matches the zero pressure gradient value downstream, 


but increases to approximately twice that value as the leading 


dge is approached 


SYMBOLS 


= speed of sound 


( = constant in viscosity relation w/us = C T7 
( = local skin-friction coefficient 
Cry = average skin-friction coefficient 
K = similarity parameter .1/6@ 
= mean free path 
17 = Mach Number 
p = static pressure 
= velocity 
Re = Reynolds Number per in 
Re, = Reynolds Number based on 4 
> = entropy 
= leading edge thickness 
1 = static temperature 


velocity in x-direction 
velocity in y-direction 
\ = distance from leading edge 
distance normal to plate surface 
= ratio of specific heats 
6 = boundary-layer thickness 


Presented at the Hypersonics Session, Twenty-Fourth Annual 
Meeting, IAS, New York, January 23-26, 1956 
[he work discussed in this paper was carried out under the 
sponsorship and with the financial support of the Office, Chief 
4 Ordnance, and the Office of Ordnance Research, U.S. Army 
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6* = boundary-layer displacement thickness 

6 = flow deflection angle 

6 = momentum deficiency (pit/pou ] u/u 
(1/yM pip | 

Mu = viscosity 

v = kinematic viscosity 

p == density 

x = interaction parameter J/*° YY C/V Re 

w = vorticity 7 u 


Subscripts 


= outside boundary-laver edge 
s = refers to shock-wave location 
= value at wall 
Q) = isentropic stagnation value 


= undisturbed stream 


INTRODUCTION 


HE KINETIC ENERGY of a hypersonic flow may be 
"T indefinitety larger than its thermal energy. Thus, 
when kinetic energy is converted into heat, small 
changes in speed such as occur in the inviscid flow field 
about a slender body are reflected as proportionally 
large changes in temperature, while at a stagnation point 
or in an uncooled boundary layer the temperature may 
rise to many times the free-stream value. These large 
temperature variations serve to distinguish hypersonic 
flows from supersonic ones with respect to both non- 
viscous and viscous effects. For example: (1) small 
velocity perturbations make such large changes in the 
local speed of sound that the linearized equations of 
motion of a supersonic flow are accurate only when the 
parameter .1/@ is everywhere much less than unity; (2 
when the Mach Number of a flow with fixed static tem 
perature becomes high enough, imperfect gas effects, 
including dissociation and ionization, will occur because 
the temperatures become large; (3) the low density of 
the fluid within the compressible boundary layer en 
hances the possibility of slip-flow; and (4) the thickness 
of the hypersonic laminar boundary layer tends to 
increase like .\/* for a given Reynolds Number, and 
may be so great that considerable interaction with the 
free stream results. This phenomenon is the subject 
of the present experimental investigation. 

So-called shock-wave - boundary-layer interaction oc- 
curs at supersonic speeds when the rapid growth of the 
boundary-layer displacement thickness near the lead- 
ing edge produces an effective thickening of the body. 
The resulting disturbance to the oncoming stream adds 
an induced pressure with a negative gradient to the 
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inviscid pressure. ‘The induced pressure has a tendency 
to thin the boundary layer and increase the local skin 
friction. Simultaneously, the shock wave emanating 
from the leading edge must increase in strength to 
account for the added flow deflection, so that stream- 
lines crossing it carry increased values of vorticity and 
entropy. In this regard, even for bodies which are 
slender in the usual sense, the leading-edge thickness is 
a parameter of great importance because of its strong 
influence on the shock strength. The streamlines cross- 
ing the shock travel very large distances downstream 
before entering the boundary layer because the relative 
mass flow in the boundary layer is very small for hyper- 
On blunt bodies the pressure induced by 
because the 


sonic flow. 
interaction is generally insignificant 
boundary layer is kept thin by the high pressures and 
gradients of the inviscid flow, but on sharp-nosed 
slender bodies the induced pressure may amount to 
several times the inviscid pressure. 

Several theoretical studies of this problem, princi- 
pally for a flat plate with an infinitely sharp leading 
edge, have been made by Bertram,' Kuo,’ Lees* # 
Lees and Probstein,® Li and Nagamatsu,® Maslen,’ 
and Shen.’ These theories have followed two general 
approaches. Shen® and Li and Nagamatsu® have 
applied the Karman momentum integral technique to 
the boundary-layer equations in the entire region be- 
tween the plate surface and the shock wave, not recog- 
nizing a boundary layer edge. The other solutions 
have considered the shock wave and boundary layet 
as separate, and have calculated the inviscid flow over 
a body represented by the boundary-layer displace- 
ment thickness using the concepts of hypersonic simi- 
larity and certain approximations. Much of this work 
is succinctly reviewed in references 9 and 10. Experi- 
ments have been performed by Baldwin,'' Bertram,'” 
and Hammitt and Bogdonoff.* '* These have been 
mainly concerned with the measurement of the induced 
pressure, although an optical determination of the 
shock wave location has been included in reference 15. 

It is found that, in general, the theory does not ac- 
curately predict the measured induced pressure dis- 


As an extreme example, the induced pres- 
91 


tribution. 
sures measured in the helium tunnel’ are about 
times those expected from the theory. Therefore, it 
appears that the theoretical fluid-mechanical model is 
either deficient in some respect, or is inconsistent with 
the physical model, especially with regard to the blunt 
leading edge. These discrepancies have prompted the 
present experiment. The aim of this experiment has 
been to measure and map the flow field for the case 
where the leading edge is so sharp as to be of vanishing 


importance. 

It will be of interest to interpret some of the experi- 
mental results in connection with those theories which 
treat the shock wave and boundary layer as distinct. 
Therefore, some of the concepts involved will be dis- 
cussed and certain significant parameters deduced 
by means of order of magnitude estimates. 

There exist two regions of flow behind the leading- 
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edge shock wave—-one near the body surface where 
viscous or boundary-layer effects predominate because 
high temperatures cause the viscosity to be large and 
the density to be small, and an external one where the 
density is large and the viscosity small so that dynamic 
forces predominate. The mass flow deficiency of the 
hot layer is balanced by an excess mass flow in the cool 
layer. These regions are distinguished very clearly 
for the laminar hypersonic boundary layer with no inter- 
action, but near the leading edge where interaction 
occurs, such a separation of predominating effects is 
subject to restrictions arising from the boundary-layer 
equations. In references 3 and 6 it has been shown 
that the Prandtl boundary-layer equations remain 
valid in the presence of interaction as long as (6/x 

|. This restriction essentially demands that the pres 
sure variation across the thickness of the boundary 
layer be small. 

Two characteristics of the flow outside the boundary 
layer are particularly important. (1) The flow is prac 
tically inviscid. A comparison of the kinematic vis 
cosity here with that in the boundary layer by means 
of the linear viscosity relation p/y. = C 77. gives 

»./b = [1/C(T/T.)?] ~ (1/CM,') < 1 

where the bars denote average values in the boundary 
layer. Also, the shear introduced by the curvature 
of the leading-edge shock is much smaller than in the 
boundary layer. (2) The expansion waves of the flow 
over the convex boundary-layer edge which are inct- 
dent upon the leading-edge shock are practically ab- 
sorbed there and, thus, serve to weaken it. For ex- 
ample, according to reference 14, the reflected strength 
of a wave is less than four per cent at Mach 10 for flow 
deflections under 40 degrees. It follows from these 
two characteristics that the flow in this region approxi 
mates a simple wave state, so that there exists a relation 
between the pressure and the flow deflection at every 
point. 

This pressure-angle approximation divides inter 
action into the categories of weak and strong. For 
the weak case, the values of A = 1/6 are small enough 
so that the pressure is proportional to 0, corresponding 
in two dimensions to the first term of the tangent- 


wedge formula for A <1, 
(p/po) — 1 = yK + [y(1 + y)/4] A? + 


(See reference 10.) The parameter upon which the 


induced pressure depends may then be deduced. The 


mass flow in the boundary layer is 
m = p.u,(6 — 6*) 
and dm/dx = p,u,((db/dx) — 6] 
These give 
6 = (dé*/dx) — [(6 — 6*)/p.u.] [d(p.u,.)/dx | 
On a flat plate, the boundary layer grows like 


6 ~ (ax/pus) = [(B/uo)/(B/P@)) > (UaoX/ Pollo) 
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LEADING-EDGE SHOCK-WAVE 


and p p. take into account the viscous 
heating. Using the linear viscosity law, and taking 


p p. tobe 1 i. 


where @ wu 


6~ (IVC V Re,) x 
If # is taken to be dé*/dx, the induced pressure becomes 


b.) — 1 = yM(db*/dx) ~ (APV CV Re,) = x 


In the strong interaction case, A is large and the 
pressure varies more nearly like 6°, as in the tangent 


wedge formula 


$]A? + 
yK | 


b.) 1 = [ly +1) 


+ [iy +1) 42K}? 


/ 


The same estimate of the boundary-layer growth may 


be made by taking 


3/p. = (T./T)- (p/p and b/po ™ M*(b/x 
Then 

b/po ~ IPVCWVR 
and 6 ~ (VAMC */ Re,'*)x 


The restriction for the Prandtl boundary-layer equa- 


tions to hold becomes 
x/MP <1 


The occurrence of slip-flow is most probable near 
the leading edge where the shear is highest, and in the 
boundary layer where the mean free path, /, is largest. 
It is assumed that the no-slip condition is satisfied if 


l(du/dy) <u 


This restriction is to be applied to the strong interaction 


region. In the free stream, yu ~ Pola@.. In the 


boundary layer / ~ 1/pso that 


lle = The/pT. ~ V Re:/MVC 


du dy is taken to be u.. 6, or may be obtained from the 


theoretical expression for C;. The restriction becomes 


6/x)/M = 1/V M Re,'“C*™* < 1 


For a fixed Reynolds Number, this is more easily satis- 
fied as the Mach Number increases because / decreases 
faster than the shear increases. Jn the weak interaction 


region, the restriction is 


(6/x)/M = M/VC WV Re, <1 


DESCRIPTION OF THE EXPERIMENT 


Model and Equipment 


The experiment was performed in the GALCIT 5- by 
5-in. Mach 5.8 hypersonic wind tunnel, which is of the 
All tests 
were made under steady-state conditions. Ahead of 
the nozzle the air stream was heated to 225°F. to pre- 
The supply pressure ranged 


closed-circuit continuously operating type. 


vent air condensation. 
between 25 and 95 psia, depending on the test-section 


Reynolds Number per inch desired. Extensive facili- 


BOUNDARY 


LAYER INTERACTION 1 


ties are provided for filtering and drying the air. More 
complete descriptions of the tunnel and its associated 
instrumentation are to be found in references 15 and 16 

The high carbon steel flat plate model shown in 
Fig. 1 was 7 in. long, 5 in. wide, and 3/8 in. thick. 
The bottom of the leading edge was beveled 15 
By comparison with a reference under a high-power 
microscope, the leading-edge thickness was estimated 
to be 0.0002 in., except at a few places where it was 
as large as 0.0005 in. For the measurement of the 
pressure induced by the growth of the boundary layer, 
there were 10 orifices in the top surface between x 
0.050 and 4.00 in. spaced approximately equidistantly 
in the scale 1/+/x. The plate was located in a region 
of the test section where the static pressure was con 
stant to within 2 per cent, and was aligned with the 
measured flow direction. The flow was regarded as 
two-dimensional. 

Three impact pressure probes with frontal heights 
0.0025 in. (seen in Fig. 1), 0.005 in., and 0.010 in. were 
The 


tips of the probes were made from stainless steel hypo 


used in surveys of the flow field about the plate. 


dermic tubing by flattening one end and honing to 
proper shape. These tips were soldered into off-set 
holders. 

The plate surface pressures were measured by means 
32-tube vacuum referenced silicone oil manometer 
The impact 


ofa’ 
The tubes were read directly to 0.1 mm. 
pressures were read on a vacuum referenced mercury 


micromanometer to 0.1 mm. 


Test Program 

The plate surface pressures were measured several 
times at each of three supply pressures. 

Two preliminary tests concerning impact pressure 
probes were made. First, because a flow deflection of 
up to 10° behind the leading-edge shock was expected, 
a test of probe sensitivity to angle of attack was made. 
As the angle of attack of the 0.0025-in. probe was in 
creased from zero, it read 0.1 per cent lower at 3°, 0.5 
per cent lower at 6°, and 1.2 per cent lower at 9°. Ac- 
cordingly, during surveys the probe was maintained 
at about 4° with respect to the free stream so that the 
angle of attack error should have been small at every 
point of measurement. Second, it is well known that 
when an impact probe being used in boundary-layer 
measurements is not small compared to the boundary- 
layer thickness, considerable disturbance to the bound- 
ary layer results. Fig. 2 shows the resulting change 
in the impact pressure profile normal to the plate sur- 
face as the size of the probe is increased. The largest 
of the three probes caused the boundary layer to 
thicken and to show a pressure bump at its outer edge 
characteristic of probe interference. This is tenta 
tively explained as follows: the static pressure in the 
vicinity of the probe is increased by the presence of 
the probe, separating the low dynamic pressure region 
of the boundary layer for some distance upstream. 
The resulting deflection of the flow is equivalent to a 
thickening of the boundary laver, and at the same time 
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LEADING-EDGE SHOCK-WAVE 


increases the efficiency of the shock system in front of 


the probe, so that it may read higher pressures 


\ detailed examination of the flow over the flat plate 


was undertaken by measuring the impact pressure pro- 
files normal to the plate surface at 10 stations between 
\ 0.1 and 3.00 in. for each of three supply pressures. 
The forward limit, x = 0.1 in., 
boundary layer being so thin that even the 0.0025-in. 


probe caused excessive interference. 


Data Reduction 


It was desired to present the induced pressure data 
as Pp px versus x in accordance with the interaction 
theory. was given by the pressure orifice readings, 
and p.. was taken as the empty-tunnel pressure at the 
point in space corresponding to each orifice location. 
The values of \/.., calculated from the impact pressure 
ahead of the leading edge shock, were 5.66, 5.79, and 
5.80, respectively, for the low, medium, and high Reyn- 
olds Number conditions. The viscosities used in 
calculating the Reynolds Numbers and 


C= gol o/tel< 


were determined from the Sutherland formula. For the 


present case, C = 0.90. 

The location of the boundary-layer edge, 6, and the 
shock wave, y,, were taken from the impact pressure 
profiles. The artificial definition of 6 shown in Fig. 3 is 
not subject to the usual ambiguity. , was taken to 
be that distance above the plate where the impact 
pressure reached its maximum value. 

Using the impact pressures and an independently 
determined estimate of either the total pressure—the 
pressure recovered by an isentropic compression to 
zero speed—or the static pressure at the same points 
in space, the following quantities were computed for 
the intermediate Reynolds Number condition: Mach 
Number, velocity, mass flow, momentum deficiency, 
and either static or total pressure, whichever was not 
estimated. Isoenergetic flow was assumed. 

Outside the boundary layer, the total pressure along 
a streamline was estimated by determining where that 
streamline crossed the shock wave by consideration 
of the mass flow. The measured strength of the shock 
along its length then determined the total pressure. 
The details are as follows. (1) The impact pressure 
ratio across a shock gives an accurate measure of its 
strength. From the impact pressure profiles a plot of 
total pressure ratio across the shock as a function of the 
shock location, y,, was made. (2) In Fig. 4, several 
streamlines are shown crossing a curved shock. The 
y,S represent the points at which the impact pressure 
was measured, and the y,,’s were sought. If pu is the 
average mass flow between yp and 4y;, then by continuity, 


Vs, = Veo — (Ptt/ Palko) (Yo — Vs) 


3) u was taken to be the speed of the flow with negli- 
gible error because the inclination of the flow was 
small. For a first approximation, py was taken equal 
to the value of pu at yo, which was known from the 
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was imposed by the 


oblique shock relations, and an approximation to y 
was obtained. (4) The total pressure ratio for this 
value of y,, and the impact pressure at y,; were used 
to compute pu at vy. (5) An improved estimate of 
pu Was made, and the process repeated until a con- 
sistent value of y,, was obtained. (6) By similar 
means, y,, Was computed using y,,, and so on until the 
streamline being computed crossed the leading-edge 
shock at a value of y, lower than for which there was 
information about the shock strength. (7) The other 
quantities listed above were computed with the aid of 
the compressible flow tables. 

Using the static pressures thus calculated and the 
measured plate surface pressure, the static pressure was 
estimated for the remainder of the flow field. The 
calculations were continued using the impact and 
Static pressures. 

6’, the “momentum” deficiency, has the interpre 
tation of being one-half the contribution of a stream 


tube to the total drag of a body. 6’ is defined by 


(pu/powt.) [1 — (wu.)|] — Ay \(p/p — |] 


For the present case, the average drag coefficient at any 


x-station is given by 
ey 
Cr = (2x | 6’ dy 
ZW 


Outside the boundary layer, 6’ is proportional to the 
mass flow times the entropy rise through the shock 
wave and, hence, for weak shocks, vanishes like (.1/,,7 — 
1)*, where JJ, is the component of J. normal to the 
shock. 

The mass flow profiles were integrated with respect 
to y by means of a mechanical integrator, and several 
streamlines were traced by locating points of equal 
value of the integral at each x-station. Several Mach 
lines were traced by a simple quadrature of the sum of 
the local Mach angle and streamline inclination. 


Discussion of the Accuracy 

The plate surface pressures were measured with less 
than 1 per cent error. The scatter of any one orifice 
reading was occasionally as high as 3 per cent, but the 
repeatability of the best straight line in a plot of p p 
versus 1/+7/xX was within 2 per cent. The viscosity 
used in calculating the Reynolds Numbers and C was 
known to within 5 per cent, and J... was known to 0.2 
per cent, so that x is accurate to within 6 per cent. 

In the boundary layer, the flow quantity profiles re- 
flect errors due to (1) impact probe location, (2) probe 
error—low Reynolds number effects, etc., (3) estimated 
static pressure error, and (4) assumption of isoener- 
getic flow. Near the plate surface, errors (2) and (4) 
are expected to produce large errors in each of the pro- 
files. These errors may, however, be circumvented to 
some extent by estimating values which are consistent 
with « u.. being zero at the plate surface and nearly 


linear in y. At the value of y where u/u. = 0.5, 
U/Uo, Pu; Ppota, M, and po/p 
affected by errors (2), (3), and (4), but 6’, which reaches 


may be somewhat 
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approximately its maximum there, is relatively un- 
affected. In the outer boundary layer and in the re- 
mainder of the region where the calculations were 
based on the estimated V/s, 
pu/p. u., and M are reasonably accurate, but 6’ is not 


static pressure, 
reliable because of its dependence on [1 — (u u.)]. 
Since this region extends over at most one-third of y,, 
and because 6’ is small compared to its maximum, the 
error contributed to the integral of 0’ is small. In the 
range of y where the calculations were based on the 
estimated total pressure, all quantities are believed to 
be nearly as accurate as the measured impact pressures. 

One check of the overall accuracy is given by how 
well the mass flow profiles satisfy the continuity rela- 


PVs 
(1 vo) | (pu/pottn) dy = 1 


0 


tion 


The value of the integral was 0.97 at x = 0.157 in., 


and was between 0.98 and 1.00 for all other stations. 
It is felt that the integrals of 6’ are probably accurate 
to within 5 per cent. 


RESULTS AND DISCUSSION 


Induced Pressures 


The induced pressure ratios are shown versus | +x 


in Fig. 5, and versus the interaction parameter x in 
Fig. 6. In these Figures, the leading edge lies to the 


right, the trailing edge to the left. Because the Mach 
Number was nearly constant, the present experiment 
offers no justification for the use of the interaction 
parameter other than permitting a direct comparison 
with theory. However, the correlation was improved 
by including the slight variations of Mach Number. 
Also shown in the latter figure are the first-order weak 
and strong interaction theory results, 

1 + 0.35x and = 0.92 + 0.52x 


p P« = Pp Pp 


The data are seen to differ from 
the weak interaction 


(See reference 10.) 
the expected results in two ways 
theory should be adequate for x < 1, but the data lie 
25 per cent higher; and the data are strikingly linear, 
failing to show any second-order variations or any 
transition between weak and strong forms. 

The discrepancy between the data and the weak inter- 
action theory might be due to one of the following 
causes. 

(1) The finite leading-edge size in the experiment. 
This possibility is to be dismissed for two reasons. 
First, in measuring the induced pressures on a flat 
plate, Bertram!” has found that the variation in the 
plate surface pressure as the leading-edge thickness, 
t, is changed is approximately equal to the correspond- 
ing variation of the theoretically calculated pressure 
of an inviscid flow over a blunt wedge-plate body. 
In the present experiment, the first pressure orifice was 
250 leading-edge thicknesses downstream, and accord- 
ing to Bertram’s calculation the leading edge makes a 
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negligible contribution to the measured 


Second, no systematic change in the slope of the curves 


pressures. 


of p/p. versus x is observed although Re, has varied 
by a factor of four. 

(2) An inaccuracy of the pressure-angle approxi- 
mation used in the theory. This will be checked by 
comparing the tangent-wedge pressure with the meas- 
ured pressure along a streamline. 

(3) An inaccuracy in the theoretical boundary-layer 
growth. This will be checked by comparison with the 
experimentally determined boundary-layer thickness. 

(4) Error resulting from the neglect of the second 
term of the expression 

6 = (d6*/dx) — {(6 — 6*) ‘p,u,] [d(p.u.) dx] 


Taking 


p/po = + 0.45x 


(1/p.u.) [d(p.u,)/dx] = (1/yp,) (dp, ‘dx 


5* = 0.756 ~ vx 


the ratio of the second term to the first becomes 
r = (0.11x)/(1 + 0.45x) 


The flow deflection is underestimated by about 7 per 
cent for x = 1 and by 14 per cent for x = 3. 

(5) A pressure variation across the boundary-layer 
thickness. The flow quantity profiles support this 
possibility. 


Boundary-Layer Growth 


The impact pressure surveys were found to define a 
boundary-layer edge very distinctly for even the low- 
est Re,, about 6,000, where the probe frontal height 
was 15 per cent of 6. These data are presented in 
Fig. 7. The straight lines faired through them are 
proportional to x°-*, x°-°, and x°-*!, respectively, for 
the high, medium, and low Reynolds Number per inch 
conditions. In Fig. 8, 
6/x versus M2/C/V Re, in accordance with theory 
The satisfactory correlation is not surprising in view 
The 
thinning of the boundary layer near the leading edge 


these data are presented as 


of the good correlation of the induced pressures. 


by the induced pressure is clearly seen, and suggests 
that 6/x does not exceed some maximum value, in con- 
trast to the infinite value predicted by a simple power 
law. Thus the requirement (6/x)? < 1 for the Prandtl 
boundary-layer equations to hold may be satisfied. 

The boundary-layer displacement thickness, 6*, as 
defined by 


| pu dy = p,u.(d — 6*) 
/70 
may be written 
5* 1 put co pu, 
—_ ] — av 
6 oe te £0 Baa 


This quantity has been evaluated from the mass flow 
profiles (which are to be discussed) and is approxi- 
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LEADING-EDGE SHOCK-WAVE 


mately 0.75. According to Fig. 8, the displacement 
thickness for large Reynolds Numbers may then be 


written 
6* x = 0.50) IPA C/V Re,) 


The value used in the first-order weak interaction theory 


1S 
6*/x = 0.48(AM?y c V Re, ) 


The method of Cohen and Reshotko” for calculating 
the laminar boundary layer with arbitrary pressure 
gradient was used to predict a boundary-layer thick- 
ness for the case Re/in. = 117,000. Taking the ob- 
served value of 6 at x = 0.100 in. as the starting point 
of the calculation in order to avoid making assump- 
tions about the flow conditions very near the leading 
edge, and using the measured plate surface pressures, 
the resulting boundary-layer growth was found to 
agree with the experimental values to within 5 per cent 


up tox = 31n. 


Shock-Wave Location and Strength 


The shock location, y,, is given for three Reynolds 
Numbers in Fig. 7. When these data are presented 
as Re, versus Re,, as in Fig. 9, it is found that they may 
be empirically represented as Ke, ~ RKe,'. It is 
interesting to note that in the Princeton helium tunnel 
at Mach Numbers of the order of 14, the shock shape 
for a flat plate with Ke, = 119 may be expressed as 
Re, ~ Re,’ (see reference 13). 

However, it is perhaps more appropriate to present 
these data in another way. The shock-wave location 
on a body equivalent to the observed parabolic bound- 


ary-layer shape 
6* = 0.50( MPV C/V Re) Vx 


has been calculated according to the Friedrichs theory, 
as outlined in 397-402. In this 
method, the flow behind the shock is considered to be 


reference 1S, pp. 


a simple wave, and the equation of a Mach line may 


be written 
y — x tan [(1 8) + 0] = YO 


8 = VAM." — 1. Along the shock, dy dx 
obtained from this expression must be set equal to 


tan!n(@)], giving the linear differential equation 


where 


}tan [(1 8) + 6] — tan [n(6)]} (dx dé) + 
x(d/d@) tan [(1/8) + 6] = —Y'(@) 


To the first order one has 


tan [(1/8) + 6] = (1/8) + 2H6 
tan [n(@)] = (1/8) + 16 
where 
A= [(y + 1)/4] (W.4/8') 
Integration gives 
] bd] : 
= , 6 dY(@) 
A@ Jeo 
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where 6) = 6 at x = 0. Together with the equation of 
the Mach lines, this expression gives a parametric rep- 
resentation of the shock shape. Taking 6 = dé6* dx, 
one obtains for the present case 


Rely, — (x’8)] = 3.13Re?3 


This result is shown in Fig. 10, and the experimental 
data are seen to compare favorably with the theoretical 
prediction. 

In Fig. 11, the flow deflection behind the shock wave 
This 


measure of the shock strength has been determined 


is given for the intermediate value of Re in. 


from the impact pressure ratio across the shock and is 
considered to be quite accurate. The same quantity 
has been calculated from the empirical relation y, ~ 
x°-s, which represents the best fit of the shock location 
for this value of Re/in. The agreement is seen to be 
rather poor. 

As a check, the reverse procedure was followed. The 
shock-wave angle determined from the impact pressure 
ratio was integrated by a simple quadrature, starting 
from the observed shock location at x = 0.1 in. The 
shock location determined by this means was virtually 
identical with that measured. It is to be concluded 
that the determination of the strength of a weak shock 
by differentiation of its measured location is a matter 


requiring some care. 


Flow Quantity Profiles 


Certain interesting features of the flow field are dis- 
closed by the profiles of the calculated flow quantities 
at the various x-stations. One such set of profiles is 
shown in Fig. 12. For example, the velocity in the 
viscous region appears to retain the form usual for the 
laminar compressible boundary layer, while in the ex- 
ternal region the velocity is near the free-stream value, 
and the mass flow ratio is larger than unity. The ex- 
tent to which there exist separate viscous and non- 
viscous regions does not diminish in the profiles nearer 
the leading edge. In contrast to the usual viscous flow 
over a flat plate, the momentum deficiency does not 
vanish outside the boundary-layer edge, but extends 
to the shock because of the wave drag of the boundary 
layer. 

A substantial pressure gradient in the y-direction due 
to the streamline curvature occurs in the inviscid flow 
region. The mass flow, Mach Number, velocity, and 
momentum deficiency show variations which are con 
sistent. Above a certain value of y the velocity and 
Mach Number decrease with increasing y because of the 
increasing static pressure and slowly varying total 
pressure, while the momentum deficiency increases due 
to the increasing mass flow. Below this value of y these 
quantities show the opposite variation because the total 
pressure change with y is rapid. 

A measure of the vorticity outside the boundary-layer 
edge is afforded by the total pressure profiles. Crocco’s 
theorem is 
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where g is the velocity, w the vorticity 7, — u,, 7 the 


temperature, and VS the entropy gradient. The entropy 


may be written 


S — S. = C, In (po,/fo)"~' 
and its gradient is approximately 
VS = dS/dy = 
i . _ | 7 C (d dy) |( Po — Dp p a 


Rid dy) (Ap Po) 


rhe vorticity at the plate surface is nearly w, = u, ‘6, 
so that 


w / wy = (6RT/u,”) (d/dy) (Apo/po) = 
(6 7 M, : (d dy) (Ap Pr ) 


At x = 0.25 in., for example w/w, = 0.1 outside the 
boundary-layer edge, but is only 0.005 near the shock. 
The rate of shearing strain v, + uw, near the shock will 
be approximately 2u, < 0 since w = 0. 

Some anomalies remain to be considered. (1) The 
behaviors of each of the flow quantity profiles near y = 
0 reflect the errors suffered by an impact probe at low 
Reynolds Numbers. (2) The profiles at x = 0.157 
and 0.250 in. (not shown) show the effects of impact 
probe interference near y = 6. (3) According to the 
profiles at x = 0.64 in. and aft, the pressure is about 
) per cent higher in the boundary layer than in the 
nearby external region. Although such an effect has 
not been explained, no reason to discount this result 
has been found. 


Flow Field 


Fig. 13 gives an experimentally determined repre- 
sentation of the complete flow field. It is seen that 
the streamlines cross the shock quite far ahead of where 
they enter the boundary layer, and the Mach lines 
intersect the shock at very small angles, so that the 
rather slow decay of the shock strength, corresponding 
to its x’- shape, is accounted for. The boundary- 
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layer growth in the first 0.1 in. determines the shock 
strength up to about x = 0.7 in. 

In Fig. 14 the static pressure along each of the three 
most forward Mach lines of Fig. 13 shows only a slight 
increase with x. Apparently the expansion waves 
incident upon the shock are reflected as very weak com 


pression waves. The flow deflection along the first 


Mach line is calculated to be 4.56° at 0.157 1n 
and 4.25° at x = 0.64 1n. This Mach line is nearly 
straight, having an inclination of 15.40° at a 0.157 


in. and 15.21° at x 0.64 in. The other Mach lines 
show smaller variations along their lengths 

The static pressure along the streamline shown cross 
ing the shock at 4 0.1 in. in Fig. 13 has been calcu- 
lated from the local streamline inclination by means 
of the tangent-wedge approximation and compared 
with the pressures taken from the calculated profiles 


The results are shown in Fig. 15. 


Total Drag 


Near the leading edge the shear is expected to be 
higher than the zero pressure gradient value because 
the boundary layer is thinned by, and Ov Oy at y 0 
increased by, the negative induced pressure gradient 
The velocity profiles are not considered to be suffi 
ciently accurate near y = O to permit a direct calcula 
tion of the local skin friction. Instead, the average 
skin-friction coefficient, Cy, has been obtained from the 


6’ profiles and is shown in Fig. 16. The drag of the 
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leading edge, approximated as the impact pressure 
times half the leading-edge thickness, is about 7 per 
cent of the measured Cr at x = 0.157 in., and is less 
at succeeding stations. Eimer’ has measured the local 
skin friction on a flat plate in the same tunnel, and the 
average skin friction corresponding to his result is 
shown for comparison. His measurements were made 
by means of a floating element balance located suffi- 
ciently far from the leading edge that interaction 
effects may be neglected. 

Maslen’ has given a theoretical expression for the 


local skin friction of the form 
C, V Re, = f(M) + [g(M)/V Re] 


It is not directly integrable to give Cr. Instead, the 


experimental value of Cy at x = 3 in. was taken, and 
its variation ahead of this station computed from 


Maslen’s result according to the expression 


Cy, = (1/x) | C,dx = (3/x)Cr, — (1/x) | C, dx 


0 J x 


The values of Cy thus determined agreed with the ex- 
perimental results to within 2 per cent. 


SUMMARY OF RESULTS 
The following results have been found. The induced 
pressures are linear in 1°VC V Re, and are about 
25 per cent higher than the first-order weak interaction 
theory. 

The boundary-layer edge is well defined for Ke, at 
least as small as 6,000. The flow outside the edge is 
inviscid, but vorticity is present. 

Away from the leading edge, the boundary-layer 
thickness is in agreement with the zero pressure gra- 
dient theory, but is thinner near the leading edge. 

The shock-wave location may be represented by 
Re[y, — (x 8)| ~ Re,*’* and is in good agreement with 
the theoretical shock location for a parabolic body. 

The paths of the Mach lines indicate that the shock 
location ahead of x = 0.7 in. is determined by the 
boundary-layer growth ahead of x = 0.1 in. The 
Mach lines are so nearly parallel to the shock that the 
strength of the shock along its length decays slowly. 

The static pressure and flow deflection are nearly 
constant along the length of a Mach line, indicating 
that waves reflected from the shock are weak. 

The average skin-friction coefficient tends toward 
and nearly matches the zero pressure gradient value 
downstream, but increases to approximately twice 
that value as the leading edge is approached. 


CONCLUSIONS 


An experiment intended to indicate the behavior of 
the interaction of the viscid and inviscid flow regions 
near the leading edge of a flat plate has been performed 
without consideration of any theory for the case where 
the leading-edge thickness is negligibly small and at a 
Mach Number. The case 


moderately hypersonic 
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where the leading-edge thickness is an important pa- 
rameter has been discussed in reference 13. 

The experiment yields a picture of the flow which is 
consistent with the weak interaction theory of Lees, 
with the exception that the induced 
The 
described by simple wave theory since entropy gradients 


pressures are 


somewhat higher. external flow field is well 
are small and hence waves reflected from the shock are 
weak. The boundary-layer growth is just what one 


expects from the observed pressure distribution. 
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An Investigation of Effects of Certain Tvpes 
of Structural Nonlinearities on Wing and 
Control Surtace Flutter 


RUNYAN,* anp ROBERT E. ANDREWS* 


DONALD S. WOOLSTON,* HARRY L. 
Langley Aeronautical Laboratory, 


SUMMARY 


his paper presents some effects of nonlinear structural terms 
on the flutter of a wing capable of bending and twisting and, also, 
of a system including a control surface. Several type of non- 
linearities in the stiffness are investigated, including free play, a 
hysteresis loop, and cubic variations. These are introduced in the 
torsional degree of freedom of the wing and in the aileron stiffness 
in the three degree of freedom system. Calculations have been 
made on an analog computer. A wind-tunnel investigation of one 
case having free play in the torsional degree of freedom has been 
made, and good correlation between theory and experiment is 
shown. 

Ihe results indicate that the stability of a nonlinear system 
is highly dependent on the magnitude of initial displacements of 
the system from equilibrium. It is shown that in many cases the 
flutter speed is decreased by increasing the initial disturbance. 
The results also indicate that when a nonlinear system becomes 


unstable its flutter may become self-limited. 


SYMBOLS 


= coordinate of pitch axis measured from midchord, 
referred to semichord b, positive for pitch axis 
behind midchord 

= wing semichord, ft 

= coordinate of aileron hinge measured from mid- 
chord, referred to semichord 6 


h = displacement of system in vertical translation 

a = angular displacement of wing in torsion 

8 = angular rotation of control surface about hinge 

hy), ay, Bo = initial displacements in degrees of freedom h, a, 
and 8, respectively 

H = height of hysteresis box 

i = moment of inertia of wing or wing-control surface 
about a per unit span length 

Ig = moment of inertia of control surface about c¢ per 
unit span length 

cio. me spring constants in equation of cubic spring 

m = mass of wing or wing-aileron system per unit span 
length 

V = moment 


radius of gyration of wing or wing-control surface 
system referred to a, V [q/mb? 
r = reduced radius of gyration of aileron referred to 


C, av ‘Ta mb? 


R/ Ruin = ratio of amplitude of response of system with 
nonlinearity present to amplitude of response 
of linear system 

as = static moment of wing or wing-control surface 
system per unit span length referred to a 

S38 = static moment of control surface per unit span 
length referred to ¢ 

J = velocity of main stream 
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Vin flutter velocity of linear svstem 

V = flutter velocity of linear system for control surface 
frequency wg = 0 

7 = location of center of gravity of wing or wing 
control surface system referred to a, Sq / mil 

XZ = reduced location of center of gravity of control 
surface referred to c, yg mod 

6 = angular free play in tersion or control surface 
rotation 

h = wing-mass air-density ratio m/xpb* 

p = mass density of air 

w = circular frequency of driving force (see Fig. 2 

Wh = natural circular frequency of wing in vertical 
translation 

Wa = natural circular frequency of torsional vibrations 
about a 

wg = natural circular frequency of torsional vibrations 


of control surface about ¢ 


INTRODUCTION 


_- ALL PROBLEMS in mechanics contain some 
elements of nonlinearity. For the most part, how 
ever, the great fund of knowledge now available to the 
physical scientist or engineer has been made on the 
the 
In many cases, the approximations intro 


grounds of linear approximations of governing 
equations. 
duced by linearization are quite valid, and the results 
are sufficiently accurate for their purposes. On. the 
other hand, the introduction of nonlinearities into a 
problem may not only change the linearized solution 
but may even introduce certain new types of phe 
nomena which are not even indicated by the linear 
theory. 

In the aeronautical field, the demand for faster and 
still faster airplanes has resulted in the absolute neces 
sity for the more accurate prediction of the various air 
plane characteristics than has been required in the past 
It has become necessary to examine very closely the 
new problems which have arisen in order to ensure that 
no important nonlinear effects are being neglected. Also 
it becomes necessary to re-examine and recast some of 
the more classical problems in the light of the possible 
effects of nonlinearities which have not been studied in 
the past. 

The purpose of this paper is to examine some effects 
of the introduction of several types of nonlinear terms 
into the classical two and three degree of freedom flutter 
problem. 

In considering the flutter problem, three types of 
involved—namely, inertia, damping, and 


forces are 
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elastic. Each of these forces may arise from either the 
aerodynamics of the problem or from the structure. 
Nonlinearities may arise from the aerodynamics, for 
instance, from the effect of airfoil thickness, or from 
flow separation or breakdown. The structural side of 
the problem may also give rise to nonlinearities in the 
inertial, damping, and elastic forces, some examples of 
which will be indicated later. This paper is concerned 
essentially with structural nonlinearities and constt- 
tutes an extension of the work done in reference 1, in 
that a third degree of freedom, a control surface, has 
been added to the bending-torsion degrees which were 
studied in reference 1. 

In examining the literature on nonlinear vibrations it 
was soon evident that only the single degree of freedom 
cases had received extensive treatment and that the 
state of the art had not progressed to the point of 
handling analytically more than one degree of freedom. 
It was therefore decided to study the problem on an 
analog computer. The type of analog used in the 
analysis involved the principle of the electronic differen- 
tial analyzer, wherein successive integrations of voltages 
are performed by the use of feedback amplifiers. 

In order to gain an understanding of some effects of 
nonlinearities, the physical systems to be dealt with 
have been intentionally limited to simple ones. The 
flutter characteristics of two configurations will be ex- 
amined; one, a wing, or all-movable control, free to 
translate and pitch about some spanwise axis, and the 
other, a wing-control surface system with freedoms in 
translation and pitch of the wing and rotation of the 
control about its hinge. Both systems are treated as 
two-dimensional, and, in the flutter calculations, in- 
compressible indicial air forces are used. 

With regard to the nonlinear springs to be considered, 
little is known of the exact nature and magnitude of 
nonlinearities which exist in actual aircraft. In the 
study to be discussed, the nonlinearities, too, have been 
kept simple in form and were selected as representative 
of types which might occur. 
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TYPES OF NONLINEAR SPRINGS INVESTIGATED 


The types of nonlinear springs considered are il- 
lustrated in Fig. | and are represented by plots of the 
force or moment required to produce a given displace- 
ment. Shown at the left of the Figure is a flat spot type 
of nonlinearity which represents free play in the hinge 
or linkage of a control system. The effects of free play 
will be considered for two different conditions: in one, 
the equilibrium condition will be taken as in the center 
of the free play, as shown by the solid lines in the 
sketch; in the other condition, a preload will be as- 
sumed to act on the system so that the equilibrium 
position lies at some point on the linear arm, as shown 
by the dashed lines. 

The second nonlinearity to be dealt with is a type of 
hysteresis, shown at the right of the Figure. As force 
or moment is increased, displacement varies linearly 
until a point is reached at which a jump occurs, after 
which the system is again linear. On the return path, 
a corresponding jump occurs at another value of the 
force. A nonlinearity of this type may occur in the 
case of a control surface with free play if friction exists 
at some point in the linkage. In the case of a wing, this 
type of nonlinearity might represent the effect of rivet 
slip. 

The third type of nonlinearity treated has been 
termed a cubic spring, where the force exerted depends 
on the usual linear power of the displacement and, in ad- 
dition, on a term containing the cube of the displace- 
This can be considered either as a hard spring, 
which stiffer 
shown by the solid line, or as a soft spring, which be- 


ment. 
becomes as displacement increases, as 
comes weaker as displacement increases, as shown by 


the dashed line. In the case of a control surface, a 


hard-spring effect might be associated with power 


controls. In the case of a structure, a hardening effect is 
found when a thin wing, or perhaps a propeller, is sub- 
jected to increasing amplitudes of torsion. A soft-spring 
effect may be associated with panel buckling. 


EFFECTS OF NONLINEAR SPRINGS 
ON AIRCRAFT FLUTTER 


The three types of nonlinear springs illustrated in Fig. 
1 have been considered in flutter investigations con- 
ducted on the analog computer. For the case of free 
play, experimental wind-tunnel results have also been 
obtained for a two degree of freedom system. The free 


play nonlinearity is considered first. 


Effects of Free Play 

Before actually presenting flutter results showing 
effects of free play, it is instructive and of interest to 
examine the zero air-speed vibration characteristics of a 
simple single degree of freedom system with this type 
of nonlinearity. 

Vibration Characteristics of a System with Free Play— 
A usual preliminary step in the flutter investigation of 
an aircraft is a ground vibration survey to determine 


natural frequencies and modes of oscillation. Such a 
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survey is generally conducted by measuring amplitudes 
of the response of a structure to a sinusoidal forcing func- 
tion over a wide range of frequencies. The results of 
such a vibration survey of a single degree of freedom 
system with free play are shown in Fig. 2. In obtaining 
these results a small amount of structural damping was 
included. 

Results are shown in Fig. 2(a) in the form of an ampli- 
tude ratio RR 
free play to the response of the linear system 


of the response of the system with 
that 1s, 
the same syste:n with zero free play—against the fre- 
, of the frequency of the driving force 
The re- 


quency ratio w w 
to the natural frequency of the linear system. 
sponse curve for the linear system with damping present 
is shown by the dashed line. The solid curves show re- 
sults obtained with a constant amount of free play but 
with successively smaller amplitudes of the driving 
force. For the curve with maximum response ampli- 
tude, the driving force was the same as that used with 
the linear system. To obtain the remaining curves, the 
driving force was successively reduced to 1 2, 1/4, 
1 Sof this value. With each reduction in driving force, 
the amplitude of the oscillation relative to the region of 


and 


free play is reduced and the response peak occurs at a 
lower value of the frequency. In addition, for the lower 
driving forces, the response curve is flattened. Since 
the shape of the response curve is frequently used to 
measure structural damping, it is seen that the value 
obtained in this example would be highly dependent on 
the amount of driving force. Thus, both natural fre- 
quency and damping would be difficult to determine. 

In Fig. 2(b) these results are presented in another 
form, representing the effect of a variation in free play 
with the amplitude of the driving force held constant. 
Results are shown for the linear case, for an amount of 
free play 6, and for successively greater amounts of 


f 2 With each increase in the 


free play of 26, 46, and 86. 
amount of free play in the system the effective spring 
constant is reduced and the amplitude of the response 
is increased. This indicates that in an actual aircraft, if 
the amount of free play increases with increasing flight 


time, the effective spring constant will be reduced. It is 
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free play in the spring restoring force—driving force constant, 


free play variable 


of interest to note in this connection that the response 


amplitude and effective spring constant for a given 


amount of free play and a given driving force can be 
approximately determined by consideration of the work 
performed per cycle and comparison with the work per 
formed per cycle in an equivalent linear system. It will 
be shown, however, that the effects of free play on the 
flutter characteristics of a system will not necessarily be 
predicted by considering the equivalent linear system 
The effects of free play on flutter are considered next 
for a two degree of freedom system. 

Flutter of a Two Degree of Freedom System with Free 
Play 


sults and analog results, which were given in reference 1, 


For this case both experimental wind-tunnel re 


will be compared again in order to establish the ac 
curacy of the analog results. In Fig. 5 the configuration 
used in the wind-tunnel tests ts shown. 

The sketch represents a two-dimensional model which 
is free to oscillate in pitch and translation. The pitch 
ing degree of freedom is provided by a bearing-sup 
ported shaft. This system, in turn, is suspended be 
tween a pair of leaf springs on either side of the test 
section, so that the entire mechanism is free to translate 

In translation the system is linear. The nonlinearity 
is introduced in the torsional degree of freedom as 


shown in the detail sketch. A leaf spring is clamped to 
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the end of the torsion axis and its free end extends up- 
ward between two setscrews. 
screws can be closed completely to give a linear torsion 
spring, or opened to provide any desired amount of 
angular free play, giving the spring characteristic shown 
at the right of the Figure. Physical properties of the 
system are given in Table |. 

Flutter results obtained both in the wind tunnel and 
by analog calculations are given in Fig. 4. In this ex- 
ample a preload acted on the system so that at equi- 
librium the model was not centered in the flat spot, but 
rather at some point on the linear curve of the spring 
diagram. 

Flutter results are shown on the lower part of the 
Figure. In obtaining these results, the system was 
given an initial torsional displacement ay and suddenly 
released. In the Figure a is plotted as ordinate 
against a flutter speed ratio |’ |’);,, where I’,;, is the 
flutter speed of the system without free play. At each 
velocity the value of a» was increased until the model 
became unstable. Throughout the test the amount of 
free play 6 was kept constant at 0.5.° 

Analog results are shown by the solid curves, the 
circles show experimental points. It is seen that flutter 
occurred well below the flutter speed of the linear sys- 
tem and that three distinct regions were found—a stable 
region, a region of mild flutter, and a region of violent 


flutter. 
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TABLE 1 
PROPERTIES OF SYSTEMS 


Three Degrees 


Parameter Two Degrees of Freedom of Freedom 


Free play Cubic spring 


a —0.406 -(). 54 0.2 
Le 0.042 0.052 0 

Yor? 0.1603 0.24 0.5 
wr, rad. /sec 58.94 94.25 Bi .7 
Wa, rad./sec 81.24 125.66 136.0 
wp, rad. /sec variable 
( 0.6 
xB 0 

rg- 0.002 
1/« 156.405 79.724 12.0 
b, ft 0.5 0.5 2.375 


Because of the presence of the preload, the system 
oscillated on the linear arm of the spring for small 
initial displacements and, at velocities below V;;,, was 
stable. For slightly larger displacements the oscilla- 
tions entered the region of free play so that the effective 
torsional stiffness was reduced. In this example a de- 
crease in torsional stiffness was destabilizing and a 
limited amplitude flutter occurred. (This limited ampli- 
tude flutter was found both experimentally and on the 
analog. Experimental points for this mild flutter are 
not shown. ) 

As initial displacements were increased still further 
the oscillations covered more and more of the region of 
free play. <A further destabilizing reduction in stiffness 
occurred, and the flutter oscillations became violent. 
Both experimental and analog results for the violent 
flutter are shown. 

The results indicate good agreement between experi- 
ment and the analog. They also show that, with pre- 
load acting on a system with free play, for small initial 
displacements the system will act as a linear system. 

In the next example a three degree of freedom system 
(wing with control surface) will be treated. First, the 
effects of free play without preload will be considered. 

Wing-Control Surface Flutter with Free Play; Zero 
Preload—Results for a wing-control surface with free 
play and zero preload are given in Fig. 5. At the upper 
right of the Figure flutter characteristics of the cor- 
This 


linear case is one of those treated by Theodorsen and 


responding system without free play are shown. 


Garrick in reference 2. Properties of the system are 
included in Table 1. 

Results for the linear system are shown in the form of 
a flutter speed ratio |’ 1) against the ratio wg wa, of 
control surface frequency to wing torsional frequency. 
In the ordinate !y is the flutter speed for control sur- 
face frequency ws = 0. It is seen that, as the ratio 
wg Wa increases from zero to about 0.8, the flutter speed 
rapidly decreases. With further increases in the fre 
quency ratio the flutter speed increases sharply. 

To examine the effects of free play in the control sur 
face, two cases were studied: in one, Case A, the fre- 
quency ratio wg/wa in the absence of free play was 0.6; 
in the other, Case B, wg/ wa for zero free play was 0.95. 
It is noted in the linear system that, in Case A, the flut 
ter speed for ws = 0 is above that for w3 wa = 0.6, while, 
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EFPFPECTS OF CERTAIN TYPES 


jor Case B, the flutter speed for w,; = 0 is about half of 


F the speed for wg wa = 0.95. 


Flutter results for Cases A and B with free play are 
shown at the bottom of Fig. 5. Results are presented 
in the form of a ratio 8» 6, of initial angular displace- 
ment of the control surface to the angular free play, 
against the flutter speed ratio |’ 1’,;,. In the two sets 
Vii, used is the linear flutter 
The 


of results the value of 
speed appropriate to the case being treated. 
ymount of free play 6 was kept constant. 

In the results for Case A, the effects of free play seem 
to be beneficial in that, for small displacements, flutter is 
first encountered well above the flutter speed of the 
corresponding linear system. When the initial displace- 
ments are increased, the flutter boundary shifts to lower 
Vin. For Case B, the presence of 
For small values of 89/6 flutter 


values of the ratio 
free play is detrimental. 
is encountered at about | 2 of the flutter speed of the 
corresponding linear the flutter 
boundary moves to lower speeds with increasing initial 


system. Again, 
displacements. 

Actually, results for the two cases are quite similar. 
In the absence of any preload the control surface is 
initially floating in the region of free play, and, for small 
displacements, the oscillations are confined within this 
region where the effective value of wg is zero. For both 
cases, therefore, flutter occurs first at the same velocity 
|), the flutter speed of the linear system at we = 0. It is 
recalled that this velocity is above !”,;, for Case A and 
below I’,,, for Case B. 

As the initial displacements and the resulting flutter 
amplitudes are increased, the oscillations extend beyond 
the region of free play a proportionate amount, and the 
effective value of wg is increased. In the linear system 
an increase in wz; means a reduction in flutter speed, so 
that in both Cases A and B the flutter boundary shifts to 
lower speeds. If sufficiently high initial displacements 
could be attained, the flutter boundaries in the two cases 
would become asymptotic to the respective linear flutter 
speeds. 

Wing-Control Surface Flutter with Free Play; Preload 
Included—-Results for Case A of the system just con- 
sidered are given in Fig. 6, with a preload added to the 
control surface spring. For ease of reference, the curve 
showing variation of flutter speed with wg wa for the 
linear system is included at the upper left. It is re- 
called that Case A corresponds toa value of wg wa = 0.6 
ina system without free play. 

Results for this case with fixed amounts of free play, 6, 
and preload are shown at the lower left of the Figure. 
Results are again given in the form of 8) 6, the ratio of 
initial displacement to free play, against the flutter 
speed ratio VV jin. 

With preload on the spring the system is stable only 
it velocities below the flutter speed of the corresponding 
linear system. For values of |” J’,;, between 1.0 and 
1.3 a region of limited amplitude flutter occurs. At 
higher values of |’ ’,;, the flutter is divergent. This is 
iresult quite different from that found for Case A in Fig. 
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», Where no preload was considered and where the sys- 
tem was flutter-free at velocities well above |” 

When preload is applied the control surface oscillates 
about a point on the linear arm of the spring diagram. 
For small amplitude oscillations the system is linear and, 
therefore, flutters at = J’,,,. If, at the linear flutter 
speed, the flutter oscillation grows into the region of 
This 
means that the flutter speed is raised, as shown by the 
flutter curve for the linear system. Without an in 
crease in velocity, the flutter amplitude is therefore 
,» up to 


free play, the effective value of ws is reduced. 


limited. Through the range of values of |) V 
1.3 the flutter amplitude increases and takes in more 
and more of the free play. The effective w, is therefore 
continually dropping and the flutter speed increasing so 
that the flutter amplitude remains limited. Finally, at 
values of | V),,, greater than 1.3, the control surface 
oscillates completely through the free play and onto the 
upper linear arm of the spring diagram. The effective 
value of wg is immediately increased, the flutter speed is 
decreased, and the flutter oscillation becomes divergent. 

It should be pointed out that, if Case B had been con- 
sidered, divergent flutter would probably have been 
encountered immediately at the value of |’,,, for Case B. 
This system lies on the side of the flutter valley for 
which any reduction in wg is accompanied by a sharp 
reduction in flutter speed. 

It is interesting to note that, in performing flutter 
studies on the analog computer, the flutter amplitudes 
in each of the degrees of freedom can be continuously 
recorded. A plot of the variation with velocity of a 6 
(the ratio of torsional displacement to the amount 
of free play) and 6 6 (the ratio of control surface 
rotation to angular the 
limited amplitude flutter is shown at the lower right of 
the Figure. The value of 8 6 at |’ I’,;, = 1.8 is that 
for which the control surface oscillates completely 


free play) in region of 


through the free play. 

In the next examples to be considered, some effects of 
hysteresis are examined. 
Effects of Hysteresis 

If, in addition to the free play just discussed, some 
source of static friction exists in the system, the hys- 
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Response of a system with hysteresis. 


Fic. 7. 


teresis type of nonlinearity described earlier is en- 
countered. Before considering the flutter results, the 
response of a single degree of freedom system with 
hysteresis is examined, as was done in the case of free 
play. 

Response of a System with Hysteresis 
sents the response of a system with hysteresis after re- 


Fig. 7 repre- 


lease from an initial displacement at zero air speed. The 
upper trace shows the displacement of the system; the 
lower trace shows the corresponding variation in the 
moment. The indentations on the lower trace occur 
when the system passes through the flat spots of the 
hysteresis box. 

In the upper trace the lines drawn tangent to the 
peaks indicate the rate of decay of the oscillation and 
show that the damping in the hysteresis system varies 
with amplitude. At high amplitudes the rate of decay 
is relatively small. As the amplitude decreases the 
damping is considerably greater. For very small ampli- 
tudes the system oscillates on a line through the center 
of the box with a linear spring constant. When the 
amplitude of the oscillation is less than the height of the 
box the structural damping is zero. 

The second effect of the hysteresis on the structural 
properties of the system amounts to the free play effect 
of introducing an effectively weaker spring. At ampli- 
tudes of oscillation for which the system passes com- 
pletely through the hysteresis box the frequency of os- 
cillation is less than the frequency at low amplitudes 
where the system is linear. 

‘lutter Results with Hysteresis—The effects of this 
type of hysteresis on the flutter of a wing-control surface 
Again, the nonlinearity has 
The 
variation of flutter speed with frequency ratio for the 


system are shown in Fig. 8. 
been introduced in the control surface spring. 


saine system without hysteresis is included for ref- 
erence. The nonlinearity is again considered in a system 
corresponding to Case A (wg wa = 0.6). 

In this example it was desirable to excite the system 
in the translation, or h, degree of freedom. The analog 
results are therefore presented in the form of a ratio, 
h,/6, of the initial displacement to the width of the 
hysteresis box, against the flutter speed ratio V/V jip. 
(In this example, for 6 = 1.0°, a value of the ordinate 
h, 6 = 10 corresponds to an initial displacement in trans- 
lation equal to approximately 1/15 of a chord length.) 
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Results are shown for two values of the ratio, // 6, of 
height to width of the hysteresis box. In the two cases 
the height, 17, is held constant. The Figure at the lower 
left applies to a relatively narrow hysteresis box for 
which 7/6 = 2. Flutter occurs first at the flutter speed 
of the corresponding linear system. 
V/Vim up to 1.7 and for fairly high initial displace- 
ments the flutter is of limited amplitude. For higher 
velocities and displacements the flutter is divergent. 

No flutter occurs below V 
displacements, the system oscillates on a linear path 
through the the box. For higher dis- 
placements, at speeds below V);;,, hysteresis damping 
effects enter, and the free play lowers the spring con- 
The in- 


For values of 


center of 


stant, combining to give a higher flutter speed. 


teraction of these two effects is also responsible for the | 


region of limited amplitude flutter. For high displace- 
ments and for values of V/V, 
extends well beyond the hysteresis region. 
ing is reduced and the effective value of wg increases so 
that the flutter speed is lowered and the flutter becomes 


The damp- 


divergent. 
Results for a wider hysteresis box, /7/6 = 1, 
The results are similar 


are shown 
at the lower right of the Figure. 
to those just discussed, except that the limited ampli- 
tude flutter extends to much higher values of I’ V;,,. 
This is probably due to the fact that the wider hysteresis 
loop introduces a greater amount of damping and that 
the increased amount of free play produces a lower 


effective value of we. 


Flutter Results with the Cubic Spring 


The final type of nonlinearity to be considered is the 
cubic spring. In Fig. 9, analog results are given fora 
cubic spring in the torsional degree of freedom of a 
system free to translate and pitch. Physical properties 
of the system treated in this case are included in Table | 
the relation between moment and dis- 
140.5a 1000’, 


moment in foot-pounds and a the tor- 


For this case, 
placement was assumed to be J = 
where .\/ is the 
sional displacement in radians. 
linear term in @ represents the spring constant in foot- 
pounds per radian of the linear system. 
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Results for this case are presented in the form of the 
initial angular displacement ap (in degrees) against the 
velocity ratio V/ V;;,._ Both hard and soft cubic springs 
were considered. Flutter boundaries are shown by the 
solid curve for the soft spring and by the dashed line 
for the hard spring. In both cases the flutter region lies 
to the right of the boundary. 


tion in torsional stiffness in the linear system was de- 


In this example, a reduc- 


stabilizing. 

For the hard spring, the flutter boundary is a straight 
The soft 
spring in this case had a destabilizing effect, in that 


line at the flutter speed of the linear system. 


flutter could be induced below the flutter speed of the 
linear system by making the initial displacement suf- 
ficiently large. For the system treated here, however, 
the deviation from J’;,;, occurs only at fairly high initial 
displacements. A similar result was noted when a 
cubic spring was introduced in the control surface of the 
three degree of freedom system. 

In concluding the discussion of the cubic spring it is 


also of interest to consider the flutter amplitudes asso- 
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ciated with the results shown here. For the soft spring 
the flutter oscillation was highly divergent at any 
velocity above that of the flutter boundary. With the 
hard spring in the system, however, the fiutter ampli- 
tude is self-limited. In general, it appears that the 
effects produced by a cubic spring depend on the effects 
of a variation in stiffness in the linear system. 


CONCLUDING REMARKS 


This paper is concerned with a study of some effects 
of simplified nonlinear structural terms on the classical 
flutter of two and three degree of freedom systems. The 
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STRUCTURAL NONLINEARITIES 63 
results were obtained by the use of an electronic differ 
ential analyzer type of analog computing machine, 
which is well suited to the study of problems of this 
type. 

The form of the nonlinearities encountered on actual 
aircraft structures is in general not very well known and 
is an area worthy of further research. In the absence 
of more definite information, three relatively simple 
characteristic types of structural nonlinearities were 
chosen and studied—namely, free play, hysteresis, and 
a cubic spring. 

As would be expected for nonlinear systems, the re- 
sults are highly dependent on the amplitude of the initial 
disturbance. In many cases the flutter speed was 
lowered as the initial disturbance was increased. Some 
of the inconsistencies with regard to flutter speed which 
have appeared between airplanes of the same type may 
be due to differences in the excitations experienced by 
the different airplanes. In addition, variations in the 
extent of the nonlinearities in the different planes may 
be significant. Variations in the amount of free play, 
for example, have been shown to greatly change the 
effective stiffness. 

The existence of a preload on the nonlinear system 
has been shown to have a very marked effect on the 
flutter speed. In actual practice, the amount of preload 
on a control surface would depend on the control surface 
deflection. Thus, the flutter speed could be a function 
of the amount of control surface deflection. 

The importance of proper initial design and subse 
quent careful maintenance cannot be overemphasized. 
In one of the cases presented here, the flutter speed was 
reduced by a factor of two when free play was intro 
However, the flutter was of a 
On the 


other hand, the existence of this mild limited-amplitude 


duced into the system. 
very mild nature and limited in amplitude. 


flutter could result in wear of certain parts, and a more 
severe flutter could be encountered later in the life of the 
airplane. 

The results of this investigation show that nonlinear 
effects can influence the flutter speed significantly and 
suggest that structural nonlinearities may often be 
responsible for some of the inconsistencies which have 
occurred in the past, particularly with regard to control 


surface flutter. 
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On the Evaluation of Two Functions Occurring | 


in Maslen and Moore’s Theory of Strong 
Transverse Waves in a Circular Cylinder 


Henry E. Fettis 
Aeronautical Research Laboratory, Wright Air Development 
Center, ARDC, Wright-Patterson AFB, Ohio 


July 23, 1956 


T IS THE PURPOSE of this note to derive simple expressions for 


the functions f(a) and g(a)! defined by the equations 


(d?f da®) + (1/a) (df/da) + 4[1 (n?/a?)|f = J,7(a@) l 
(d?g/da*) + (1/a) (dg/da) + 4g = J,2(a 2 
subject to the conditions 
f=ay7=) ata =U 
df/da = dg/da = 0 ata =8 
Particular integrals of these equations are 
‘a l ex Jo,(2a@ cos 0) — Jo,(2a) 
Jn = : dé Bf 
tirJ (0 sin? @ 
(—)" (£°7 J,(2a@ cos 6) J (2a) 
dn = J aor cos (2n0) dé } 
Lor 0 sin? @ 


That Eqs. (3) and (4) are solutions of Eqs. (1) and (2) is readily 


verified by direct substitution after making use of the relation? 


1 ({°7 
Jia) = J J 2a cos 8) dé § 
T 0 


in Eq. (1) and the relation 
Ll fs 
J Ji)(2 a@ cos 6) cos (2n6) ds 6 
T 0 
in Eq. (2). 


Eqs. (3) and (4) may also be derived directly by applying La 


grange’s method of variation of parameters to Eqs. (1) and (2), 
] 


inserting the integral representation [Eqs. (5) and (6)], and 


carrying out che integrations with the aid of Lommel’s Formula 


er x 
J xJ (ax) J p(Bx) dx = ad (ax) J (Bx 
a? — p? 
ht 1. ax : 2 Bx) 7 


Expressions for Eqs. (3) 
merical calculation may be obtained as Neumann series,‘ starting 


and (4) which are suitable for nu- 


with the expansion 


Jon(2a cos 0) = J,27(a) + 2 2. Juaetaas (a) cos (2r6 8 
? 1 
which, when inserted in Eqs. (3) and (4), gives 
= ] ex 1 — cos 276 
i, = - Juche J la = d9 +] 
2rr=] . 0 sin? @ 





valu 


tabl 


taini 


plete 


wher 








l ¥ ex] cos 2r4 
d Jhla)J a) — cos 2n6 dé 10) 
2rr=1 0 sin? 6 
id, ce (see Appendix 
ex, | cos 276 om 
| eo cos 2n6 dé = 0) ifr < nf 
Jo sin? @ > 1] 
= 2r(r n ifr > a 
the final series become 
i, =--— 7 9 J a) J a 12) 


——_____ 
] 
i 
A= =) Tnx Ha 13 
r=] 
rhese series are absolutely and uniformly convergent for all 
values of m and a,® and, because of the availability of extensive 
curring tables of Bessel functions, provide a very rapid means of ob- 
Strong taining numerical values 
linder Since the functions f and @ satisfy f(0) = g(0) = 0, the com- 
plete solution is obtained as 
f = AJa,(2a) , J 
nent 
g = BJ(2a) +9 
where 4 and B are determined from the conditions f’(8) = g’(8) 
= (), and f and g are found by differentiating Eqs. (12) and (13). 
°SSiONS for 
os \PPENDIX 
2 a) l 
, x1 — cos 2r9 
9 Evaluation of : cos 2n6 dé 
7 0 sin? 6 
Px] cos 279 1 f*1 — cos 276 
| - cos 2n6 dé = J cos 2n6 + 
JI sin? 6 2 0 1 — cos @ 
1 (= 1 — cos 2r6 
cos 2n6 dé 
2/70 1 + cos 6 
*z 1 — cos 2ré 
2) | = cos 2n6 dé 
: ) 0 1 — cos 6 
From the known relation (see reference 6, formula 860.1), 
10 j 
x cos mé ma” 
| dé = = 
: , J9 | + a? + 2a cos 0 1 — a? 
is readily 
lation? man integer >0, la) < 1 
ve find 
” 
* 1 — cos m6 (1 —a™) 
d@ = 2lim | z e 
0 1 — cos @ e-!] 1 — a* 
= mr, m-20 
6 ”, 
Thus 
Pe] cos 270 1 (71 —cos2n + 1/0 
; cos 2n6 d6 = d6 + 
ying La- } 0 ] cos # 270 1 — cos @ 
9 
and (2), ] ™ 1 — cos 2(n r)0 ex 1 — cos 2né 
|], and iia 1 d6 
2 0 l cos 6 0 1 — cos 6 
rmula 
=-i[(n +r) 4+ n—r)| —2n 
y ind therefore 
C) 7 } *;7 1— cos 270 \' if r < n 
cos 26 dé = )o tr > 
0 nz Jair » r n 
for nu- e sin’ 6 on are 
starting | 
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Supersonic Flow Past an Oscillating Cascade 
with Supersonic Leading-Edge Locus 


Frank Lane 

Research Associate Professor of Aeronautical Engineering 
Daniel Guggenheim Schoo! of Aeronautics, New York 
University, New York 

August 3, 1956 


— LAPLACE TRANSFORM METHOD utilized by John Miles! in 
treating the wind-tunnel wall interference problem for the 
oscillating wing in supersonic flow may, by a slight extension, 
be employed in the study of oscillatory supersonic air forces on 
two-dimensional cascades having nonzero stagger and arbitrary 
interblade phase-lag angle, provided the locus of blade leading 
edges is not swept behind the Mach lines. The condition of 
subsonic leading-edge locus (leading-edge locus swept behind the 
Mach lines), with its attendant nonvanishing disturbance exist 
ing infinitely far upstream of every blade, would require Fourier 
rather than Laplace transform and subsequent formulation of the 
problem as an integral equation in much the same fashion as 
utilized by Miles in the subsonic wind-tunnel interference case 
The apparently restrictive assumption of identical blade modes 
and amplitudes but with equal phase lag adjacent 
blades has been shown by Lane? to be completely general with 


respect to the flutter problem, the ultimate determination of the 


between 


critical interblade phase-lag angle constituting one part of the 
solution of the blade-row critical flutter condition 

The cascade geometry is as shown in Fig. 1 Symbols are the 
same as those used by Miles, lengths being scaled with respect 
to blade chord c, the represents 
dimensionless stagger distance while o represents interblade phase 
The dimen 


and additional variable 2x 


angle, each blade leading the blade below it by o@ 


normal) gap distance 


sionless length 2y,; now represents blade 
rather than tunnel height 

The boundary-value problem is formulated for the strip domain 
shown shaded in Fig. 1. Under the system-mode shape char 


angle o, conditions at y = 2y,~* are identi 


acterized by the phase 


cal to e'® times those at y = 0* at a point shifted a distance 
2x, to the left. This connection provides sufficient data for a 


complete formulation and solution of the problem. The re 
striction to supersonic leading-edge locus ensures the vanishing 
of all disturbances upstream of all blade leading edges and hence 
for x () in the shaded strip of Fig. 1. This permits use of 
Laplace transform and leads to Miles’ equation {[Eq. (3.2)| in 
the transform ¢ (s, y) of the dimensionless disturbance potential, 
g(x, ¥) 


where A? = B2l(s + uM)? + x? 


In accordance with Miles, the dimensionless downwash is given 
The boundary 


by a(x)e™, a(s) being the transform of &@(x) 
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conditions are oo 
¢ x O) = a(x) es 
¢ x dy = e°%q (x 2x oF where 
In terms of ¢(s, y), the conditions of Eqs. (2) become j cc 
: (x,y) =L A 
- { 
o.(s, 0) = a(s) ; 
Gy(s, 2y e(te — 2x18) & (5) 
The solution of Eq 1) subject to the boundary conditions ; 
. 2) ; letting 
of Eqs. (3) is 
- io —2x1s = tii 


cosh AY cosh X(y 2y1 


and utilizing the fact th 


IN sinh 2A); 
This reduces to Miles’ result for the case of wall interference L cf 
tting o = 7.x, = 0. Byconvolution, : 
upon letting o , COR VORETIO there results, letting 
ie, f 2 2! o+2 
o (x, 0+) = Boe ~ixM ) 2 « Tn (x V x? — B%4ny:)?§.1 (x — BAny,) — Qei(et2x0M 
? 0 
. f / a 
Ly Je tk W 2x1)? B2(2y tn 
n=0 
— . ¥ we - : 
g(x, 2y,7) = Bu'e eMx <9 pa Jo vay x B2(2y, + 4ny1)?4 1 [x B(2y iny 
n 0 
eile 2xixM) ad KY x 2x,)2 
n=U 


These expressions, substituted into Eq. (5), give the dimensionless disturbance potential amplituc 


then provides sufficient data to compute oscillatory pressure-jump amplitude. 


As a consequence of the interblade phase-lag condition, the pressure jump across the lower blade of the shaded strip of Fig 


p(x, 0-, 1) — p(x, 0%, 1) = ee? 2y17) 


P(X -- 2x 


= 2p U2y( x) ¢ 


the latter relation serving to define y(x) in accordance with Miles’ notation. 


(x) = (1/2) [(0/dx) + ik] {e~*o(x 4 


2x1, 291) g(x, 


= é 2K M 
if 

f = f(r 2) 

€ = ] n=O0 

= 2 n 0 
( 

24 l|x 2x 

9 ,1 
B? ny; ct | \ 
le g(x, y) at y = 07 


e. @ 


it 
f 


Using the exponential expansion of csch 2Ay, and, for 


2 aM 
1 (x 2 
B(2y; + 40 
2x B 40H 
and y 2y 


=} 


The linearized expression for pressure then gives 


where it should be noted that the & utilized is that of reference 1, based on chord, and is twice the usual k based on semichord 


there results 


OX TEX 
— B~! 9) : { / { / 
y(x) = - + tk 7 eo alge M(x-+2m1—$)| 9 > Joyn V (x 2x i 2By, ] 
y4 Ox J0 n=0 
. o+2 / x ; 
i[x + 2x; — ¢ — 2By, (1 + 2n)] — efet2) Dg Tbe A(x — 02 — (4Bny,)2} U(x — ¢ — 4Bn 
n=0 
‘ 
: ») Ye f / * ; . 
ale )e yale (ef 2. énJo KV (x — &) $Bny,)*4 (x —¢ tBny, 
J 0 n=0 

‘ / p Ve > . a ' - 

pelo 2uxiM Jor Vv (x ¢ — 2x,)? — (2By,)2 (1 + 2n)24 1[x c 2 2By,(1 


7 

2n) 24 
V1 d 
1 On 


ie ' 
ay 


Phis 


l is given 


(10) 


1] 


Thus 


12 


The presence of the unit step functions and the supersonic leading-edge locus condition serve to cut off all contributions to the first in- 


tegral for values of ¢ greater than x, thereby permitting replacement of the upper limit (x + 2x,) 
series involving the e, become identical, and Eq. (12) may be simplified to the form 

Q Fy wy 

R i re) » § / ° 

y(x) = —B + tk ) alfje ixM(x—§ eg = > Jo VK VY (x + 2x ¢}* 1 + 

Ox J 0 n=0 

iO / " ° ot 
1 + 2n)A] +. a oe ba (2 2x ¢)? 12(1 + 2m)?51 [x 2x; 


wherein 2 and A are defined, respectively, as 


+ 2 Wx, A 


QO = = 


o 2By, 


1 Miles, John W., The Co 
Wind Tunnel, Journal of the 


In performing the 0/O0x operation, in the step involving differ- 
entiation of the integrand, derivatives of the unit step functions 


” 


by x. 
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The Warping of Triangular Wings for Minimum ; 7 
Drag at Supersonic Speeds . a 
—— Optimum, @ = Ago+ Agi ly! + Aiox 
Doris Cohen os +Aooy® + Ai xlyl + Agox® 
sh Av ]} Aeronautical Research Scientist, Ames Aeronautical Laboratory, j= ‘| jf __.. Optimum, conical camber (Ref 
( NACA, Moffett Field, Calif. : 
& August 6, 1956 
) 6 
simplicity, H. TsteEn! has investigated the minimum drag due to lift 
S. to be obtained with triangular wings of conical shape and 
6 loading. Lower drag may be obtained, however, if the restric- 
tion to conical form is removed. Using higher order homogeneous Go 
solutions? * already available,* * ® the drag of laterally sym- Ct le —o 
7 metrical triangular wings having camber and twist developable | T TT 
in the form of the double power series Flat plate 
— r —_ + _ + 
a(x, y) = p ie Ax" |y|° 7 > () 1 —— T * T 
a ae 
as been computed for VW = v2, with the following results  — | —— 2 } 
Fig. 1): “Lower bound i eae ~ 
Q Whereas—by optimizing a five-term series of a purely conical (Ref. 12) 
| family —less than 1 per cent reduction in C,/C,? from the flat- i-—_—+ T T T t 1 
plate value was achieved! for wings with m = 0.5 and 0.8, and 
81 per cent when m = 1, gains of 8.2, 6.1, and 8.9 per cent, 
respectively, can be obtained at these same values of m by using % 5) 4 3 S Te 12 
} :pproximately the same number of terms (r + s < 2) of the type ” 
y ( 9 given in Eq. (1). With additional terms (r + s = 3), the gain : + 
when m = 1 can be increased to 9.6 per cent (The figures 0 e 3 
This — cited, from both reference 1 and the present work, include the on a 
full theoretical benefit of leading-edge suction Fic. 1 Drag due to lift, triangular wings, MV = V2. 
| is given The shapes derived (Fig. 2) have generally negative camber 
; when the wing lies entirely within the Mach cone from the apex, a(x, ¥) = Aw + Aaiy) + Awe + Any? + Anx'y, + Asox® : 
(10) except for a reversal of the curvature in the vicinity of the rear <) 
third of the root section. As the leading edge approaches the As previously indicated, the gains obtained by increasing the num 
Mach cone, the region of positive camber spreads spanwise and ber of terms to ten (r + s < 3) were not very great (of the order 
11) — then forward so that, when m = 1.2, the indicated camber is for of 1 per cent), and, in the case of narrow wings, there was littl 
i Thins the most part positive. As in other investigations of the prob- change in shape 
lem of drag reduction (e.g., reference 7), washout of the angle of 
attack toward the tips is consistently indicated, particularly for 
the wider wings. _—— ym 
oid : Half wing ) oe 
It may be noted that the shape for m = 0.3 differs considerably Pion view 
ear the leading edge from that of a low-drag wing of the same oT ;* 
sweep derived by Tucker. The reason for the difference is that sie . 
the present calculations have as their basis a set of assumed een 1 
shapes, each of which, like the flat plate, causes infinite velocities . 
round the leading edge, resulting in a suction force acting in 
opposition to the drag. (On airfoils of finite thickness, an equal 8 
| force results from the very high velocities around the nose acting : \ pd 
( 12) on the small but finite frontal area of the section In Tucker’s S ” ? 
example, the loading is stipulated and in such a way as to pre- y P E . 
‘ first in- | clude any such leading-edge singularities. A substitute for the N P . 6 
, the two | leading-edge suction is therefore obtained by bending the surface ° ——— > 
downward into the stream. * . = SS 
The low value of drag indicated by the present theoretical cal- 0 " eo Wer Sakae ae 
culations can be approached in practice only if it is possible to 
maintain the type of unseparated flow assumed in the analysis 4 
The indicated decrease of angle of attack and pressure toward oy 
the tips is, of course, beneficial from the standpoint of preventing > a} y/m : 
flow separation. Nevertheless, it seems that a certain amount 0 <4 _—— 8 
t 13) of downward camber should be introduced in the immediate ~ o! —— 
{ Vicinity of the leading edge. A highly localized modification 
could be used to prevent flow separation while permitting the 
Werall lift and drag to remain as calculated “ 4 = . 
S y/m= O-——~ 
The shapes shown in Fig. 2 were all derived by combining six = 2 z = 
terms of Eq. (1 that is, = 4 = — 
foil ine  - *. 6 —_ 
pp. 671- “It should be noted, however, that the value of the drag computed by 0 2 4 6 8 - 
Tucker in his example is lower than the minimum possible value later cal eo 
sor Blade § culated by Grant? for all wings having loadings of the same form. Grant’s Fic. 2. Shape of the surface for minimum drag due to lift, 
66, Janu- — values, in fact, are essentially those for the flat plate with leading-edge M = V2. (z* is the elevation of the trailing edge and may vary 
Suctior along the span 
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Fic. 3. Lengthwise loading for minimum drag 


The lengthwise loading 


Kx) = 2 | Ap(x, y)dy (3) 


0 


is plotted for two values of m in Fig. 3. Simple incidence in the 
case of a flat plate or any conical camber, as in Tsien’s investiga- 
tion, will result in a triangular lengthwise loading. The opti- 
mum loading appears to be fuller over the forward portion of the 
wing, falling off toward the rear on the narrower wing. Only 
on the narrow wing is any considerable tendency seen toward the 
elliptical lengthwise load which is optimum in the limiting case 
of an infinitely slender wing.’ If the wing is of finite span, the 
drag depends on the load distribution as it appears in certain 
other projections as well.!! The lower bound shown in Fig. 1 is 
derived'? by assuming the impossible condition that all such pro- 
jected loadings, as well as the span loading, are elliptical over the 
entire projected length of the wing. The condition is hardest 
to approach in the case of the wide delta wing, and the minimum 
drag is consequently farthest from approaching the lower bound 
at that end of the aspect ratio range. At the lower end, how- 
ever, (m < 0.5, A 2) the lower bound is very nearly attained 
in the present calculations, explaining the apparent stationary 


state of the m = 0.3 solution. 
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Comments on ‘‘Analysis of a Rolling Pull-Out 
Maneuver’?! 


W. J. G. Pinsker 
Royal Aircraft Establishment, Bedford, England 
August 6, 1956 


r I SHE COMPUTATION OF fin loads in rolling pull-out maneuvers 
Oakes and P. T. Holliday! is clearly a 


substantial improvement on the crude estimates derived so far 


suggested by W. J 


However, the method again makes an arbitrary assumption 
that of constant angle of attack. This assumption offers no sub- 
stantial simplification of the computation nor can it be justified 
from piloting considerations, and the results may be still fairl 
inaccurate. 

In violent maneuvers, such as the rolling pull-out, pilots cannot 
be expected to discriminate more than the approximate position 
of the aircraft with reference to the horizon. Any success he 
may occasionally have in relieving normal or, for that matter, 
lateral loads by coordinated control movements can only be 
purely accidental. If, in fact, such action were feasible, the need 


for predicting such loads would be largely nonexistent. Mr 
Oakes and Mr. Holliday argue now 


elevator constitutes the more probable case, the incidence vari- 


that, although constant 


ations resulting from constant elevator have little influence on 
the magnitude of the final peak amplitude in vertical tail load 

i.e., on sideslip. This is certainly not generally true as the sta- 
bility of the lateral and longitudinal motion of a rolling aircraft 
is dramatically altered by the free interplay of inertial cross 
coupling between pitch and yaw. In fact, W. H. Phillips? has 
shown that either of these motions may become divergent at 
certain critical rolling velocities. Even before this condition is 
reached, the damping of one of the two oscillations is reduced, 
and the frequencies are substantially altered. Taking the typica 
case of an aircraft on which the period of the directional oscil 





lation is longer than the period of the high-frequency pitching | 


oscillation, the lateral oscillation in rolling will increase in period 
and lose damping, and the inverse will happen to the longitudinal 
oscillation. As a consequence, in a rolling pull out the first peak 
in sideslip will be delaved and the pitching response speeded uy 
so that large changes in incidence occur before sideslip has de 
veloped a peak value and will, of course, then have had sufficient 
time to influence the lateral motion. Also the second peak in 6 
may be considerably in excess of the first peak so that it cannot 
be considered safe to discontinue the computation until the air 
craft has settled back to a steady-flight condition. In fact, the 
maximum load is very likely to occur after the cessation of the 
aileron application, as can be seen from two typical examples of 
fully simulated rolling maneuvers in Figs. 1 and 2 

To predict aircraft motions of such complex nature from any 
simplified computations requires very scrupulous justification 
On the other hand, to represent the full pitching-moment equa 
tion and angle of attack as variable in the analysis would increas 
computational labor by not more than 20 per cent, and the gain 
would be an incomparable increase in reliability. As most de 
signers will have electronic simulators available for this work 
simplifications of the type proposed by Mr. Oakes and Mr 
Holliday seem hardly worth considering 
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Simulated rolling pull out 


On the other hand, for fast maneuvers, gravity has been found 
to have relatively little effect on the motion, and it is suggested 
that the gravity components may be omitted with very little loss 
in accuracy This will facilitate electronic simulation in par- 
ticular, where the representation of the trigonometric functions 
would be exceedingly difficult 

Using principal inertia axes instead of stability axes and drop 
ping minor aerodynamic terms, an adequate solution may be 


obtained by solving the equations 
a = QO — PB + (qaS/mUt)Craa 
8 = —R+ Pact (quad m Uo) CygB 


Q = PRU, — 1,)/1,| + (qaSl/I,) X 


[Cmax + Cyo(Ql/l 


+ Cme(ad / Uo) 


R = PQ. — 1,)/1,) + (qaSb/I,) [CugB + C,, X 


(Pb/2Uy)) + C,( Rb/21 


1 Cc 5.) 
no 4%A 


Pp = (gaSb/TIz) [Cr,(Pb/2U0) + CipB + Cis 464) 


using the initial conditions that, at ¢ = 0, 


P=P=8B=8=Q=R=a=0 


a = a = initial incidence of the principal inertia axis 
R = Ryo = (g/Uo) sin @ 

YQ = = (g Uo) (n — cos do) 

od = cos! (1/n) 
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Fig. 2. Simulated rolling pull out 
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Mathematical Analysis and Analog Simulation 
of Atmospheric Turbulence Gust Velocities 
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Control Systems Division 
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INTRODUCTION 


F  eiereswers MEASUREMENTS of atmospheric turbulence gust 
velocities by Clementson,' Douglas Aircraft Company, 
Inc., personnel,?,) NACA _ personnel,* Summers, and Chilton 
indicate that, over a considerable range of frequencies, the 
power spectral density function G(w) may be represented by a 


c/w? type of curve. For lower frequencies, G(w) is more com 


plicated, approaching a constant value as w approaches zero 


ANALYSIS 
\ power spectral density function which is able to explain 
much of the observed results in a continuous way is given by 
one whose related autocorrelation function R(r) has a simple 
damped exponential-cosine form, 


R(r) = Ae 


k 


COS Cr A2-@G, #@#>0;. «26 l 


where 4 has units of power—e.g., (ft./sec.)*, c has units of rad 


sec., and & has units of sec (see Fig. 1 


The corresponding power spectral density function is 


9 
G(w) = =f. R(r) cos wr dr 
rJ 0 


= [2Ak/x| | [w? + (Rk? + c2)|/lwot + 2(R? - eat +) 
k? + c2)2]} 


Note that for small w 
G(w) ~ 2AR lr k? + ¢? = constant 3 


while for large w 


























R(t) 
\ a. — 
b / amen 
Fic. | Autocorrelation function 
seen? seta? 
Ge) Glew) 
ie) my > _ * - 
Fic. 2. Individual power spectra 








70 JOURNAL OF THE AERONAUTIC 





| 











 % 
' / 
G(w)! — Glew \ \ 
_ | a 
Fic. 3. Typical composite power spectra 
WHITE NOISE + “ 
2 (p) 
SOURCE : se i 
4 
} 
p+b 


Fic. 4 


Generator of gust velocity power spectra 


G(w) ™ c/w? where c = 24k/r (4) 


If 3c? < k?, no maximum of G(w) occurs, and G(w) decreases 
monotonically to zero as w approaches infinity. If 3c? > k?, a 


single maximum occurs at 

w, = [—(k®? + c*) + 2c(k? +4 c?)! 2)! a (5) 
with G(w)) G(O) (see Fig. 2). The combination of two of 
these situations together, using different sets of parameters, will 
generate many composite curves (see Fig. 3). 


curves, as pictured in Figs. 2 and 3, compare favorably with 


Power spectra 
much observed data in references 1-6. Similar power spectra 
curves have been observed in radar fading and error records.’ 


ANALOG SIMULATION 


A method for generating input turbulence gust velocity power 
spectra, of the continuous form discussed above, will now be de- 
veloped. Consider white noise possessing a constant power 
If this white 
noise is passed through a linear constant parameter system, with 


spectral density function of N—i.e., G,(w) = N. 


frequency response function ¥(jw), then the output power spec- 
tral density function is 


Gow) = N|Y(jw) |? 6) 


where |} (jw) ? = square of the absolute value of V(jw). 

The desired power spectral density function has the form shown 
in Eq. (2). Equivalence between Eqs. (2) and (6) occurs if one 
constructs a system with transfer function, ¥(p) = Y(jw) where 
p = jw, given by 


) 


V(p) = [2AR/Nx]!/? {[p + (k2 + c2)!/7]/[p2 + 2kp 4 


(p) are physically realizable since no poles or 


\(p) and also Y 
zeros lie in the right half of the complex p-plane. A simple 
mechanization for Y(p) is shown in the feedback circuit of Fig. 
+ which has the transfer function 
V(p) = y(p)/x(p) = [AAalp + B)1/[p? + (a + dp 4 

ab + A,A2A3) (8) 
With one parameter in the design arbitrary [since there are five 
adjustable parameters (4;, A», A3, a, 6), while only four con- 
namely (A, k, c, N) 
between Eqs. (7) and (8) results by setting 


stants are specified in Eq. (9) , equivalence 


A\Ag = [2AR/Na]'/? b = (kh? +022) 
a +b = 2k ab + AjA2Az = k? + C2? as 


In the special case where c < k, the desired transfer function 
Eq. (7)] takes the approximate form 
> r 


V(p) = [2Ak/Nr]}/2[1/(p + b) (10 


Combinations of such circuits in parallel will generate composite 


power spectra like those shown in Fig. 3. 


' 
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Further Comments on ‘‘Supersonic 
Diffuser Instability’’ 


Charles L. Dailey 

Engineering Center, University of Southern California, Los 
Angeles 

August 13, 1956 


L. TrimptI' has taken exception to the writer’s omission? of 
e reference to NACA work on supersonic inlet instability per 
formed by Trimpi, Ferri, Nucci, Sterbentz, Evvard, and Davids 
Although this work was available in an unclassified form early 
in 1955, as Trimpi observes, it was, unfortunately, not available 
when the writer’s article was submitted (October, 1954) Se 
curity regulations also precluded mention of several other papers 
on the subject by authors at many universities and aircraft and 
engine companies in this country and in England, including work 
done by the writer as early as 1947. 

Since the NACA reports mentioned by Trimpi are now unclassi 
fied and readily available, the reader can explore the various 
opinions expressed therein as thoroughly as he wishes. The 
specific references to Trimpi’s views on buzz and to points in 
reference 2 are sufficiently pointed to require some discussion 

It is observed, first, that Trimpi’s theoretical formulation 
of the problem is physically unconvincing. He represents a 
supersonic diffuser as a one-dimensional duct of variable cross 
sectional area. This is a customary approximation, valid as 
long as the flow is stationary up to the inlet. However, when th 
subcritical shock is in front of the inlet, during the buzz cycle, 
the cone surface boundary layer of a conical spike diffuser inter 
acts with the moving shock to produce a separation strong enough 
to by-pass the inlet entirely. At this moment, air, which has 
already entered the inlet, reverses and spills out around the cowl 
lip. As this continues and the high-pressure air stored in the 
plenum chamber empties out, a relatively high-frequency oscilla 
tion of the expelled shock is observed, due apparently to the by 
passed air flowing obliquely over the sharp cowl lip. Thes¢ 
points have been observed and described in detail in reference 2 
and are regarded as part of the basic explanation of the buzz 
phenomenon. They clearly cannot be explained by use of a one 
dimensional theory. 

To avoid this difficulty in his calculation, Trimpi simply uses 
experimentally determined pressures for the part of the cycle in 
which the shock is expelled and the characteristic relations of 
constant area channel flow for the ‘‘fill-up’’ phase of the cycle 
In this way, he uses the experimental ‘‘answer”’ to fill in the inter- 
val over which his theory is not valid and employs the theory of 
constant-area channel flow where it does apply to calculate a 
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fairly accurate variation of pressure with time as the plenum 


chamber fills up due to the excess mass inflow in the supercritical 


phase of the cycle 

It is important to distinguish two things here: (1) his use of 
the measured pressure.time history when the shock is expelled 
is artificial in terms of his one-dimensional theory; (2) his use of 


the one-dimensional theory to describe the pressure-time history 
after the artificial starting condition has been imposed in the 
fill-up phase of the relaxation cycle is sound. Trimpi’s models 
had sufficiently long plenum chambers that the simple constant- 
area theory could be used For more realistic cases, the tedious 
method of step by-ste p construction would have to be employ ed. 

It should be mentioned that the writer was aware of Trimpi’s 
work at the time reference 2 was written, although reference to it 
could not be made at that time because of security regulations 
In fact, the reference on page 748 of reference 2 to the use of the 
characteristic method for the fill-up phase was made with Trimpi’s 
excellent result in mind. As mentioned there, it was not felt 
that the additional accuracy afforded by this method was neces 
sary in establishing the essential relaxation aspect of the buzz 
cycle 

It is the writer’s impression that Trimpi’s work is an unusual 
mixture of good analytical judgment and inaccurate physical 
description. For example, the statements! concerning (1) the 
thought that a buzz cycle can be generated by the spontaneous 
“outward”’ motion of a shock in a one-dimensional channel and 
(2) that the subsequent wave patterns can result in reflected 
compressions which, by the simple mechanism of reflection, are 
somehow able to pump the plenum chamber back up to its initial 
pressure are not entirely convincing. 

For the first point to be true, Trimpi’s one-dimensional channel 
would have to have a pulsing area in the oncoming ‘‘free’’ stream 
which would permit it to adjust to the instantaneous mass flow 
required when it reached the inlet. In fact, the stream tube 
would have to pulse through negative values to permit outflow 
of air which enters the inlet just before the buzz cycle starts 
It is not clear what laws of nature would have to be invoked to 
permit such a behavior. It is difficult to find a meaning for the 
second point mentioned above in view of Trimpi’s apparent in- 
sistence that the pressure rise in the plenum chamber is not due 
to the excess supercritical mass inflow over the outflow corre- 
sponding to a subcritical exit throttle setting—i.e., a relax- 
This same point appears again in the dis- 
He apparently 


ation phenomenon 
cussion of the exit throttle effect on frequency 
does not recognize that the excess inflow during fill-up increases 
as the exit flow decreases with a decreasing throttle area. 

One point which Trimpi misses completely is that buzz con- 
sists of a sequence of individual relaxation cycles. Although they 
occur at almost periodic intervals for average mass flows below the 
transition point, the individual nature of the cycles and the ran- 
dom time interval between cycles is clearly demonstrated for oper- 
ation very near the transition point. This is mentioned on page 
737 of reference 2 and is shown in Fig. 9 of that reference. It is 
obvious from this consideration alone that the theory of duct 
resonance employed by Trimpi is qualitatively incorrect. As 
mentioned ahove, the only validity it has is in prediction of the 
pressure-time history of the fill-up phase of the buzz cycle. 

Trimpi objects to the writer’s omission of upstream traveling 
compression waves in his consideration of the possibility that 
plenum-chamber separation might cause inlet instability It 
was not intended that those qualitative considerations serve 
any purpose other than to show that plenum-chamber effects must 
be unimportant and that attention must be focused on the inlet 
region to find the cause of instability. It is interesting that 
Trimpi attaches such importance to this point immediately 
after mentioning an instability theory of his own concerning 
conditions in ‘‘the supersonic diffuser and immediate neighbor- 
ing subsonic diffuser’? and also mentions work of Ferri and Nucci 
concerning ‘‘the initial rate of subsonic diffusion,’ which seems 
to indicate that he also feels the inlet region is important with 
respect to stability. 

Trimpi also objects to the writer’s statement that instability 


FORUM 71 


is a result of choking of the inlet flow. Experimental work at the 
U.S.C. Engineering Center currently being conducted by A. C 
Brown, recently of the Bristol Aeroplane Company, and th 


by the Office of 


writer under a research contract supported 
Scientific Research is shedding light on this aspect of the insta 
bility problem. This work is broader than that of reference 2 
in that it includes the ramp and scoop type as well as the spik 


type It is unclassified, and the results will be published soon 
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Erratum—‘‘Note on the Unsteady Heat 
Regenerator’’* 
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September 18, 1956 


ips AUTHORS are indebted to H. Riedel for drawing their 
attention to an error in Eq. (2 rhe correct equation is 


hs oT /ot t Cy O pl ol T Cy O pul On ) 


and leads to the result that there is no compressibility correc 
tion to the total mass of gas available at the outlet at constant 
temperature when the pressure is kept constant 


* Meyer, R. E., and d’Ews Thomson, T. A., Readers’ Forum, Journal 
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On a Third Mode of Longitudinal Control- 
Fixed Oscillation 


Bernard Etkin 

Research Associate and Associate Professor of Aeronautica 
Engineering, Institute of Aerophysics, University of Toronto, 
Toronto, Canada 

August 14, 1956 


SYMBOLS 


I period 


time to damp to half amplitude 


u speed perturbation 
a angle of attack perturbation 
6 angle of pitch perturbation 
Cr static stability derivative, o¢ Jae 
a 
DISCUSSION 
Bie SHORT-PERIOD and phugoid modes of longitudinal oscilla 
tion of airplanes have been known for a long time Phese 
are characteristic of aircraft over a wide range of configuration 


and flight variables. If the variations in uv, 6, and @ be depicted 
by rotating vectors, these modes appear typi illy as in Figs. 1 
and 2. Fig. 1 shows the phugoid mode The important char 
acteristics are the phase difference of about 90° between u and @ 
and the nearly zero value of the a perturbation rhe short- 
period mode, shown in Fig. 2, is characterized by the nearly zero 
value of the speed variation 
Some recent calculations have revealed what appears to be 

new mode of oscillation, fundamentally different from the two 


classical ones. The appearance of this mode is best shown by a 
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Fic. 1. Typical vector diagram of phugoid mode. 


plot of the loci of the roots of the stability quartic for varying 
c.g. position. Such a plot is illustrated (not to scale) in Fig. 3 

The calculations described were for a conventional configura 
derivatives 


values of the stability gross 


altitude = 30,000 ft., speed = 500 m.p.h 


tion with typical 
weight HW’ = 100,000 Ibs., 

For forward c.g. positions, the conjugate roots giving the short 
period mode are at D,D', and those giving the phugoid are at 
E,E'. As the c.g. moves aft, 
reals at A and B. Point A corresponds to a c.g. position 1 per 
cent of the M A.C. ahead of the stick-fixed neutral point, and point 
B corresponds to the ¢.g. almost on the neutral point. As the 
c.g. moves still farther aft, the complex roots disappear and are 
Of the four real roots, 


these roots approach the axis of 


replaced by pairs of real roots as shown 
three are negative, corresponding to subsidences, and one, on 
the branch BG, is positive, indicating a static instability 

One branch of the locus coming from A and a second branch 


coming from B meet at C. This is at a c.g. position between 
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Fic. 3. Variation of roots with c.g. position. 
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Fic. 2. Typical vector diagram of short-period mode 


2 and 38 per cent aft of the neutral point. From C there springs 


a new pair of complex branches, CF and CF', which define a 
mode, essentially different from the phugoid 


Analysis of this mode with the c.g. 3 per cent 


third oscillatory 
and short-period 
aft of the neutral point has shown it to be of long period and 
heavily damped: 

i = 76.5 sec 


thay = 0.65 see 


(For more aft c.g. positions, the damping decreases rapidly 
The vector diagram of Fig 
It is seen to be entirely different from 


$ shows the relationships among the 


variables in this mode 
those of Figs. 1 and 2; in particular, all three degrees of freedom 
are significantly excited. No two degree of freedom approxima 
tion to this mode is possible 

between the phugoid and 


“cross” 


In a sense, the mode is a 
short-period. Its period is comparable to that of the phugoid, 
but the number of cycles to damp is characteristic of the short 


period mode. The long period is responsible for the significant 
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inges, and the positive Cma which exists with the c.g. aft 


neutral point is the cause of the large variations in a 
One of the problems of airplane design is to determine the per 
issible aft limit of the c.g The solution of this problem requires 
lerstanding of the behavior of the airplane for unstable as 
for stable positions of the c.g rhis is particularly so 


en the air-frame dynamics are neutrally stable, or even un 


auto 


nd acceptable flying qualities are provided by an 
The investigation reported here throws additional 


of the neutral 


tic pilot 


ht on the dvnamic behavior when the c.g. is aft 


The Effect of an Arbitrary Distribution of 
Surface Temperatures on Heat Transfer in 
the Compressible Laminar Boundary Layer 

Craig M. Schmidt and Arnold J. Hanawalt 

Engineers, Bell Aircraft Corporation 

195¢ 


Buffalo, N.Y 


Aerodynamics 


ootember 2 


[: IS THE PURPOSE of this note to discuss the effect of a realistic 
t-transfer coefficient 


wall-temperature distribution on the he: 
The examples that have been chosen by Chapman and Rubesin 
heat 


ind by Lighthill* point out a very marked change in the 


transfer coefficient However, it is felt that these examples ar« 
not very realistic for many high-speed flight conditions. Gen 


erally, in supersonic or hypersonic flight, the leading-edge tem 


perature is below the recovery or adiabatic wall temperature, anc 
the wall-temperature profile is such that it is ste« p at the leading 
edge and becomes more flat as the trailing edge is approached 
rhe following consists of four examples that tend to point out 
the degree of accuracy attained by the constant wall-temperature 
theory The first assumes a cubic wall-temperature variation 
ind is analyzed by the theory of Chapman and Rubesin. Since 
this theory is limited to a power series expression for the wall 
examples use the theory of 


temperature, the remaining three 


Lighthill, modified to include compressibility effects. Equilib 
flight 
while, for the last case, the wall temperature is assumed to vary 
The S¢ 


point out that, for many design purposes, it is not 


rium temperatures are calculated for two conditions 


is the square root of the distance examples tend to 
necessary 
to use the more complex variable wall-temperature theory 

For the first example, using the method of Chapman and 
Rubesin and their notation, let us consider the following, which is 
representative of the temperatures associated with high-speed, 


igh-altitude flight 


i $00°R 7. = 6,000°R 
t(x* 7 — 7,* = T,*(—0.667 — 0.646x* + 
0.792 x*? — 0.3125 x* 
T.* = 2:94 C = 0.758 


The ratio of variable wall temperature g and f to their counter 
parts for a constant wall temperature are given in the following 


equations: 


ob — (172°; 0)/T. 
At un fi — 
. £~— (T/T, V’(0 
V'(1)x* + ao ¥"(2)x*? + a3 ¥'(3)X*5] 
1 — (T,*/T,* { 


rhe results of these computations are shown in Fig. 1 It can 
be seen that the variation of 4/h is not particularly large, show 


ing a decrease (from the constant wall-temperature theory) at 


the leading edge of 16 per cent and a maximum increase of about 
12 per cent occurring at about x* 0.6 

Introducing a higher degree polynominal than in the previous 
example in order to obtain a steeper gradient at the leading edge 
followed by a relatively constant temperature shows the effect of 
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characteristics for square-root variation of 7, 


Fic. 3 


the temperature gradient on the heat-transfer coefficient to be 


even smaller. This also is generally what has been found using 
Lighthill’s method, and illustrations of this are listed below 
Lighthill? has shown the effect of variable wall temperature on 
equilibrium temperatures (i.e., convective heat in equals radi 
ative heat out) to be very significant for the several examples 
chosen. However, the physical assumptions used in this com 
parison are inconsistent—namely, that the surface conductivity 
for the variable wall-temperature case was zero and for the con 
stant wall-temperature case was infinite. If the conductivity 
is assumed to be zero for both cases, then the results are much 
different than what Lighthill has shown, as ean be seen for the 
two examples shown in Fig. 2. These two examples show that 
the simpler constant wall-temperature theory is accurate to 
within 3 per cent for the calculation of equilibrium temperature 
Lighthill also states that his theory agrees to within 3 per 
cent with the more exact theory of Chapman and Rubesin 
This is only true if his expression for heat transfer is altered by 
factor C 
account for the evaluation of the physical properties of the air 
The addition of this factor, Lighthill's 


the equilibrium temperatures and neces 


the multiplicative VC (notation of reference 1) to 


however, invalidates 
analytic solution for 
However, 


sitates the use of an iteration method of solution 


since Cy/7/C is quite close to 1.0 and since it enters approxi 


mately to the 1/4 power for the equilibrium temperatures, the 
results shown in Fig. 2 are felt to be essentially correct 


\s an 


perature on heat transfer, the following case was examined: 


Tw = (1/3)T. 


idditional example of the effect of variable wall tem- 


(4/6)7. WX 2, 100°R., 7, 6,000°R. 


This is felt to be somewhat representative of accelerating or 


quasi-steady-state flight at high speeds. The results are shown 


in Fig. 38. A maximum of 16 per cent increase over constant 


wall temperature theory is observed. As the leading-edge wall 
temperature approaches 7,, this increase becomes larger and may 


reach a value of 50 per cent for a square-root variation of 7, 
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However, 7 ‘y 


istic except possibly for low supersonic speeds 


it the leading edge is not felt to be very 
Examples given by Chapman and Rubesin and by Light 
tend to overemphasize the importance of variable wall temper 


ature on heat transfer while more realistic examples show 


heat flux to be less dependent on the space history of the 
temperature In a good many design cases, the use of the sir 
constant wall-temperature theory would introduce error 
greater than about 10 to 20 per cent \s far as equilibrium 


temperatures are concerned, the error would be much less t 


this. However, in some cases such as wind-tunnel tests wl 
the temperature levels are much different than that experien 
temperature gradients 


in free flight or in local areas of large 


(e.g., ends of fuel tanks), the variable wall temperature would 


be much more important 
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A Note on Missile Launching 


Telford W. Oswald 
Engineering Staff As 
Calif 


A ugust 21 


stant, Hughes Aircraft Co ilver y, 


1956 
: A RECENT PAPER,! a method of computing launching traje 
tories for the case of a missile 


velocity has been presented 


whose natural frequency is p1 


portional to In practice, such mis 


siles are char icterized bv small ch inges 1n inertia and stability 


margin. The following analysis extends the method to account 


for large changes in inertia and stability margin such as occur 


with a solid-propellent-rocket-powered missile launched at sub 
sonic spec ds and accelerated to supersonic speeds Such missiles 
experience simultaneously a forward shift of the center of gravity, 
inges in 


an aft shift of center of pressure, together with large ch 


the lift curve slope 
The outlined method provides an approximate solution which 


is readily obtained by quadr iture It affords a large time saving 
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Fic. 2 
ver the method of numerical integration of the differential 
juations of motion and avoids the troublesome question of con 


It applies equally well to missiles which are unstable 
t launch, providing the stability increases with velocity. Typi 


presented compared with “exact”? numerical 


that 


examples are 
These 


d approximation 


tegrations indicate the method provides a very 


Phe longitudinal equations of motion (short-period oscillation 


re referred to a set of principal body axes Phe longitudinal 


tion is assumed to be a uniform acceleration rhe damping 


] 


ivatives are assumed to be proportional to velocit ind the 


queney in pitch is approximated by the expression 


ra A?/U l 


validity of these latter two assumptions is shown in Fig. 1 
Under these conditions the equations of motion are combined 
ir differential equation for the angle of 
With the 


is iIndepe ident vari ible, the equation is 


me second-order line 
ick pec 


ng the flight path, x, 


reference 1 for derivation distance 


angle of attack w/ 


Vv —UM,' + Z0'M, 


f = forcing function 


pitch frequen 


become 


Eq. (2 


ind & = wox, 


[ T 1? a* F l é 
n 1 y = 

dé iL $e? 2a 
rhe homogeneous equation obtained from E 


tion of ordet where 


1/2) V1 + (A2/a? 1 


solution to Eq. (8) in the case 





( mplete 
eger* 1s given by 


y = ¢ T (3 4 C.J 


ere the particular solution, / may be obtained by the 


ethod of variation of parameters 


*If » is an integer, J must be replaced by } Bessel function of the 


econd kind 


FORUM re 


The solution for the angle of attack for the se of inching 
with a finite velocity UL’; is given b 
a exp CiJ ( l t 
where Ww 2a ist value I tt laun eer ( 
( nd ( re ¢ tants t ve evaluated f the init 
tion 
\ exp (¢ y ( ( ‘ 


exp (i ; j 
x Z vii ‘ 
CJ Csi w \ 
and /(é, &) is determined from the equation 
. 
a 
I(é, 3 exp Cod J 
2 sl i 
v7 § 
ve » 
J (i : ( ) 
, oI ~ 
« é - J 
where 
= lexp (i 2au Vr/1 6)/1 
}(Z(6)/m] +4 T sin €)/n cos O}f 10 
Since it is assumed that V 2 where iS ¢ 
stant, the forcing function, f(£), can be determine i function 
of 
wox = ay l/?/2 1] 
In this form the equation provides for any integrable forcing 


numerical integration, Eq. (9) can 
Combining this with Eq. (6), the angle of attack 
With the angle of attack determined, 


determined from. the equation 


function, and, by means of 
be evaluated 
i function of € is known 


the flight-path tangent can be 


rhe trajectory can then be determined from the equati 


| 
| 
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In applying this method of calculation it is first necessary to 


determine vy from Eq. (4), after approximating the missile sta 


bility by Eq. (1 In doing this, it is suggested that the value of 


v nN (1/2), where » is an integer, be selected and the corre 
sponding values of 4A computed from Eq. (4) and those of w,,, from 
Eq. (1 Use of this value is permissible since the Bessel func 


the order + It is desirable 


(1/2) are 


tions are continuous functions of 


because the Bessel function of order + n tabulated 


in references 2 and 3 

Ty pical launch cases for a repres¢ ntative missile have been com 
puted, and the accuracy of the basic approximations is shown 
in Fig. 1 The results are presented in Fig. 2 in the form of flight 
path angle versus distance for a subsonic launch at high altitude 


(case 1) and a supersonic launch at low altitude (case 2 These 


results a step by step 


are compared with those obtained from 
integration of the differential equations of motion 


In order to indicate how sensitive the method is to the accu 


racy of the approximation, the solution for an initial angular ve 
locity is computed using two different approximations to the 


stability velocity relation. These two approximations and the 


results are shown in Fig. 3 
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On Aeroelastic Solutions 


William J. Gaugh* and Joseph K. Slap** 
Northrop Aircraft, Inc., Hawthorne, Calif 


August 27, 1956 


_ ERENCE | discusses a procedure for estimating aeroelastic 
effects on aircraft lateral stability and control. This paper 
is one of a number in the literature which, more or less rigorously, 
Cer 


tain remarks pertinent to such methods might be warranted at 


derive short-cut methods for aeroelasticity calculations 
this time 

In these procedures, simplifying assumptions must, of course, 
be made. Among the more common are aerodynamic strip 
theory,? variation of stiffness parameters with fourth-power of 
the chord,* and known shape of the deflection modes.? These 
assumptions usually lead to equations which permit isolation of 
the parameters being varied in the preliminary design study 
Thus, the effect of varying one or more of the parameters can be 
analyzed. Although the trends predicted are usually correct and 
of value, the basic air-load and stiffness assumptions introduce a 
lack of confidence in the accuracy of the estimated aeroelastic 
effects for any one particular model or condition. It is therefore 
expedient to use the ‘‘quick-and-dirty”’ approach for getting the 
design into “the right ball park’’ but not for fixing on a precise 
configuration 

The advent of high-speed electronic calculators has provided 
a valuable tool for detailed engineering work. This is particu 
larly true in quasi-stationary aeroelasticity. Once a method of 
solution has been selected, it can be programed for automatic 
stiffness influence coefficient ma 


calculation If air-load and 


* Aerodynamics Engineering Specialist 


** Group Engineer, Electronic Computing and Numerical Analysis, 


Fluid Dynamics Department. 
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trices are used, their size should definitely be variable in the com 


puter program. By allowing the number of span stations 


specified, the designer can obtain solutions, at minimum cost 
which 


For example, if a large 


accuracies at the particular design stage 


number of configurations or conditions is to be perused and onl 


crude estimates of wing stiffnesses are available, then perhaps as 


few as five semispan stations should be considered In an ad 


vanced design, this number may be 20 For the smaller ma 


trices, more than 60 conditions have been calculated? during an 


S-hour night shift. It is suggested that, for a few of these cases 


10- or 20-point solutions be computed in order to estimate the 


As a rule of tl 


effect 





f loosely spacing the span stations 





doubling the number of stations approximately triples the 


puting time. With the above in mind it is now proposed that, in 


the preliminary phase, the range of possible configurations be nar 


rowed by use of approximation techniques, such as are recom 


mended by references 2-4. The optimum preliminary design 


should then be chosen using a precise theory set up for high 


speed calculatior As stated previously, a significant amount 


of computing time can be saved here by use of fewer calculation 
points than would be warranted in advanced design stages 


If a configuration exhibits small aeroelastic effects, then a rough 


estimate will probably not have serious consequences. But if 


severe losses, such as 50 per cent of Cc; or near-reversal of C 
a 


are contemplated, then it is essential that the utmost confidence 
in the analvsis be guaranteed 


Almost always, an approximate method is compared with 


more refined one for justification purposes. Frequently, however, 


as in reference 1, the comparison is made using either an extremel 


simple configuration or one having a high stiffness level. In 


either case, fair to good agreement can be expected sut extrap 


olation of the same technique to complex or moderately flexible 


models should be executed with caution Substantiation of ip 


proximating assumptions must be made via a wide variety of 


examples. If this is not done, the limitations must be stated 


precisely in order to avoid misuse 


In reference 1, for example, al! limitations of the theory are 


not apparent It is perhaps illuminating to investigate thesé 
Eq. (2) of that reference expresses the wing tip twist angle 
an infinite series of twists: 

€ € € € 2 


] 


where ¢ is the twist due to the initially applied loads and each 
successive €;+; Is caused by the loads resulting from ¢ Eq. (4 
of the reference assumes that the ratio of successive twist angles 
isaconstant: 
€2/€; = €:/@ = €;/€ = R } 
Phe conclusion is then drawn that 
e = ¢€,/(1 k 5 


The author neglects to discuss & other than to state ‘‘for conven 
tional wings with sweptback elastic axes this term is usually posi 
tive and of order unity so that an error in 1/(1 k) is small com- 
pared to an error in k.”’ Actually, & is of prime importance and 
is equal to the ratio g , Where g is the flight dynamic pres 
sure and gyi, is the divergence dynamic pressure. If k < } 


then not only is Eq. (5) invalid but the series in Eq. (2) does not 


converge and the wing in fact diverges physically Eq. (4 
however, gives a finite answer for «€ even when g > > 0 It is 
precisely here, in obtaining an unconsciously incorrect result, 
that the danger lies in applying an equation without a clear ex 


ae > i, 
cannot represent Eq. (2) for certain formula 
That is, the 


position of its limitations then the wing is physically 
stable, but Eq. (5 
tions of the problem since the series again diverges 
from ¢; must have certain 


operator which yields e properties 


in order for Eq. (5) to hold if g > |g »q 0. This can im 
pose a severe restriction because, even if g,;, were a mathematical 
+ These figures apply for the IBM Type 776 Electronic Calculator, using 


a matrix formulation of the theory of reference 5 as outlined in reference 6 


ire consistent with his estimate of the air-load and stiffness 
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in the com jmit on g, it would by no means be a physical limit for lighly Continuity also requires that 
tions to be swept wing his series convergence problem is discussed in 
; a j : ; , (0 O Op:/Ol + V-(¢ 

mum cost, more detail in reference 5. It 1s adequate here to emphasize that : q 
nd stiffness the foundation of a simplified theory may be invalidated by re Now, in order that the displacement effect of the surface ; 
e, if a large ns not uncovered because of insufficient stud Reference nonviscous flow be the same as that actually produced bv the 
d and onl unfortunately selected only a nondivergent unswept wing for boundary layer, one must impose the conditi 

perhaps as rrelation purposes. Consequently, the inherent limitations were 

In an ad not encountered 0 
maller ma , = , : 
, py Eqs. (2), (4), and (5), when substituted in Ec t d, after 
r during an REFERENCES Re 

simplification, 

these « iSes, Rodden, W im P 1 1 / / Y 
timate the : / bilif he I 11 Stat ind ( 1 Journa 

of thumb {the Aeronautical Science Vol. 23, No. 7, pp. 660-662, July, 1956 A ng A nig pq) dy | 
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sre ‘TeCOND HW irbit Plan Form and Stiffness, NACA TN 1876, 1949 

r\ ] 4 He epeth, John M., Waner, Paul G., and Kell, Robert J iS fied wherein / is taken to be infinite for purpos« f integratior Cl 

iry design : =~ 

Vv ( 1 { frect he Rol raft, N N : 
, 7 ‘ ; i ‘ ACA 7 the boundary layer 

» for high 37) 1055 
nt amount Gaugh, William J., and Slap, Joseph K., Detern iio I 1 Wins Eq 7) is the defining relation for the displacemen urface ot 
calculation A¢ Paramete Journal of the Aeronautical Sciences, Vol. 19, No i three-dimensional unsteady boundary laver For steady flow 

= 3 73-182, March, 1952 the second term of Eq. (7) is zero, and the previous result! re 
iges : 

. Ernest A WwW 4M VW ; 1 f Effect } Ia lime? ] 
’ ol . mains; and, further, if the boundary laver is two-dimensional, 
en a rough Northrop Aircraft Report No. NAI-56-249, March, 195¢ 


‘'s. But if the appropriate solution of Eq. (7) is 


rsal of C, ; 
, + - 9% Pa J, 
confidence A=6 (1 - ) ay x 


ed witha | Displacement Thickness of the Unsteady which is the usual two-dimensional result 
, however, Boundary Layer The second term of Eq. (7) relates to the piston-like effect of 


extremel aerodynamic heating discussed, for example, by Van Dvyke 
level In Frank n K Moore* and § mon ( strach ** Van Dvke considered the case of a doubly infinite flat plat 
ut extrap Cornell Aeronautical Laboratory, Inc., and Lewis Flight Propulsior suddenly set in motion at a sufficiently high Mach Number s 
lv flexible Laboratory, NACA, Respectively that pressure waves are generated by the boundary-layer growth 
ion of ap August 29, 195<¢ For this problem, the boundary layer depends only on 4 rhus, 
variety of the first term of Eq. (7) vanishes, and a suitable solution 1 
be stated 5 he PURPOSE of this note is to extend a previous result! for R 
the displacement effect of a three-dimensional boundary A= (: : ) d\ q 
heory are laver to apply to the unsteady boundary layer and to apply the J0 , 
ate these result to certain flows for which unsteady boundary-layer solu a . 
which may be evaluated from Van Dyke's solution for the profil 
t angle as tol ire available j 
; jf , : of p/p; to vield 
rhe general definition of displacement thickness may be stated 
s follows At the outer limit of the boundary laver (savy at A V 2vit/n 1)? 10 


= h, where y is measured normal to the surface) a certain dis 
where v is kinematic viscosity, JJ is Mach number, and is the 


ind each tribution of normal velocity v appears owing to the displacement 
5 } . Oo ots . ats ‘ (10 onto s to ) < sult: 

Eq. (4 efect exerted bv the boundarv laver: if the flow were entirel ratio of specific heat Eq. (10) conforms to van Dyke's result 

ee : ‘ ‘ Sas ‘ted as a dis nt tl 1ess eff, 
ist angles nonviscous, the same distribution of v at y = h could be pr when interpreted as a displacement thickness effect 

; , , . O ns s« ons te hich Eq. (7) mav 1 simply 

luced by the presence of an impermeable surface 4 A (x, z, t ther unsteady solutions to which Eq. (7) may be simp 

where - are coordinates in the surface. This fictitious dis applied are available* *— namely, those for the compressible lami 

‘ ; “ ; 1 — c on . 1 . . 1 
' placement surface y = A is useful for determining the effect of the nar boundary layer of a semi-infinite insulated or isothermal flat 
1 , x . . ] S1-S lo y so ons Oo 

boundary layer on the outer flow. It is required to determine plate in nearly quasi-steady flow Phe solutions were obtained 
is expansions in the various time derivatives of a differentiable 


A(x, s, f) in terms of the boundary-layer profiles which are pre 
5 but otherwise arbitrary flight velocity | 





a sumed to be known 
. . , ; hs . ae « ‘ny } necsel j 
In the fictitious nonviscous flow (subscript o.f.), we have the Insulated Plat Phe required quantities for the isulated 
r conven , : : ; 
} boundary condition at the displacement surface plate are 


1ally posi | 





mall com | i , = qgi‘grad A + OA/Ol ] aa (1 0g ) 1794 xCmr ¢ 
0 = ay (oi )1 
wasup and fF Consequently, at y = h, if the boundary layer is thin, Ja pigi/ Vi 
mic pres- | é —_ i ve a 
< o , = Qi:grad A + OA/Ot + (h A) (0 oy 2 0.645) ( 1) Af? 0.592) |1 + (2.938 1)M 
does not Where g is the projection of the velocity vector on the surface vl''/U?) + (0.694) [1 + (2.702 1) MW ‘ 
¢ 5 nd ‘TI 1 note e Te) oO t y of t] . ‘ ) 
Eq. (9), and subscripts ] and 0 denote evaluation at the outer edge of the 0.445) [1 (1.361 : iy er ] F 1 
0. tts boundary layer and the surface, respectively 
. 
ct result, In the boundary layer (subscript 6./.), continuity requires that ( 1 f ) as \ay rCM 671 1410 2 992 
a = i Mi ) ~- é 
clear ex = ) m y 
; O( pr’; ov = Op/dt — V-( pig 3 /J0 , \ 

yhysically | 

; : » ( j"2 1 (2 9% 2] i j 1 1) yy2)2 4 i ) 
formula ntegrating, one finds, after rearrangement of terms, ruse 3.229) (x?l l 1.042) (xl + (12 
jat is, tne p h = h\{Op,/ot + V-( pig + where x is measured aft from the leading edge and C is a constant 
SOpere 2) °) required to adjust an assumed linear viscosity relation 
+ conan v- rig pq) dy + p p)dy } If it is assumed that the solution for A has the same form as 
ra “0 Of Jo that of Eq. (11), with different numerical coefficients, then sub 
ator, using Head, Aerodynamic Research Department stitution of Eqs 11) and (12) into Eq. (7) and collection of lik« 
erence 6 hicf, Applied Mechanics Branch terms yields the following result for A 
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vU'/ U2) + 0.310 [1 (3.322 1) M?| (x#U"/l 
0.776 [1 + (2.056) ( 1) A?| (xl U2)? + ; 13 


Eq. (13) is written in such a way that the significant differences 
between the true displacement thickness A and the thickness 6* 
defined in the usual way by Eq. (8) are immediately evident 


Isothermal Plate 


a 4 4 
+ 0.167 (5 1)M? + 1.126 ( ) 
Viel ’ 


0.592] 1 + 0.486 1) 1s? 


For the isothermal plate,’ 


0.694 1 1 + 0.259 (- 1 ).\/? 


[ d “# 
soa 0.0280 1)1\/2 4 0.128 ) 
49 


where 9 denotes temperature and the subscript w denotes surface 


conditions, and 


(: : ) ay > ne & 1)? [0.288 
Jo p Vi 
0.444 (x l"/U? 0.3095 (xl'"/l 0.0713(xl U2)? 
2 |(0,, A,)/A,) 10.969 0.0443 (xU'/l 0.316(U7" /l 
= , >\e i se 
0.1603(x/ (°)? + \ 15 


Substitution of Eqs. (14) and (15) into Eq. (7) vields the displace 


ment surface for the unsteady motion of an isothermal plate 


A = 6* + 0.5735 VV xCy,/l } 1 + 0.669 (5 1) jl? 
(xl''/U2) + 0.310 [1 0.863 ( 1)M?2] (xl’"/L 
0.776 [1 0.0644 (y — 1)IP] (xU'/L { 16 


Eq. (16) shows that the difference between A and 6* for an iso- 


thermal plate does not depend on wall temperature 
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On Some Concepts of Nonlinear Oscillation 
Theory 


K. Klotter 
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Stanford 


IT A RECENT PAPER in this Journal,! F. H. Clauser has set him 
self the task of promoting an understanding of nonlinear 
phenomena by presenting them ‘in terms of fundamental and 
well-understood concepts.’ He may be credited with having 
accomplished his difficult task, in general, very well; some of 


the remarks in the paper are truly elucidating 
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' 
In a field which is still in rapid evolution and where concept 
in many instances have not vet been formulated to complet | 
satisfaction, any author of an expository paper is burdenc 
the two responsibilities of helping such an evolutior 
not doing anything which endangers that process. The | 
feels that in a few cases the author of the paper quoted vi t 
stipulation (b) because of either improper generalization 
fusing terminology | 


rhree items will be singled out for discussi 


1) The phenomenon of the sudden changes in tl 


of forced oscillations (Fig. 24 of the paper), the so-called 
phenomenon,” is linked to and called ‘hysteresis’ by the aut 
is unfortunately has been done by other authors before I 


similarity between a true hysteresis loop (e.g., in a displacement 
force plane for elastic hysteresis or in a B-H-plane for magnet 
hysteresis) and the “loop” in an amplitude-frequency plane is } 
purely superficial, however, because (a) the latter lacks t 
essential property of a true hysteresis loop (as defined in 
field of physics) that the loop area represents the energy dissi 
pated; (b) in response diagrams (of the amplitude-frequen 
plane) even if jumps do occur they need not necessarily lead t 
closed ‘‘loops’’—e.g., in the case of an undamped system; (c) phys 

ical hysteresis curves, on the other hand, do not show any | R 


jumps at all Hence, the term 





sudden changes—1i.e 
teresis’’ ought to be carefully avoided in describing the 
phenomenon. Conversely, using it certainly does not help 


clarify the concepts 





2) Clauser uses the important terms “autonomous” and ‘‘n 
autonomous” in strict contrast to the established usage H 
starts from the generally accepted criterion (which he dos ov 
cite literally but rather in a paraphrased form) that if the inde 9 





pendent variable ¢ does not appear explicitly in the differenti 


equation, the system is called autonomous, and that if ¢ d 


appear explicitly, the system is called nonautonomous het 
however, Clauser goes on to apply the criterion to the 
side of the differential equation exclusively. Consequently, he 


calls all systems autonomous except those which have varying 





coefficients. Cases where a driving function appears on the right the 
hand side of the equation he calls ‘‘stimulated autonomous sys- } I 
tems,’’ as contrasted with “isolated autonomous systems’ 
there is no driving term 
his being a matter of definitions, the issue, of course, cannot 
be decided logically However, new definitions or usages shoul 
not be introduced without this fact being mentioned openl 
and without giving excellent reasons for the departure fro! 
generally accepted terminology } 
Even judging Clauser’s terminology on its own merits, 
” 


would have to raise objections to calling ‘isolated’’ a system like 


the one described by van der Pol’s differential equation 


because there is a continuous exchange of energy; mechani 
energy has to be fed into the system repe ited] period ill 
In steady state Hence, energ\ wise at least, the svstem cer 
tainly is not “isolated.” 


3) The writer maintains that a departure from recently intt 


duced terminology is necessary in order to comply with a muc 


older and very well established one 
Everybody, engineer or not, calls a spring hard if its spring 
constant (slope of the force-deflection curve) is large and soft 1 


the spring constant is small. In 19388 M. Rauscher? started a us 


age where a spring is called hard if the slope of the force deflec 
tion curve is increasing and soft if that slope is decreasing Phis 
usage, unfortunately, has spread, and it has been adopted b 
Clauser, although it conflicts with the much older terminology 


and hence is bound to create confusion The writer wants t 





repeat here his earlier suggestions** of resolving the diffict 
by retaining the words “hard”? and “soft’’ for describing the 
slope (of the foree-deflection curve) and by using the words This 
the change in sl ype 


“hardening” and ‘‘softening’’ for describing 


unplitu 
illed 

the aut 
fore H 
placement 
TY inagnet! 


plane 


ned 
lerg\ 
»-frequen 


In an 


dissi 


rilv lead t 


show 

erm 

the ju 

ot help t 

and ‘‘n 

isage. He 
} 


e agoes tl 
f the 


differenti 


inde 





uently, he 
ye varying 
1 the right 
NOUS SYS 
ystems” if 
se, cannot 
ges shoul 
ed openl 
ture fror 
lerits, one 


vstem like 
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The Shear Flow Along a Flat Plate With Uniform 
Suction 


Department of Physics, Faculty of Science, University of Tokyc 


kyo, Japan 


noer 5, 1956 


R: CENTLY, several authors have investigated the boundary 

laver in a shear flow In this note, an exact solution of the 
Navier-Stokes equations will be presented, which represents the 
youndary layer along an infinite flat plate with uniform suction 


situated in a shear flow 
Let 


to the plate, respective ly 


ixis and y-axis be taken parallel and perpendicular 
If the flow 


the x 


condition is assumed to 


steady and independent of the x-coordinate, then the Navier 
Stokes equations and the equation of continuity for incompressible 
fluid become 
au dy — (1/p) (Op/Ox + vid-u/ay" l 
dv dy = —(1/p Op/oyv + r(d dy 2 
dv/dy = 0 3 
where (1, ire the velocity components, p is the pressure, p is 
the density, and v is the kinematic viscosity 
Eq. (3) gives 
const 7 } 
Hence, Eqs 1) and (2) become 
Op/Ox p|v(d?u/dy*) — du/dy 5) 
Op/ov = 0 6 
respectively From Eq. (6), p is found to be a function of 4 
ly, while Eq. (5) shows that 0f/O0x is a function of y only 
rherefore 0p/Ox should be a constant (a, say Thus, 
a const 7 
nd Eq. (5) reduces to 
v(d2u /dy? au ay a/p bad 


The solution of Eq. (8), subject to the appropriate boundary 
ndition, 
{ 0 at 1 0 Qg 
easily found to be 
t wow v+tl — 1. exp v/s 10 
ere 1 is an arbitrary constant and 
w —a/p 1] 
If it is assumed that 
0) (12 
that there is a uniform suction along the plate, then the last 
erm of Eq. (10) shows the exponential decay as y increases 
| 
it 


a> ony + I = 13 


as v 


indicates that the above obtained solution corresponds to 


he boundary layer over a flat plate in a shear flow defined by 
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the thickness of the boundary laver being of order of 0 

It is interesting to note that the solution [E« LO)! « ) 
constructed by sing the shear flow d the usu 
suction flow®; this is because of the li rit f the basic equ 
tions, Eqs. (1), (2), and (3 \lso, it is to be noted that the pre 
sure gradient a, the vorticit , and the sucti velocity re 
related to each other by Eq. (11 \ccording if the usu 
boundary condition, u —~ | sy=> instead of Eq. (12) i 


issumed, then the solution reduces to the 


no pressure gradient 


Finally, it should be emphasized that the boundar iver 
nature appears only in the case of suctior 0. Indeed, if 
0, the solution [Eq. (10)] represents an exponentially di 
verging velocity as y — o, and if O—-1e., no suction is 
sumed—Eq. (10), subject to the boundary condition of Eq. (9 
has only the plane Poiseuille flow as its solution 
u = (a/2pyv)y* + By, 8 const 
In conclusion, the writer would like to express his cordial thank 
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buckling of a 
is vielded useful 


rb 


originally derived for 


ONNELL’S 


EQUATION, 


il loading! h 


D 


results in a large variety 


circular cvlinder under torsion 


of load applications for cylinder pré 


lems. This note contains a discussion of an equation in the 
Donnell form which has application to thin shells of constant 
thickness his equation is applicable to plates and cylinders 


over the entire surface and may be applied, with certain limita 


tions, to more general shells. The equation has been used herein 


to determine the classical expression for the buckling stress of a 
sphere under uniform external pressur¢ 
In 


are 


this analysis, orthogonal curvilinear coordinates, x and vy, 


used in the manner described by Reissner.2 This permits 


writing the metric relation for a thin shell in the form 


where ad 


The 


lations are 


stress-strain-displacement 
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Xx 


Phe ensuing development will be restricted to surfaces with 


constant but not necessarily equal radii of curvature for which 


Oa) / VE 
OJa»/ O§ Oa 


Utilizing Eqs. (2), this permits writing the first th 


in the form 


(Ou/ OX 

Ov/Ov) + (a 

(Ou/OV) + (Ov/OX t/G 
When these are employed in deriving the compatibility equation 
for the stresses in a small plate-like element at the origin, is 


found that (using the Airy stress function, F 
1/E)V4*F = (0?/dx?) (w/R,) + (0 
If Donnell’s simplification for cylinders, 
R 


is extended to shells, 


u/R,< 1/a;) (OW/OE8), 7 <. | 2) (Ow / OE 7 
then the curvature relations of Eqs. (3) will lead to the following 
equilibrium equation for transverse forces acting on the plate 
like element: 

DV tw + tlozr (0?w/Ox?) + 27(07w/Oxdy) + o,(07?w/O: 
(1/R,) (0?F/oy?) + (1/R,) (0? F/ox?)] +p =0 (8 
A form of Donnell’s equation may be obtained by eliminating 
F between Eqs. (6) and (8 
DV8w + tV4[o.(07w/Ox?) + 27(02w/Oxdy) 4 j(0?w / Oy? 
Vip + (Et/R,R,) }R,[(04/0,1) (w/R 


/ 
(0*/Ox?0y?) (w/R R,[(04/Ox*) (« 
(0*/Ox?0y?) (w/R 


u/aiR 


This equation is applicable to plates and cylinders. For other 
surfaces, such as a sphere, it is restricted to a small region in the 
neighborhood of the origin of coordin nsequently, it 
may be limited for purposes of stress analysi However, for 
determination of buckling stresses, which illy involve small 


regions of a shell, it could have useful application 


When R, is chosen at a constant value, R, and RX, is equal te 
i 


aR (@ = arbitrary constant), then Eq. (9) reduces 


DV'w + t¥4 [o Ox?) + 27 (0?w/Oxdy) + a, (07% 


Vip + (Et 0?/Ox?) + (1/a) (02/04 w) = 10) 


When a = , Eq. (9) reduces to Donnell’s equation for a 


‘vlinder, and, when a = 1, the expression for a sphere is obtained: 


vil DV 


As an illustration of the simplicity gained through the use of 
Eq 11 , it is employed to find the classical buckling stress fora 
thin-walled spherical shell under uniform external pressure, 
The solution to Eq. (11) is assumed to be 

= Wy Sin mx sin my 12) 


in which the buckle wave lengths are required to be small relative 
toR. Since o, = c, = ¢ (= pR/2t), it is found that 
0 = 2(m?D/t) 4+ t/(2m?R? 13) 


Minimization with respect to m leads to the well-kn 
result 
- (Et/R 
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. . «» Response of Simple Panels to 
Intense Acoustic Loading 


(Continued f 


laboratory siren studies have shown the fatigue life 
to be markedly shortened for these panels by a nominal 


increase in the noise intensity level. 
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